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MASS-FLOW RATES AND TOTAL- PRESSURE RECOVERIES 
FOR A SEMICIRCULAR OPEN-NOSE INLET 
AT MACH 6 

By John S. Dennard and Lillian R. Boney 
Langley Research Center 

SUMMARY 

An all- internal- compression inlet has been tested in a Mach 6 tunnel. The inlet 
was semicircular with provision for radiation cooling over much of the compression 
region; boundary-layer bleed was incorporated slightly upstream of the throat. A 
FORTRAN IV computer program employing the appropriate shock and characteristic 
equations was used to determine the inlet internal aerodynamic flow field. Pressure 
distributions in the supersonic diffuser were in good agreement with those predicted by 
the method of characteristics. The maximum total pressure in the subsonic diffuser was 
18.7 percent of the free-stream total pressure, with a corresponding bleed flow rate of 
11.3 percent. Maximum kinetic -energy efficiency was 91.5 percent, a level which is in 
the range of interest of practical engine applications. 

INTRODUCTION 

The use of air-breathing engines in the hypersonic flow regime is being seriously 
considered because of the high theoretical specific impulse. Only chemical rockets have 
previously been used in the hypersonic flow regime. The inlet which supplies air to a 
hypersonic air-breathing engine will be required to operate at very high altitudes and 
velocities. Temperatures will reach very high values on leading edges and in regions 
such as the throat where the cross-sectional area is restricted and high pressures occur. 
Thus, the conventional requirements of high-pressure recovery and low drag are com- 
pounded by the cooling requirements in hypersonic atmospheric flight. 

Two-dimensional all-internal-compression inlets have been tested and reported for 
Mach numbers from 1 to 6 in reference 1, for Mach numbers from 1.6 to 4.1 in refer- 
ence 2, and for Mach number 6.9 in reference 3. In reference 4 is presented results of 
a three-dimensional axially symmetric all-internal-compression inlet tested at Mach 
number 5. Boundary- layer bleed flow rates for these investigations range between 5 and 
30 percent of the total inlet mass flow. Results of reference 1 indicate that the 


supersonic boundary- layer growth and large wetted areas are major problems in 
achieving high-pressure recovery in two-dimensional inlets and that the large wetted 
areas also present a major cooling problem for the internal surfaces of the inlet. 

The internal flow performance of a three-dimensional semicircular inlet designed 
for radiation cooling throughout most of the compression region is discussed in the pres- 
ent paper. The inlet has essentially zero wave drag, and compared with other inlets, has 
a small wetted-area throat which is suitable for a simple wall- cooling arrangement. It 
can be incorporated in an engine pod or integrated with the airframe. A computer pro- 
gram written in FORTRAN IV language to determine the supersonic rotational flow field 
by the three-dimensional method of characteristics with the appropriate shock equations 
was used in designing the inlet. The computer program is written in two parts and is 
presented as an appendix. Data included are longitudinal static-pressure profiles 
throughout the inlet and subsonic diffuser, total-pressure recovery in the subsonic dif- 
fuser, and mass-flow ratios. 


SYMBOLS 


Moo free- stream Mach number 

w mass-flow rate, kilograms/second (measured by calibrated plug-venturi) 

Woo free- stream mass-flow rate in a cross-sectional area equal to frontal area 

of inlet lips, kilograms/ second 

p static pressure, newtons/meter^ 

total pressure, newtons/meter^ 
average total pressure, newtons/meter^ 
free-stream total pressure, newtons/meter^ 
radius, centimeters 
X distance from inlet lip, meters 


Pt 

Pt 

R 


2 



y 


ratio of specific heats 




kinetic- energy efficiency, 



7-1 

2 


7-1 

i Pt ) J 

Moo^ 


APPARATUS AND MODELS 
Apparatus 

The tests were run in the Langley 20-inch Mach 6 tunnel, which is described in 
reference 5. The free-stream total temperature is 455° K (820° R) and the total pres- 
sure is 2655 kilonewtons per square meter. The test Reynolds number per meter is 
24.58 X 10®. The unit Reynolds number per meter at Mach 6 and 24 384 meters altitude 
is 5.5 X 10®. Thus, an inlet 0.47 meter in diameter at 24 384 meters altitude would have 
a Reynolds number equivalent to the test Reynolds number. Photographs of the basic 
model are presented in figures 1 and 2, and a sketch is presented in figure 3(a). The 
model features all-internal compression, radiation cooling, boundary-layer bleed, and a 
subsonic diffuser having an area ratio of 3.22 and an equivalent half-angle of 3°. The 
aerodynamic design was chosen after a study of several possible inviscid designs was 
made. The supersonic diffuser is conical and has a half- angle of 9°. This design angle 
was chosen as a reasonable compromise between an isentropic inlet, which has excessive 
length coupled with ideal pressure recovery, and a steeper initial angle, which provides 
a desirable structural form but has poor aerodynamic characteristics. The estimated 
boundary- layer growth effectively increased the inviscid 9° half-angle to 9.126°; this 
value was used as a basis for a final set of inviscid calculations. The inlet was designed 
by using a computer program written in FORTRAN IV language to compute the super- 
sonic rotational flow field by the three-dimensional method of characteristics. The 
FORTRAN program is an essential feature of the investigations, inasmuch as accurate 
and rapid calculations of such features as the curved shocks and the inviscid total- 
pressure recovery are not easily made by hand. Therefore, for the use of other investi- 
gators, the program is included herein as an appendix. The computer program indicates 
that the shock from the inlet lip intersects the model center line at station 14.72 centi- 
meters, a point slightly upstream of the apex of a large notch in the flat side of the inlet. 
Starting the supersonic flow through the inlet is accomplished by opening a flap built into 
the flat lower surface of the model upstream of the boundary-layer bleed. Figure 2 is a 
photograph of the model with the starting flap opened. 
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The conical compression surface of the basic model is split to form a semicircular 
inlet and a flat plate is secured to the open side. The flat plate is notched in a symmet- 
rical V-shape for a distance of 14.801 centimeters from the leading edge to provide an 
opening for radiation cooling of a large portion of the compression surface. (See fig. 1.) 
The notch does not affect the theoretical three-dimensional flow except at the apex, where 
the theoretical internal shock becomes a normal shock; a small amount of spillage is 
expected to occur at that point. 

At station 20.927 centimeters, a boundary- layer-bleed scoop approximately 
0.061 centimeter high surrounds the internal flow passage and the bleed air exits through 
the flat-plate side of the model. The Mach number at the edge of the boundary layer is 
estimated to be 3.10 at this station. Mass-flow integrations of the boundary- layer flow 
indicate that the bleed mass flow is 8.7 percent of the total mass flow. 

The flow undergoes an additional conical compression downstream of the bleed sta- 
tion to the throat, where the Mach number is reduced to 2.95 and the total-pressure ratio 
is 0.727. The throat length is 5.31 times the initial throat radius and the throat expands 
at an equivalent conical half-angle of 0.5° to allow for boundary- layer growth. A normal 
(or terminal) shock at the throat should reduce the Mach number to 0.478 and the total- 
pressure ratio to 0.249. 

The subsonic diffuser further reduces the Mach number to 0.115. At this low Mach 
number, the static pressure and the total pressure have almost the same value; 0.247 and 
0.249 times the free-stream total pressure, respectively. The static-pressure rise is 
113.31 times the free-stream static pressure. A transition from a semicircular to a 
circular duct is located in the duct downstream of the diffuser, and a combination plug- 
venturi is provided for mass-flow measurement and shock-position control. 

Static pressures are measured at orifices located in a longitudinal row on the model 
from the inlet lip to the diffuser-measuring station. The diffuser is considered to be 
divided by concentric semicircles into 12 sections equal in area. One tube is placed in 
each of the 12 sections midway between the semicircles to measure the total-pressure 
recovery. A schlieren system is employed to determine visually when starting occurs 
and to see the flow disturbances occurring at the bleed flow exits. 

The area distribution within the model is presented in figure 4. The capture area 
is 43.89 square centimeters and the length of the initial supersonic compression region 
(fig. 3(a)) is 20.927 centimeters, with a downstream area of 6.11 square centimeters at 
the boundary-layer-bleed station. 

The minimum throat area is 4,19 square centimeters at station 26,327 centimeters; 
the area increases gradually to 4.77 square centimeters at station 35.334 centimeters. 
Downstream of this point, the area increases to 15.35 square centimeters in the subsonic 
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diffuser. There are no struts or steps in the internal flow passages other than the 
boundary- layer- bleed slot; changes in the slope of the area curve are all moderate. 

Models 

A few preliminary runs of the basic configuration (fig. 3(a)) showed that the model 
throat was choked, that a shock was standing in the supersonic diffuser, and that air was 
spilling from the apex of the V-shaped notch. In order to increase the bleed mass flow 
and relieve the choking, the basic configuration was modified several different ways to 
make four new configurations. 

The modifications which produced configuration I involved three changes of the 
basic configuration: the flap was shortened by 1.27 centimeters; the approach to the 
boimdary- layer-bleed slot was cut away at a 15° angle starting at station 19.812 centi- 
meters; and the size of the bleed exits was doubled and two new bleed exits 2.54 centi- 
meters long were made adjacent to the opening at the flap trailing edge. A sketch of con- 
figuration I is shown in figure 3(b). 

Configuration II was formed by rounding the approach to the boundary- layer-bleed 
slot to a radius of 5.08 millimeters and thus reducing the effective size of the bleed flow 
slot. This modification, which was accomplished by placing a new insert at the entrance 
to the bleed slot, was made in an attempt to decrease the strength of the expansion 
entering the throat. At the same time, the downstream edge of the boundary- layer-bleed 
slot was notched to break up any shock concentrations within the throat. This configura- 
tion is shown in figure 3(c). 

Two additional modifications were made to reduce the boundary- layer thickness 
and the extension of the normal shock into the subsonic diffuser. The bleed exits were 
further enlarged by cutting them 1.91 centimeters farther forward and perforating the 
throat with 10 holes 1.59 millimeters in diameter. A sketch of this model, designated 
configuration in, is shown in figure 3(d). 

Configuration IV was made by closing the throat bleed holes and doubling the size 
of the serrations at the downstream ecfee of the boundary- layer-bleed slot. A sketch is 
shown in figure 3(e). 

The semicircular model used in the present investigation is not suitable for simple 
variable-geometry configurations. Therefore, the tests were made with all the compon- 
ent parts except the starting flap and the exit plug in a fixed position. 
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RESULTS AND DISCUSSION 


Static-Pressure Distributions 

The theoretical static-pressure distribution is presented in figure 5(a), and the 
modified experimental pressure distributions for configuration I are presented in fig- 
ure 5(b). Each curve represents a different setting of the area at the plug-venturi. 
Upstream of the bleed station, the predicted and measured pressures are in good agree- 
ment. Downstream of the bleed station, however, the measured pressures indicate a 
pressure drop which probably originates as an expansion on the cutaway portion of the 
flap and the cutaway approach to the boundary-layer-bleed slot. This expansion results 
in high Mach numbers and, consequently, in large reductions of total pressure in the ter- 
minal shock that follows. This condition should be avoided whenever possible by using 
scoops rather than flush slots for boundary-layer removal. The terminal- shock pres- 
sure rise is distributed throughout the throat and much of the subsonic diffuser rather 
than concentrated in the throat as it should be for maximum recovery. The terminal 
shock is followed by subsonic compression to the station at which the total pressure is 
measured. At the condition at which the static -pressure rise is a maximum, the bleed 
mass-flow rate is 13.3 percent of the capture mass flow. The values of the average total 
pressure measured at the end of the subsonic diffuser and the mass-flow ratio as meas- 
ured by the plug- venturi are shown in the right-hand margin of the figure. In each case 
for comparison, the average-total-pressure ratio is also plotted by using the same sym- 
bol that was used for the corresponding static-pressure-ratio curve. Static -pressure 
measurements for configuration II are presented in figure 5(c). The results are similar 
to those shown in figure 5(b). 

Static -pressure measurements for configuration III are presented in figure 5(d). 

The static -pressure distributions reveal the strong expansion immediately downstream 
of the boundary- layer-bleed slot noted previously in figures 5(b) and 5(c) and terminal- 
shock and subsonic -diffuser characteristics similar to those of figures 5(b) and 5(c). The 
maximum static pressure is about 1 percent of p. ^ greater than the level obtained for 
any previous test condition. This maximum pressure occurs when the terminal shock 
originates in the converging duct upstream of the throat. Stable operation at this condi- 
tion suggests that most of the boundary layer has been removed by the enlarged bleed 
exits. The bleed flow rate, defined by 1 - -^, is 2.7 percent greater at the point of 
maximum total-pressure recovery than the corresponding rate for configuration II. 

Static -pressure distributions for configuration IV are presented in figure 5(e). The 
expansion downstream of the boundary-layer bleed is still apparent, but a distributed 
pressure rise occurs in the throat region for all flow conditions. The terminal shock is 
still extended downstream into the subsonic diffuser, but the maximum, diffuser static 
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pressure has increased to 0,183 times the free-stream total pressure and the corre- 
sponding bleed mass-flow rate is 4.4 percent below the rate for configuration III. The 
terminal shock for maximum total-pressure recovery starts slightly upstream of the 
throat section but appears to be quite stable in this configuration. This stability indi- 
cates that the boundary layer has been effectively removed. 

Total- Pressure Ratios 

Total-pressure- ratio distributions within the constant-area duct downstream of the 
subsonic diffuser are presented in figure 6. Lines of constant-total-pressure ratio are 
mapped at intervals of 0.0010, The more constant-total-pressure distributions occur for 
the critical and subcritical operating conditions (i.e., the terminal shock at the throat or 
upstream of the throat, respectively), while large total-pressure gradients exist for 
extreme supercritical operating conditions (i.e,, terminal shock downstream of the 
throat). These large total-pressure gradients resulted because the normal shock was in 
a region of high velocities where thick boundary layers with lambda-type shocks would 
be likely to exist. These shocks extend into the constant- area duct downstream of the 
subsonic diffuser. The results of this shock extension may be seen in the static-pressure 
distributions presented in figure 5. 

The variation of total-pressure ratio with mass-flow ratio is presented in figure 7. 
These data are similar for all four configurations in that none of them have any useful 
subcritical flow regions for the normal, closed-flap mode of operation. This character- 
istic is typical of internal-compression inlets. Configuration IV with large serrations 
on the downstream side of the bleed flow slots provides the peak total-pressure recovery 
of 0.187 (design value is 0.249) at a value of the mass-flow ratio of 0.887. This mass- 
flow ratio corresponds to a bleed mass-flow rate of 11.3 percent, which is 2,6 percent 
larger than the design value of 8.7 percent because of the enlarged bleed flow entrance 
and exit areas. Although the total-pressure recovery and mass-flow ratio drop sharply 
for subcritical operation with the flap closed, opening the flap slightly permits excess 
mass flow to bypass the throat and the dashed-line portion of the curve results. Opening 
the flap improves the total-pressure recovery because the position of the normal shock 
can be located in the region near the boundary-layer-bleed slot rather than the higher 
Mach number region near the apex of the V-shaped notch. The results obtained from 
configuration IV indicate that the mass-flow rate may be reduced to a level of about 
60 percent, while the total-pressure ratio is reduced only about 1.5 percent. Similar 
stable subcritical-mass-flow operation is reported in reference 6 for a fixed-geometry 
asymmetric inlet operating at Mach number 3,85, The inlet of reference 6 at several 
angles of attack also provided mass-flow and pressure- recovery performance character- 
istics superior to those of axially symmetric inlets. It is believed that the inlet of the 
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present investigation would have similar performance characteristics at a nonzero angle 
of attack because both inlets have increasing frontal areas as the angle of attack is 
increased. 

Although the inlet of the present investigation was designed with a subsonic dif- 
fuser, similar inlets may be used to supply air to supersonic combustion ramjets. Sur- 
veys were attempted at the throat to determine the recovery and to compare the experi- 
mental and theoretical flow conditions. In every attempt, the rake choked the throat and 
endeavors to start the flow destroyed the rake. Despite these discouraging results, the 
excellent correlation of static pressures in the supersonic diffuser indicates that the flow 
up to the throat behaved as theoretically predicted. The distributed static-pressure rise 
in the subsonic diffuser probably is a primary cause of the discrepancy between predicted 
and experimental total pressures. Thus, this type of inlet would probably be satisfactory 
for supplying air to a supersonic combustion ramjet. 

Kinetic- Energy Efficiency 

The pressure recovery data of figure 7 are presented in figure 8 in the form of the 
variation of kinetic -energy efficiency with mass-flow ratio. The maximum value of 7]y^ 
is 0.915. This value is 0.017 less than the design value of 0.932. Reducing the mass flow 
by opening the flap and spilling 30 percent of the flow reduces the value of 77 ^^ to 0.909, 
a reduction of only 0.006. These experimental values of kinetic -energy efficiency are 
considered to be within the range of interest for practical engine applications. 

Design-point recoveries and efficiencies were not obtained because expansion waves 
generated at the boundary-layer-bleed slot increased the local Mach numbers and shock 
losses at the throat. 


SUMMARY OF RESULTS 

An investigation of a semicircular all-internal-compression inlet has been made at 
a Mach number of 6 . The inlet was designed by using a computer program written in 
FORTRAN language; the program employed the three-dimensional method of character- 
istics. The investigation indicated the following results: 

1. The experimental pressures in the supersonic diffuser were in good agreement 
with those predicted by the method of characteristics. 

2. The terminal shock was distributed throughout the throat region and much of the 
subsonic diffuser. 
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3. The terminal shock originated in the converging duct upstream of the throat. 

The stability of this flow structure indicated that the boundary layer had been effectively 
removed. 

4. The maximum total-pressure recovery of 0.187, which was 0.062 less than the 
design value of 0.249, occurred with a boundary-layer bleed mass-flow rate of 

11.3 percent. 

5. Subcritical mass-flow rates of 60 percent were achieved with a decrease of only 
1.5 percent in total-pressure recovery. 

6. The maximum value of kinetic- energy efficiency was 0.915, a value which was 
0.017 less than the design value. This value is considered to be in the range of interest 
for practical engine applications. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 30, 1965. 
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APPENDIX 


PROGRAM FOR SUPERSONIC DIFFUSER FLOW FIELD 


A computer program has been written in FORTRAN IV language (ref. 7) in two 
parts to determine the flow field within the supersonic diffuser. The three-dimensional 
method of characteristics with the appropriate shock equations is used. The application 
of the method of characteristics to supersonic rotational flow is described in reference 8. 
The upstream flow field is assumed to be uniform with no Mach number or pressure 
gradients and the direction of flow is assumed to be parallel to the body axis of symmetry. 


The first part of the program defines the initial shock and the flow field between 
the initial shock and the reflected shock, as shown in figure A-1. The input data 
required are the Mach number, the entropy, the equation defining the wall, the flow deflec- 
tion angle on the shock, and the flow deflection angle on the wall. Calculations are made 
at the leading edge and along characteristic lines starting at the shock and going upward 
to the right until they reach the wall. As the flow field becomes filled and the shock 
approaches the center line, characteristic lines of the same family begin to cross each 
other at some distance from the shock. The data obtained in this region, of course, are 
not usable, but the calculations are made in the region near the shock until a subsonic 
Mach number downstream of the shock is reached. If the shock is too steep, subsonic 

flow may exist at large radial distances 
from the center line and the approxima- 
tion for the second part of the program 
will not be accurate. For the present 
case, subsonic flow exists only within a 
radius of 0.0077 times the lip radius and 
thus affects only 0.0059 percent of the 
mass flow. Therefore, the effect is con- 

, , sidered negligible for these calculations, 

(a) Initial shock and flow region which follows. 




(b) Reflected shock and flow region which follows. 

Figure A-1.- Sketch of initial and reflected shocks and associated flow regions. 


The second part of the program 
defines the reflected shock and the super- 
sonic flow field which follows. The posi- 
tion of the reflected shock is first esti- 
mated as follows: 

1. The first point of the reflected 
shock is assumed to be the last computed 
point (nearest the center line) of the 
initial shock. 
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2. The flow direction downstream of the reflected shock is assumed to be parallel 
to the center line. 

3. The flow at the center line ahead of the reflected shock is assumed to be at the 
same Mach number, static pressure, and total pressure as the flow at the last computed 
point of the initial shock. The region of subsonic flow at the center line downstream of 
the initial shock is ignored. 

The path of the reflected shock is estimated on the basis of these assumptions, and equa- 
tions are written for the variation of local Mach number M, flow angle 6, and entropy S 
with distance from the center line. A term is included in these equations to correct for 
slight upstream or downstream deviations of the reflected shock from the predicted path. 
The equations used in subroutine subr are as follows: 

M = 2.685947 + 0.725327Y + (-17.506468y2 + 7.870908Y - 0.152477)^/^ 


- (1.6 - 2.857Y) 


X - (Y + 0.7479) 
0.277 


_ -0.0325 
^radians yO-6475 


1.65Y^*'^^® + 0.0424 


X - (Y + 0.7479) 
0.277 

y1.068 


o _ 153 

yO.675 

where X is the distance downstream of the inlet lip and Y is the radius expressed in 
terms of the lip radius. Calculations are made at the intersection of the initial shock 
and the center line and along characteristic lines starting at the reflected shock and 
going downward to the right until they reach the center line. The calculations are con- 
tinued vuitil characteristic lines intersect the center line at a given distance downstream 
of the shock. 

The printout data of the first part of the program included the following parameters 
for each point calculated in the flow field: 

X distance downstream of inlet lip, expressed in terms of lip radius 

Y radius to point, expressed in terms of lip radius 

e flow angle, degrees 

M local Mach number 
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W ratio of local to limiting velocity 


S = R loggj 



(where R is the gas constant) 


R. 

Pi 


ratio of local static to free- stream static pressure 


— i static-pressure rise divided by free-stream dynamic pressure 

The printout data of the second part of the program differ from those of the first 
part in that total-pressure and mass-flow parameters are included and static-pressure 
ratios are omitted. This change is made to facilitate a mass-weighted integration of the 
flow field at the station where the reflected shock intersects the inlet wall. The total- 
pressure and mass-flow parameters are as follows: 


Pt 

Pt,l 


ratio of local total pressure to free-stream total pressure 




parameter used to determine mass-flow rate 



parameter used to determine mass-weighted total-pressure 
ratio 


The computer program for the IBM 7094 electronic data processing system with 
typical input data follows. 
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SAMPLE INPUT TO INITIAL SHOCK REGION 

SOATA 

J1 =0 « J2 = 9999. . J3=I .4 , J4 = l 716. « J5=l . . J6 = -. 1 606395, J7 = 6« , J8 = 0, J9= I . » 
J10=.035, J1 1 =.9943776, 312=4.7645, J13=-9. 1260 145, 314=293.013, 

315 = - 16. 7598 18, 3 16=1 .3. 

WHERE 

31 Y LIMIT ON SHOCK 

TYPICALLY=0 TO STOP CALCULATION WHEN INITIAL SHOCK 
REACHES THE CENTERLINE 

32 X LIMIT ON WALL 

TYPICALLY=9999. TO ALLOW Y LIMIT TO STOP THE CALCULATION 

33 GAMMA 

TYPI CALLY=1 .4 FOR AIR 

34 R GAS CONSTANT 

TYPI CALLY=1 71 6. FOR AIR 

35 B WHERE Y=Ax+8 IS EQUATION DEF I NI NG ■ WALL 
TYPICALLY=1 . 

36 A WHERE Y=Ax+B IS EQUATION DEFINING WALL 
TYP I CALL Y = -. 1606395 


37 MACH NUMBER AHEAD OF INITIAL SHOCK 
TYPICALLY=6. 


J8 

XB COORDINATE 
TYP I CALLY=0 

OF 

FIRST 

POINT 

ON 

WALL 

AND 

INITIAL 

SHOCK 

J9 

YB coordinate 

OF 

FIRST 

POINT 

ON 

WALL 

AND 

INITIAL 

SHOCK 


TYPICALLY=1 . 

310 XA coordinate OF SECOND POINT ON WALL 
TYPICALLY=.035 

311 YA COORDINATE OF SECOND POINT ON WALL 
TYPI CALLY= .9943776 

312 MACH number AT FIRST AND SECOND POINTS 
TYPICALLY=4. 7645 

313 FLOW DEFLECTION ANGLE AT FIRST AND SECOND POINTS — NEGATIVE 

angle degrees 

TYPI CALL Y= -9 .1260145 

314 ENTROPY AT FIRST AND SECOND POINTS 
TYPICALLY=293.013 

315 SHOCK ANGLE AT FIRST POINT NEGATIVE ANGLE DEGREES 

TYPICALLY=-16.759818 

316 +1. TO PRINT ALL POINTS, OR =0 TO PRINT ONLY WALL AND 
SHOCK POINTS 

TYPICALLY=1 . 

$ INDICATES END OF INPUT DATA 


13 



APPENDIX 


subroutine bodyx computes the wall equation y=ax+b. 

A NEW subroutine WILL BE NEEDED IF THE EQUATION OF THE 
WALL IS CHANGED. 

SAMPLE INPUT TO REFLECTED SHOCK REGION 

SDATA 

j j =9999. ,J2 = 5..J3=1.4.J4=1716. .08 = 2.785644 .J9= ,022561 ,010 =2.785644, 
011=0,012=1.1 538616 ,013=0,014=3711 .433,01 5=29.610459,016=1 . S 
WHERE 

01 Y LIMIT ON SHOCK 

TYPICALLY=9999. TO ALLOW X LIMIT TO STOP THE CALCULATION 

02 X LIMIT ON CENTERLINE 
TYP1CALLY=5. 

03 GAMMA 

TYPICALLY=1 .4 FOR AIR 

04 R GAS CONSTANT 
TYPICAlLY=1716. FOR AIR 

05 NOT USED ON REFLECTED SHOCK 

06 NOT USED ON REFLECTED SHOCK 

07 NOT USED ON REFLECTED SHOCK 

08 Xri coordinate OF FIRST POINT ON REFLECTED St-iOCK 
TYPICALLY=2,7S5644=LAST COMPUTED POINT OF INITIAL SHOCn 

09 Y6 coordinate OF FIRST POINT ON REFLECTED SHOCK 
TYP I C All Y= .02256 I 

01'- XA COORDINATE OF CENTERLINE POINT 
TYP I CALLY=2, 785644 

011 YA coordinate OF CENTERLINE POINT 

typically=o 

012 MACH number at FIRST POINT ON REFLECTED SHOCK AND 
CENTERLINE POINT 

TYP1CAlLY=1 .1338618 

013 FLOW DEFLECTION ANGLE AT F I RST POINT ON REFLECTED SHOCK 
AND CENTERi_INE POINT 

TYPICALLY=0 

014 ENTROPY AT FIRST POINT ON REFLECTED SHOCK AND 
CENTERLINE POINT 

TYPICAlLY=371 1 .433 

015 SHOCK angle AT FIRST POINT ON REFLECTED SHOCK 
— POSITIVE ANGLE DEGREES 

typically =29. 61 0459 

016 +1. TO PRINT All points, OR =0 TO PRINT ONLY 
centerline and shock POINT 

TYP I CALLY= 1 . 
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$ indicates end of input data 

SUBROUTINE SUBR COMPUTES ENE WM » T HC 1 . SC 1 AS FUNCTIONS OF 
XC AND YC. these ARE THE EQUATIONS FOR M,-THETA»S 
AHEAD OF THE REFLECTED SHOCK. A NEW SUBROUTINE 
WILL c)E NEEDED IF EQUATIONS ARE CHANGED. 
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IBFTC main ref 

MA I NPROGRAM-F I RSTSHOCKREG I ON 

COMMONXA t YA « THA . WA , SMUA , CMUA « STHA , CTHA , SA , X3 » YB » THB » 

1 WB « SMU8 « CMUB » STHB « CTH3 , S3, EP3 » TEPB , XC » YC » THC » WC « 

2SMUC . CMUC » STHC » CTHC « SC . EPC « TEPC « EMUA » EMU3 » EMUC . TH AC . 

3WAC . SMUAC » YAC . EMUAC » STHAC , EL AC * TH3C » WoC . SMU3C t Y3C » EMUBC . 

4STHBC • EPBC » TEPBC « EMBC t TANMA » TANMS . EMC * PC t CPC . EM » THCP » 

5C0S1 ,TAN1 ,C0S2,TAN2»K0UNT,TAPE»AL1GHT»3EGS»SVC0N,SVC0N1 » 
6SVCON2.SVCON3tSVCON4,SVCON5«SVCON6iSVCON7tSVCONB,SVCON9,SVCOlO» 
7SVC01 1 , EEP.EEPl « EEP2 , EEP3 t EEP4 « EEP5 t EEP6 t EEP7 , C NTRL , CNTRL 1 . 
8CNTRL2 » CNTRL3 . CNTRL4 .CNTRLB . CNTRL6 » KNTRL7 . CNTRuS , CNTRL9 , CNTR 1 0 * 
9KNTR 1 1 « CNTR 1 2 « CNTRl 3 i CNTR 1 4 » CNTR 1 5 . CNTR 1 6 » CNTR 1 7 « CNTR 1 8 t CNTR 1 9 » 
XCNTR54« ERASE* ERAS 1 , ERAS2 * ERA S3 * ERAS4 , ER AS5 , ER A56 * ERAS 7 * 

XERAS8. ERAS9* ERAS 1 0, eras 1 1 * ERAS12 *ERAS1 3 * ERAS 1 4 , ERAS 1 E * ERAS 1 6 
COMMONERAS17.ERAS18*ERAS19*ERAS20 «ERAS2 1 * ERAS22 * ERAS23 * 

1 ERAS24 ,ERAS25 .ERAS26 *ERAS27*ERAS28 , ERAS29 , ER AS30 « ERAS 3 1 * 

2ERAS32 *ERAS33,ERAS34 ,ERAS35 * DEC . DEC 1 * DEC2 , DECS , DEC4 , 
3DEC5*DEC6*DEC7*DEC8.DEC9*DEC10*D£C1 1 *K0N.C0N1 *PRT* 

4PRT1 ,PRT2*PRT3.PRT4,PRT5*PRT6.PRT7,PRT8 *PRT9,DEG* 

5TEMP * TEMPI * TEMP2 * TEMP3 * TEMP4 * X I « Y I , TH I * W I * SMU I « 

6CMUI ,STHI,CTHI,SI .KERl 0 * EX 1 * EX2 *C1 ,C2*C3,C4*C5,C6* 
7C7«C8*C9.C10*C11*C12*C13«C14*C13«C16*C17,C1B,C19*C20* 

8C21 ,C22*C23,C24,C25,C26*C27*C2d*C29*C30*SAv*SAVl * 

9SAV2 * SAV3 * SAV4 ,SAV5 . SA V6 * SA V7 * SAV8 * S A V9 » SC 1 * THC PR * ENEWST » 
XENEWCT .EKP.EKPl , EKP2 , EKP3 , EKP4 « EKP5 * ENEWTH , T HC 1 * ENEWM , 

XAA .KAC.NRET.SW2 *EN (9 ) , PRNOPT * HOLD , J , l< ♦ L * STOR (1000) 

COMMON LTEST.BSHKl 
TAPE=2. 

alight=o 

DEC=0 
DECl = 1 . 

DEC2=2. 

DEC3= 1 .570796327 

DEC4=6. 2831 85303 

DEC5=0 

DEC6=0 

DEC7=7 

DEC8=0 

DEC9=9 

DECl 0=0 

DFCl 1=11 

KON= 1 

CONI =3. 

DEG=57. 2957795 
9 CALLSLITE(O) 
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WRITE(6«98) 

98 FORMAT ( 1 H 1 , 1 HX * 1 5X . 1 HY . 1 5X *9HTHETA DEG . 7X . 1 HM . 1 5X ♦ 1 HW * 1 5X * 1 HS * 1 5X » 

14HP/PI . 12X«9H(P-P1 )/QI ) 

DO 999 NN=1 . 1000 
999 ST0R(NN)=0 
HOLD=l . 

J= 1 

k:=i 

U=1 

CALLSL I TE (2 ) 

KERI 0=0 

CALL. FINP( 16«CNTRL«CNTRL1 . CN TWL3 » CNTWL5 « CNTR 1 5 « CNTR 1 6 , EM 1 « Xd « YB ♦ 

1 XA . Y A « EMB * THB « Sd , EPB , PRNOPT ) 

EPB=EPB/DEG 
THB=THB/DEG 
CNTRJ_6 = 9999. 

KNTRL7=3 
CNTRL8= 1 . 

CNTRL9= 1 . 

<NTR1 1 =0 
EEP5=EM 1 *EM1 
EEPl =EEP5+EEP5 
EEP4=EEP1+DEC1 
EEP=EEP5*EEP5 

EEP3=( (CNTRL3+DEC1 ) /DEC2 ) *EEP5+DEC1 

EEPl =-EEP3*IZEPl 

EEP2 = - ( EEP5+DEC2 ) »-EEPd 

SW2 = 0 

W8=SQRT (DECl/ (DEC1+DEC2/( (CNTRL3-DEC1 )*EM8*EMB) ) ) 

SMUB=DEC 1 /EMB 
SMUA=SMUB 

CMU3=SQRT ( 1 .-SMUB*SMUB ) 

CMUA=CMUB 

STHd=SIN(THB) 

CTHB=COS ( THB ) 

STHA=STHB 
CTHA=CTHB 
SA = SB 
WA = WB 
THA=THB 

TEPB=SI N (EPB )/COS (EPB ) 

BEGS=EPB 
BSHK 1 =0 
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CALLSHKPT 
CALLPRINT 
CALLWRTSK 
BSHKl =1 . 

CALLMOVCB 

CALLMOVAC 

CALLMOVBA 

CALLIPT 

GO TO ( 22, 22 » 22 22. 23 ) ,NRET 

22 CALLMOVCS 
CALLWOV 1 C 
CALLWRTGN 
CALLPRINT 
CALLMOVSC 

23 CALLWRTGN 
CALLPRINT 
CALLSPACE 
CALLWRTFL 
CALLRDB 

GOTO! 10,11 ) ,NRET 
10 CALLEXIT 
1 1 CALLRDA 

GOTO (12,13) ,NRET 

12 CALLEXIT 

13 CALLSHKPT 
CALLPRINT 
CALLWRTSK 
CALLMOVCA 
CALLRDB2 

GOTO( 14,15) ,NRET 

14 CALLEXIT 

15 CALLGEN 
CALLIPT 

GOTO (24, 24, 24, 24, 24, 25) ,NRET 

24 CALLMOVBA 
CALLMOVCB 
CALLMOVIC 
CALLWRTGN 
G0T07757 

25 CALLMOVBA 
CALLMOVCB 

28 CALLWRTGN 
CALLPRINT 
CALLBODY 
CALLWRTGN 
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CALLPRINT 
CALLSPACE 
CALLWRTFL 
7 CALLRDB 

GOTO( 16*17) *NRET 

16 CALLEXIT 

17 CALLRDA 

GOTO (18*19) *NRET 
13 CALLEXIT 

19 CALLSHKPT 
CALLPRINT 
CALLWRTSK 

IF ( YC-CNTRL)5*30* 30 
30 CALLMOVCA 
CALLRDB2 

GOTO (20,21 ) *NRET 

20 CALLEXIT 

21 CALLGEN 
CALL I PT 

GOTO ( 26 * 26 *26 * 26 ,26 * 27 ) ,NRET 

26 CALLMOVCA 
CALLMOV I C 
CALLWRTGN 
CALLPR I NT 
CALLMOVAC 

27 CALL'^'RTGN 
CALLPRINT 

4 CALL(V)OVCA 
CALLRCB2 

GO TO (31*32) ,NRET 

31 GOTOl 

32 CALLGEN 

IF (DEC )361 O, 34 ,3610 
34 CALL WRTGN 

IF (PRNOPT )3,2 *3 
CALL PRINT 
GO TO 4 
1 CALLMOVBA 
CALLMOVCB 
CALLRODY 
CALLWRTGN 
CALLPR I NT 
3610 DEC=0 
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CALLWRTFL 
CALLSPACE ' 

IF(XC-CNTRL1 )6»6«5 
6 IF(SW2)8,7»8 
8 WR I TE (6 » 100) TAPE, ALIGHT 
5 G0T09 

7757 CALLPRINT 
CALLMOVBC 
G0T028 

100 F0RMAT(2E16.8) 

END 

SIBFTC EPSR REF 

SUBROUT I NEEPSR 

COMMONXA , YA , THA , WA , SMUA , CMUA , STHA , CTHA , SA , XB , YB , THB , 

1 WB,SMUB,CMUB,STH6,CTHB,SB,EPB,TEPB,XC,YC,THC,WC, 
2SMUC,CMUC,STHC,CTHC,SC,EPC,T£PC,EMUA,EMUB,EMUC,THAC, 
3WAC,SMUAC,YAC,EMUAC,STHAC,ELAC,THBC, W6C,SMUdC,YciC,EMUBC, 
4STHBC,EPBC«TEPSC ,EMBC, TANMA, TANM8, EMC* PC, CPC , EM, THCP* 

5C0S1 ,TAN1 , C0S2,TAN2,K0UNT, TAPE, alight, BEGS, SVC0N,SVC0N1 , 

6SVC0N2 , SVCON3 , SVC0N4 , SVC0N5 , SVC0N6 , SVC0N7 , SVC0N8 , SVC0N9 , SVCO 1 0 , 
7SVC01 1 , EEP,EEP1 , EEP2 , EEP3 , EEP4 , EEP5 * EEP6 , EEP7 , CNTHL , CNTHL 1 , 
8CNTRL2,CNTRL3,CNTRL4,CNTRL5,CNTRL6,KNTRL7,CNTRLU,CNTRL9,CNTH10, 
9KNTR1 1 ,CNTR12,CNTR13,CNTR14,CNTR15,CNTR16,CNTR17,CNTR18,CNTR19, 
XCNTR54, ERASE *ER AS 1 , ERAS2 , ERA S3 * ERAS4 , ER ASS , ERAS6 , ERAS7 , 
XERAS8,ERAS9,ERAS1 0,ERAS1 1 ,ERAS12 ,ERAS13,ERAS14,ERAS15,ERAS16 
C0MM0NERAS17,ERAS18,ERAS19,ERAS20,ERAS21 , ERAS22 , ERAS23 , 

1 ERAS24 , ERAS25 , ERAS26 , ERAS27 , ERAS28 , ERAS29 , ERAS30 , ERAS31 , 
2ERAS32,£RAS33,ERAS34 ,ERAS35 ,.DEC , DEC 1 , DEC2 , DEC3 , DEC4 , 
3DEC5,DEC6,DEC7,DEC8,DEC9,DEC10*DEC1 1 ,KON,CONl ,PRT, 

4PRT 1 , PRT2 , PRT3 , PRT4 , PRT5 , PRT6 * PRT7 , PRT8 , PRT9 , DEG , 

5TEMP, TEMPI , TEMP2 , TEMP3 , TEMP4 , X I , YI ,THI ,WI ,SMUI , 

6CMUI ,STHI ,CTHI ,SI ,KER10,EX1 , EX2 , C 1 , C2 , C3 , C4 , C5 , C6 , 
7C7,C8,C9,C10,C1 1 , C 1 2 , C 1 3 , C 1 4 , C 1 5 , C 1 6 , C 1 7 , C 1 8 , C 1 9 , C20 * 

8C21 ,C22,C23,C24,C25,C26,C27,C28*C29,C30*SAV«SAV1 , 

9SAV2 , SAV3 , SA V4 , SA V5 , SA V6 , SAV7 , SAV8 , SAV9 , SC 1 , THCPR , ENEWST , 
XENEWCT,EKP,EKP1 , EKP2 , EKP3 , EKP4 , EKP5 , ENEWTH , THC 1 ,ENEWM, 
XAA,KAC,NRET,SW2,EN(9) , PRNOPT , HOLD , J , K , L , STOP ( 1000 ) 
ERASE=ERASE*ERASE 
ERASl =ERASE*EEP+EEP 
ERAS2= (ERASE*EEP1+EEP2 )/CONl 
ERAS13= (EEP3*EEP3*ERASE+EEP4 )/C0Nl 
ERAS6=ERAS1 3*ERAS2+ERAS1 
ERAS9=ERAS2*ERAS2 
ERAS15=ERAS1*ERAS13-ERAS9 
ERASE=ERAS 1 5+ERAS 1 5 
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ERASl 0= (ERAS6/ERASE ) *ERASl +ERA52 

ERAS6=ERAS6*ERAS6 

ERAS7=ERAS 1 3*ERAS1 3+ERAS2 

AA= (SORT (ABS( ( ERAS7 + ERAS7 ) *ERASE + ERAS6 ) )/ERASE; )*ERAS1 
AA = ATAN2 ( AA «ERAS1 0 ) 

IF ( AA )2 .2. 1 
2 AA=AA+DEC4 
1 ERASF=AA/CONl 
ERASE = COS (ERASE ) 

2RASR=SQRT ( -ERAS 1 5 ) 

ERA5E=- ( (ERAS9+ERAS9 ) *ER ASE+ERAS2 )/ERASl 
ERAS2 = SQRT ( DECl -ERASE ) 

ERASl =-SGRT ( ERASE ) 

TEPC==ERAS1 /ERAS2 
EPC = ATAN ( TEPC ) 

RETURN 

END 

SIBFTC BODY REF 

SUBROUT I NEBODY 

COMMONXA » YA , THA » WA , SMUA, CMUA » STHA , CTHA i SA , XB * YB * THB » 

1 WB«SMUB,CMUB.STHB.CTHB,SB,EPD.TEPB.XC*YC»THC»WC» 

2 SMUG « CMUC . STHC , CTHC . SC » EPC » TEPC » EMU A , EMUB « EMUC » THAC« 

3WAC , SMUAC * YAC »EMUAC » ST HAC , EL AC » THDC * WBC » SMUBC » YdC * EMUBC « 

4STHBC « EPBC« TEPBC ,EMBC* TANMA , TANMB »EMC »PC . CPC » EM, THCP» 

5C0S1 »TAN1 ,C0S2,TAN2,K0UNT,TAPE»ALIGHT,BEGS»SVC0N»SVC0N1 , 
6SVCON2 , SVC0N3 , SVC0N4 , SVCONB , SVC0N6 , S VC0N7 , SVC0N8 ♦ SVC0N9, SVCO 1 0 
7SVC01 1 ♦ EEP.EEPl , EEP2 , EEP3 , EEP4 , EEP5 » EEP6 , EEP7 , CNTRL » CNTRL 1 , 
8CNTRL2,CNTRL3,CNTRL4 , CNTRLb , CNTRL6 , KNTRL 7 , CNTRlS , CNTRL 9 , CNTK 1 0 
9KNTR1 1 ,CNTR12,CNTR13 , CNTR 1 4 , CNTR 1 5 , CNTR 1 6 , CNTR 1 7 , CNTR 1 8 , CNTR 1 9 
XCNTR54 .ERASE , ERASl , ERA S2 , ERA S3 , ERAS4 , ERASB , ERAS6 » ERAS7 , 

XERAse, ERAS9, ERAS 1 0, ERASl 1 .ERAS 12 . ERAS 1 3 . ERAS 1 4 . ER AS 1 5 . ERAS 1 6 
COMMONFRASl 7 . ER A S 1 8 . ER AS 1 9 .ERAS2 0 . ERAS2 1 .ERAS 22 .ERAS23. 

1 ERAS24 . ERAS2B . ERAS26 . ERA S27 . ERAS28 . ERAS29 . ERAS30 . ERAS31 . 
2ERAS32 . ERAS33 . ERAS34 .ERA535 . DEC . DEC 1 . DEC2 . DECS . DEC4 . 
3DECB.CEC6.I>EC7.DEC8,.DEC9.DEC10.DEC1 1 .KON.CONl .PRT. 

4PRT1 .PRT2 .PRT3 .PRT4 .PRT5 .PRT6.PRT7 . PRT8 . PRT9 . DEG. 
bTEMP. TEMPI . TEMP2 . TEMP3. TEMP4 . X I . Y 1 . TH I . W I . SMU I . 

6CMUI .STHI .CTHl .SI .KERIO.EXI .EX2.C1 .C2.C3.C4.C5.C6. 

7C7.C8.C9. CIO. Cl 1 .C12.C13.C14.C15.C16.C17.Cl8.C19.C20. 
bC2 1 . C22 . C23 . C24 . C2B . C26 . C27 . C28 . C29 . C30 . SAV . 5AV1 . 

9SAV2 . SAV3.SAV4 . SAVE . SAV6 . S AV7 . 5AV8 . SAV9 . SC 1 . THCPR .ENEWST . 
XENEWCT . EKP .EKPl .E<P2 . EKP3 . EKP4 . EKP5 . ENEWTH . THC 1 .ENEWM. 

XAA .KAC . NRET . SW2 .EN ( 9 ) . PRNOPT .HOLD . J .K .L . STOR (1000) 
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COMMON LTEST 

WBC=WB 

K0UNT=0 

ERAS2=SMUB/CMUB 
EMUB=ATAN (ERAS2 ) 

SC=SA 

9 THCP=THC, 

CALLBODYX 
IFCCNTRL8)2. 1 *2 

2 ERASE=CTHB*CMUB 
ERAS 1 =STHB*SMuB 
ERASE=-ERAS1 +ERASE 
ERAS2=SMUB/CMUB 

ERASE= ( ERAS2*ERAS1 )/ERASE 
YBC= ( YB+YC )/DEC2 
ERASE= ( (XC-XB )/YBC )*ERASE 
G0T03 
1 ERASE=0 

3 CALLSL I TET (2 *LTEST ) 

G0T0(5*4 ) iLTEST 

5 CALLSLITE(2) 

G0T06 

4 AA=0 
G0T07 

6 IE(CNTRL9)8,7,8 

8 AA= ( ( ( (SC-SB )/CNTRL3 )*SMUB ) /CNTRL5 ) *SMUB 

7 ERASE=AA-ERASE 

WC= (- (THB-THC ) *ERAS2 -ERASE )*WBC+WB 

CALLMUSR 

IF (DEC )21 .20*21 

20 EMUB= ( EMUC+EMUB ) /DEC2 
SMU6=S I N (EMUS ) 

CMUB=COS (EMUB ) 

THBC= ( THB + THC )/DEC2 
STHB = S I N (THBC ) 

CTHB=COS ( THBC ) 

WBC= ( WB+WC ) /DEC2 
KOONT= (KOUNT.+ KON ) 

I F (K0UNT-KNTRL7 )9 . 10.10 

10 IF (ABS (THC-THCP )-CNTRL6) 11.11.9 

11 STHC=S I N ( THC ) 

CTHC=COS (THC ) 

21 RETURN 
END 

SIBFTC THASR REF 
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SUBROUT I NETHASR 

COMMONXA»YA,THAiWA«SMUA,CMUA »STHA,CTHA»SA,XB» YB»THB« 

1 WB « SMUa . CMUB » STH8 . CTHa , SB , EPB t TEPB . XC f YC . THC . WC » 

2SMUC « CMUC 1 ST HC , CTHC t SC » EPC * TEPC * EMU A , EMUB « EMUC , THAC « 

3WAC , SMUAC » YAC » EMUAC » STHAC , EL AC » T HBC » WBC » SMU6C « YoC , EMUBC « 

4STH8C «EPBC* TEPBC «£MBC» TANMA t TANMB , EMC » PC » CPC » EM • THCP* 

5COS1 ♦ TANl « C0S2 «TAN2 ,KOUNT « TAPE » AL I GHT » BEGS » S VCON « SVCONl * 

6SVC0N2 , SVC0N3 , SVCON4 . SVC0N5 « SVC0N6 . SVC0N7 » SVC0N8 , SVC0N9 , SVCO 1 0 « 
7SVC011 . EEPfEEPl « EEP2 « EEP3 . EEP4 . EEP5 . EEP6 . EEP7 . CNTRL « CNTRL 1 » 
8CNTRL2. CNTfJL3,CNTRL4 «CNTRL5» CNTRL6<I<NTRL7< CNTRL8 ,CNTRl9,CNTR10» 
9KNTR1 1 <CNTR12<CNTR13<CNTR14<CNTR15<CNTR16<CNTR17<CNTR18<CNTR19< 
XCNTR54 < ERASE < ERAS 1 < ERA S2 < ERA S3 < ERAS4 <ERAS5<ERAS6<ERAS7< 

XERAS8< ERAS9< ERAS 1 0 < ERAS 1 1 < ER AS 1 2 < ERAS 1 3 < ERAS 1 4 < ERAS 1 5 < ERAS 1 6 
COMMONER AS 1 7 < ERAS 1 8 < ERAS 1 9 < ERAS20 < ERAS2 1 < ERAS22 < ERAS23< 

1 ERAS24 < ERAS25 < ERAS26 <ERAS27 < ERAS28 < ERAS29 < ERAS30 < ERAS3 1 < 

2ERAS32 < ERAS33 <ERAS34 <ERAS35 < DEC < DECl < UEC2 < DEC3 < DEC4 < 
3DEC5<DEC6<DEC7<0FC8<DEC9<DEC10<DEC11 <I<ON<CON1 <PRT < 

4PRT1 <PRT2<PRT3<PRT4<PRT5<PRT6<PRT7<PRT8<PRT9<DEG< 

5 TEMP , TEMPI < TEMP2 <TEMP3<TEMP4<XI < Y I < THl < W I < SMU I < 

6CMUI < STH I , CTH I , S I <KER1 O <EX 1 < EX2 <C1 <C2<C3<C4<C5<C6< 
7C7<C8<C9<C10<C1 1 <C12<C13<C14<C15<C16<C17<C18<C19<C20< 

8C2 1 <C22<C23<C24<C25< C26< C27 < C28 < C29< C30 < SAV < 3AV1 < 
9SAV2<SAV3<SAV4<SAV5<SAV6<SAV7<SAV8<SAV9<SC1<THCPR< ENEWST < 
XENEWCT<EKP<E1<P1 <EKP2 < EKP3 < EKP4 < EKP5 < ENEWTH , THC 1 <ENEWM< 

XAA <KAC < NRET < SW2 <EN ( 9 ) < PRNOPT <HOLC < J <K<,_<ST0R(1000) 
ERAS6=EEP5*ERASE 
ERAS7=CNTRL3-DEC1 
ERAS8=£RAS7+DEC2 
ERAS9=SIN(THA) 

ERAS10=C0S(THA) 

ERASl 1 =ERAS9/ERAS1 0 

eras 1 2 = ERAS6-DEC 1 

ERAS 1 2 = EEP5/ERAS 1 2 

ERASl 3=DEC2*CNTRL3 

ERAS 1 3 = ERAS7/ERAS1 3 

eras 1 4 = ERAS2*ERAS2 

ERASl 5=ERAS1 /ERAS2 

ERAS16=ERAS1 *ERAS10 

ERAS16=-ERAS2^<-ERAS9 + ERAS16 

ERAS17=ERAS16*ERAS1 6 

ERAS18=ERAS2*ERAS10 

ERASl 8=ERASl *ERAS9+ERA S 1 8 

ERASl 9=ERAS1 8*ERAS1 8 


23 



APPENDIX 


ERAS20=ERAS1 6/ERASl 8 
ERAS21 =ERAS6-ERAS13 
ERAS2 1 =DEC 1 /ERAS2 1 
ERAS22=DEC1/CNTRL3 

ERAS22= ( ERAS6*ERAS1 3+ERAS22 ) *ERAS6 
ERAS23=-DEC1/ERAS22 

ERAS21=( (ERAS23+ERAS21 ) /ERAS? ) *DEC2*CNTRL5*EEP5*ERASl *ERAS2 

ERAS23- ( (ERAS8/0EC2)*ERAS12-DEC1 )/ERAS14 

ERAS23=-ERAS 1 2*ERAS 1 2*ERAS8*ERASE+ERAS23 

CALLSSR 

CONTINUE 

ERAS22=ERAS1 1*ERAS1 1 
ERAS24=-DEC1 /ERAS22 
ERAS24= (ERAS24-DEC1 )/ERAS23 
ERAS5=ERAS24*ERAS21 

ERAS26= ( (ERAS8/ERAS7)*ERAS15 )/ERAS20 
ERAS26= (ERAS26-DEC1 ) *ERAS1 7+DECl 
ERAS28 = ERAS16*ERAS1 8 , 

ERAS27= IERAS28/ERAS14 )*ERAS24 
ERAS29=ERAS24-DEC 1 

ERAS27= ( ( (ERAS19-ERAS17)*ERAS15*ERAS29+ERAS27)/DEC2)*ERAS8 

ERAS27=-ERAS29*ERAS28*ERAS7+ERAS27 

ERAS35=DEC1/ERAS26 

ERAS3 = SQRT (ERAS35 ) 

ERAS4= ( (-ERAS35*ERAS3 )/ERAS7 )*ERAS27 

RETURN 

END 

SIBFTC IPT REF 

SUBROUTINEIPT 

COMMONXA«YA,THA*WA»SMUA,CMUA ♦ ST HA , CT HA » SA , XB ♦ YB * THB * 

1 WB*SMUB»CMUB»STHB«CTHB,Sa,EPB»TEPB»XCiYC«THC»WC* 
2SMUC.CMUC»STHC»CTHC»SC»EPC»TEPC.EMUA.EMUB«EMUC»THAC» 

3WAC » SMUAC ♦ YAC « EMUAC » STHAC , EL AC * T HBC » WBC ♦ SMUBC » YBC « EMUBC « 

4STHBC ♦ EPBC « TEPBC * EMBC« TANMA » TANMB *EMC . PC * CPC * EM , THCP » 

5C0S1 .TANl ,C0S2»TAN2*K0 UNT,TAPE* alight .BEGS »SVCON.SVCON l . 

6SVC0N2 . SVC0N3. SVC0N4 . SVC0N5 . SVC0N6 . S VC0N7 » S VC0N8 « S VC0N9 » S VCO 1 0 ♦ 
7SVC01 1 . EEP.EEPl . EEP2 . EEP3 . EEP4 . EEP5 . EEP6 . EEP7 . CNTRL . CNTRL 1 « 
8CNTRL2 . CNTRL3. CNTRL4 .CNTRL5. CNTRL6 . KNTRL7 ♦ CNTRL8 ♦ CNTRL9 ♦ CNTR 1 0 . 
9KNTR I 1 . CNTRl 2 . CNTR 1 3 . CNTR 1 4 . CNTR15. CNTR 1 6.CNTR1 7« CNTR 1 8. CNTHl 9» 
XCNTR54. ERASE .ERAS I . ERAS2 . ERA S3 » ERAS4 . ERAS5 . ERAS6 . ERAS? . 

XERAS8. ERAS9 . ERAS 1C . ERAS 1 1 . ERAS 1 2 . ERASl 3 . ERAS 1 4 . ERAS 1 5 . ERAS 1 6 
C0MM0NERASI7.ERAS18.ERAS19.ERAS20.ERAS21 . ERAS22 . ERAS23 . 

1 ERAS24 .ERAS25.ERAS26 .ERAS27.ERAS28. ERAS29.ERAS30 .ERAS31 . 

2ERAS32 . ERAS33 . ERAS34 .CRAS35 . DEC . DEC 1 . DEC2 . DECS . DEC4 . 
3DEC5.DEC6.DEC7.DEC8.DEC9.DEC10.DEC1 1 .KON.CONl .PRT. 
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4PR T 1 « PRT2 . PRT3 » PRT4 , PRT5 » PRT 6 1 PRT7 » PRT8 * PRT9 « DEG f 
5TEMP« TEMPI , TEMP2 . TEMP3 * TEMP4 » X I »YI ,THI «WI » SMU I « 

6CMUI »STHI ,CTHI ,SI «KER10*EX1 «EX2«C1 , C2 « C3 , C4 « C5 . C6 t 
7C7.C8,C9«C1 0«C1 1 i C 1 2 ♦ C 1 3 » C 1 4 . C 1 5 « C 1 6 « C 1 7 » C 1 8 « C 1 9 . C20 ♦ 

8C21 «C22tC23,C24»C25,C26*C27tC28«C29.C30»SAV»SAVl « 
9SAV2»SAV3*SAV4«SAV5»SAV6»SAV7»SAV8.SAV9.SC1 ♦ THCPR « ENEWST * 
XENEWCT«EKP«EKP1 « EKP2 . EKP3 , EKP4 . EKP5 « ENEWTH , THCl »ENEWM, 

XAA »<AC«NRET,SW2.EN(9 ) ♦ PRNOPT « HOLD , J , K «L , STOR ( 1000 ) 

IF(XC-XA-,01 )2« 1 * 1 
2 NRET=6 
RETURN 
1 X1=XA+.01 

TEMP=XC-XA 
TEMP=. 01/TEMP 
YI = ( YC-YA )*TEMP+YA 
TH1= (THC-THA )*TEMP+THA 
W 1 = ( WC-WA )*TEMP+WA 
SI = (SC-SA )*TEMP+SA 
CTHl =COS (THI ) 

STHI=S1N(THI ) 

TEMPI =WI*WI 
TEMPI =DEC1 /TEMPI 
TEMP2=CNTRL3-DEC 1 

SMUI = ( (TEMPl-DECl ) /DEC2 ) *TEMP2 
CMUI =SQRT (DECl -SMUI ) 

SMUI =SQRT (SMUI ) 

NRET= 1 
RETURN 
END 

SIBFTC SHKPT REF 

SUBROUT I NESHKPT 

COMMONXA « YA , THA « WA . SMUA. CMUA « STH A » CT HA , S A « X3 » YB » THB » 

1 WB.SMUBiCMUB«STHB»CTHB,SB,EPS«TEPB,XC»YC,THC.WC» 
2SMUC,CMUC»STHC«CTHC«SC,EPC«TEPC»EMUA»EMUB,EMUC«THAC* 

3WAC , SMUAC . YAC » EMUAC ♦ STHAC . EL AC « THBC ♦ WBC « SMUBC ♦ Y3C » EMUBC , 

4STHBC » EPBC « TEPBC,, EMBC « TANM A « TANMB , EMC ♦ PC « CPC « EM « THCP » 

5C0S1 »TAN1 «C0S2»TAN2»K0UNT«TAPE» ALIGHT«BEGS«SVC0N»SVC0N1 « 

6SVC0N2 * SVCON3 « SVC0N4 ♦ SVC0N5 , SVC0N6 , SVC0N7 « SVC0N8 « SVC0N9 « SVCO 1 0 ♦ 
7SVC01 1 . EEP«EEP1 « EEP2 * EEP3 * EEP4 , EEP5 « EEP6 « EEP7 « CNTRL « CNTRL 1 » 
8CNTRL2» CNTRL3«CNTRL4 «CNTRL5«CNTRL6*KNTRL7,CNTRL8«CNTRL9* CNTRl Ot 
9KNTR 1 1 . CNTR I 2 . CNTR 1 3 t CNTR 1 4 , CNTR 1 5 ♦ CNTR 1 6 » CNTR 1 7 « CNTR 1 8 « CNTR 1 9» 
XCNTR54, ERASE* ERAS 1 . ERAS2 * ERAS3 ♦ ERAS4 .ERAS5 . ERAS6 * ERAS7 * 
XERAS8*ERAS9«ERAS10«ERAS1 1 « ER AS 1 2 « ERASl 3 * ERAS 1 4 , ERAS 1 5 , ERAS 1 6 
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COMMONERAS 1 7 « ERAS 1 8 . ERAS 1 9 « ERAS20 « ERAS2 1 » ERAS22 » ERAS23« 
1ERAS24 *ERAS25«ERAS26«ERAS27»ERAS28.ERAS29»ERAS30.ERAS31 ♦ 
2ERAS32 , ERAS33 , ERAS34 ,ERAS35 « DEC « DEC 1 ♦ DEC2 » DEC3 » DEC4 . 
3DEC5*DEC6»DEC7»DEC8»DEC9.DEC10»DEC1 1 »KON.CONl .PRT* 

4PRT1 ,PRT2«PRT3»PRT4«PRT5»PRT6«PRT7iPRT8»PRT9«DEG» 
5TEMP»TEMP1 ♦ TEMP2 ♦ TEMP3 » TEMP4 » X I ♦ YI .THI »WI »SMUI » 

6CMUI ,STHI ,CTH1 ,SI ,KER10*EX1 ♦ EX2 » C 1 , C2 * C3 ♦ C4 . C5 , C6 ♦ 
7C7.C8»C9,C1 0«C1 1 « C 1 2 , C 1 3 » C 1 4 i C 1 5 « C 1 6 « C 1 7 , C 1 8 . C 1 9 , C20 « 

8C21 .C22,C23,C24,C25»C26,C27,C2a,C29«C30«SAv*SAVl « 

9SAV2 , S AV3 , SA V4 , S AV5 » SA V6 » S A V7 . SAV8 « SAV9 ♦SCI. THCPR ♦ ENEWST ♦ 
XENEWCT.EKP.EKPl .EKP2 ♦ EKP3 . EKP4 ♦ EKP5 . ENEWTH , THC 1 .ENEWM . 

XAA .K AC ♦ NRET ♦ SW2 ♦ EN ( 9 ) ♦ PRNOPT ♦ HOLD ♦ J ♦ K . L ♦ STOP (1000) 

COMMON LTEST.6SHK1 

CONTINUE 

CONTINUE 

THAC=THA 

WAC=WA 

ERASF=SMUA/CMUA 
EMUA=ATAN (ERASE ) 

EMUAC=EMUA 

EPBC=EPB 

TEPBC=TEPB 

K0UNT=0 

ERASE=STHA/CTHA 
CALLEPSR 
CALLTHASR 
6 ERAS6 = S IN ( THAC ) 

ERAS7 = C0S ( THAC ) 

ERAS8=S I N (EMUAC ) 

ERAS9=C0S (EMUAC ) 

ERAS 1 0 = COS ( EPBC ) 

eras 1 0= (S I N (EPBC ) ) /ERASl 0 

ERASl 1 =ERAS8XERAS9 

ERASl 3=ERAS6*ERAS8 

IF (BSHKl )40.40 .41 

40 ERASl 3=ERAS7*ERAS9-ERAS1 3 
GO TO 42 

41 ERASl 3=ERAS7*ERAS9+ERAS1 3 

42 ERAS12=ERAS6*ERAS9 
IF (BSHKl )50.50.5l 

50 ERAS12=ERAS7*ERAS8+ERAS12 
GO TO 52 

51 ERAS12=-ERAS7*ERAS8+ERAS12 

52 ERAS14=ERAS1 2/ERASl 3 
XC=ERAS14-ERAS1 0 
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YC=-XB*ERAS1 0+Y8 

XC= ( XA^ERASl 4-YA+YC ) /XC 

YC=XC*ERAS10+YC 

YAC= ( YC+YA )/DEC2 

IF (CNTRL8 ) 1 » 2» 1 

2 ELAC=0 
G0T03 1 

1 ELAC=( ( (ERAS8*ERAS8 ) /ERAS9 )*ERAS6 )/ERASl 3 
ELAC=( ( XC-XA )/YAC )*ELAC 
31 ERAS15=ERAS4/WAC 
IF (BSHKl )60. 60*61 

60 ERAS15=ERAS15-ERAS1 1 
GO TO 62 

61 ERAS15=ERAS15+ERAS1 1 

62 ERAS16= (ERA53-WA )/WAC 
ERAS16=ELAC-ERAS1 6 

1 F ( CNTRL9 ) 4 * 3 , 4 

4 ERAS18=CNTRL3*CNTRL5 

ERAS I 8= (ERAS8*ERAS8 ) /ERASl 8 
ERAS15=ERAS1 8*E RA S5+ ER A S 1 5 
ERAS16= (SA-SC )*ERAS1 8+ERAS16 

3 ERAS15=ERAS16/ERAS1 5 
THCP=THC 
THC=THA+ERAS 1 5 

CONT I NUE 

WC=ERAS4*ERAS15+ERAS3 
CALLMUSR 
IF (DEC )21 *20.21 
20 EMUAC= ( EMUC+EMUA ) /DEC2 
WAC= (WA+WC )/DEC2 
THAC=(THA+THC)/DEC2 
ERAS7=C0S ( THC ) 

ERAS6=S I N ( THC ) 

ERASE =ERAS 6/ERA S7 
CALLFPSR 

EPBC= (EPC+EPB )/DEC2 

KOUNT=KOUNT+KON 

IF (K0UNT-KNTRL7 )6* 5* 5 

5 IF ( ABS ( THC-THCP )-CNTRL6 ) 7* 7* 6 

7 IF(CNTRL9)9*8*9 
9 CALLSSR 

GOTOlO 

8 SC=SA 
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GOTOl 1 

10 CONTINUE 

11 STHC=SIN(THC) 

CTHC=COS (THC ) 

CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
21 RETURN 
END 

SIBFTC MOVAC REF 

SUBROUT INEMOVAC 

COMMONXA* YA.THA » WA,SMUA,CMUA «STHA.CTKA,SA.XB« YB»THB* 

1 WB » SMUB * CMUB ♦ STHB ♦ CTHB , SB , EPB . TEPB « XC » YC » THC « WC « 

2SMUC « CMUC « STHC « CTHC » SC « EPC « TEPC . EMUA . EMUB » EMUC » THAC *' 

3 WAC , SMUAC , YAC , EMuAC « STHAC , EL AC . THBC « WBC » SMUBC » YdC . EMUBC » 

4STHBC» EPBC. TEPBC ,e.M8C» TANMA • TANMB »EMC«PC « CPC * EM » THCPt 
5C0S1 «TAN1 *C0S2.TAN2»K0UNT,TAPE«ALIGHT»BEGS»SVC0N*SVC0N1 , 
6SVCON2«SVCON3.SVCON4 »SVC0N5»SVC0N6»SVC0N7» SVCON8»SVCON9,SVC010» 
7SVC01 1 1 EEP.EEPl . EEP2 » EEP3 i EEP4 « EEP5 * EEP6 ♦ EEP7 , CNTRL « CNTRL 1 * 
6 CNTRl 2«CNTRL3«CNTRL4 «CNTRLB«CNTRL6«KNTRL7*CNTRLt).CNTRL9.CNTR10« 
9KNTR1 1 ,CNTR12.CNTR13*CNTRl4,CNTKl5*CNTR16,CNTR17f CNTR18 ,CNTk 19. 
XCNTR54. ERASE » ERAS 1 ♦ ERAS2 * ERA S3 » ERAS4 « ERAS5 » ER AS6 * ERAS7 « 
XERASS,ERAS9,ERAS10tERASl 1 , ER ASl 2 , ERAS 1 3 , ERAS 1 4 ,ER AS 1 5 ♦ ERAS 1 6 
COMMONERAS 1 7 » ERAS 1 8 * ERAS 1 9 « ERAS20 * ERAS2 1 . ERAS22 « ERAS23 « 

1 ERAS24 »ERAS25«ERAS26 *ERAS2 7«ERAS28*ERAS29»ERAS30«ERAS31 ♦ 

2ERAS32 » ERAS33 , ERA534 » ERAS35 . DEC » DEC 1 « DEC2 « DECS * DEC4 , 
3DEC5,DEC6*DEC7»DEC8.DEC9.DEC10*DEC1 1 «K0N«C0N1 «PRT» 

4PRT1 ,PRT2*PRT3fPRT4«PRT5«PRT6.PRT7«PRT8»PRT9»DEG« 

BTEMP. TEMPI » TEMP2, TEMPS, TEMP4 ,XI , Y I ,THI , W I , SMU I , 

6CMU1 fSTHI ,CTH1 ,SI ,I<ER10,EX1 , EX2 , C 1 , C2 , C3 , C4 , C5 , C6 , 
7C7,C8,C9,C10,C1 1 ,Cl2,C13,C14,C15,C16,C17,Cl8,C19,C20, 

8C21 ,C22,C23,C24,C25,C26,C27,C28,C29,C30,SAV,SAV1 , 
9SAV2,SAV3,SAV4 ,SAV5,SAV6,SAV7,SAV8,SAV9,SC1 , THCPR , ENEWST , 
XENEWCT,EKP,EKP1 , EKP2 , EKP3 , EKP4 , EKP5 , ENEWTH , THC 1 ,ENEWM, 

XAA ,KAC , NRET , SW2 , EN ( 9 ) , PRNOPT , HOLD , J , K ,L , STOP (1000) 

XC = XA 
YC = YA 
THC=THA 
WC = WA 
SMUC=SMUA 
CMUC=CMUA 
STHC=STHA 
CTHC=CTHA 
SC = SA 
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RETURN 

END 

SIBFTC M0V3C REF 

SUBROUT I NEMOVBC 

COMMONXA * YA $THA » WA, SMUAt CMUA » STHA , CTHA . S A , X3 t YS« TH8 « 

1 W8.SMUB«CMUB*STHB.CTHB,SB,EP3.TEPB*XC»YC«THC«WC» 

2SMUC » CMUC * STHC » CTHC « SC » EPC » TEPC » EMUA « EMUB « EMUC « THAC « 

3WAC , SMUAC , YAC , EMUAC « STMAC , EL AC, THBC , WBC , SMuBC , YoC , EMUBC , 
4STHBC,EPBC,TEPBC,bMBC,TANMA, TANM8,EMC,PC,CPC,EM,THCP, 

5C0S1 ,TAN1 ,C0S2,TAN2,K0UNT, TAPE, ALIGHT, BEGS, SVCONfSVCONl , 
6SVC0N2 , SVCON3,SVCON4 ,SVC0N5,SVC0N6,SVC0N7, SVC0N6,SVC0N9,SVC01 0 
7SVC01 1 , EEP,EEP1 , EEP2 , EEP3 , EEP4 , EEP5 , EEP6 , EEP7 , CNTRL , CNTRL 1 , 
8CNTRL2, CNTRL3,CNTRL4 ,CNTRL5, CNTRL6,<NTRL7,CNTRL6,CNTRL9, CNTRl 0 
9KNTR1 1 ,CNTR12,CNTR13,CNTR14,CNTR15,CNTRI6,CNTR1 7 , CNTR 1 8 , CNTR 1 9 
XCNTR54 , ERASE,ERAS1 , ERAS2 , ERA S3 , ERAS4 , ERAS5 , ERAS6 , ERAS7 , 
XERAS8,ERAS9,ERAS1 0,ERAS1 1 , ER AS 1 2 , ER AS 1 3 , ERAS 1 4 , ER AS 1 8 , ERAS 1 6 
C0.MM0NERAS17,ERAS1B,ERAS19,ERAS20,ERAS2 1 , ERA S22 , ER AS23 , 

1 ERAS24 ,ERAS25 .ERAS26 ,ERAS27,ERAS28,ERAS29,ERAS30 ,ERAS31 , 
2ERAS32 , ERAS33 ,ERAS34 ,ERAS35 , DEC , DEC 1 , OEC2 , DECS , DEC4 , 
3DEC'D,DEC6,OEC7,DEC8,nEC9,l.)EC10,D£Cl 1 ,KON,CONl ,PRT, 

4PRT1 ,PRT2,PRT3,PRT4,PRTB,PRT6,PRT7,PRT8,PRT9,DEG, 

5TEMP , TEMPI , TEMP2 , TEMPS , l'EMP4 ,XI , Y I , TH I , W I , SMU 1 , 

6CMUI ,STHI ,CTHI ,S1 ,KER10,EX1 , EX2 , Cl , C2 , C3 , C4 , C5 , C6 , 
7C7,CS,C9,C10,C11,C12,C13,C14 , C 1 6 , C 1 6 , C 1 7 , C 1 8 , C 1 9 , C20 , 

8C2 1 , C2 2 ,C23 , C24 , C25 , C£6 , C2 7 , C28 , C29 , C30 ,SAV , SAVl , 
9SAV2,SAV3,SAV4 ,SAV5,SAV6,SAV7,SAV8,SAV9,SCl , THCPR , ENEWST , 
XENEWCT,£KP,EKP1 , EKP2 , EKP3 , EKP4 , EKP5 , ENE WTH , THC 1 ,ENEWM, 

XAA ,KAC ,NRET , SW2 ,EN ( 9 ) ,PRNOPT ,HOLD , J ,K ,L , STOP ( 1000) 

XC = XB 

YC = YB 

THC=TH3 

WC = WB 

SMUC=SMUB 

CMUC=CMUB 

STHC=STHB 

CTHC=CTHB 

SC = SB 

RETURN 

END 

SIBFTC 0ODYX REF 

SUBROUT I NEBODYX 

COMMONXA ,YA , THA , WA , SMUA , CMUA , STHA , CTHA , SA , X6 , YB , THB , 

1 WB ,SMUB ,CMUB ,STHB ,CTHB, SB , EPB , TEPB , XC , YC , THC , WC, 
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2SMUC , CMUC » STHC , CTHC » SC , EPC ♦ TEPC » EMUA . EMUB « EMUC « THAC t 
3WAC » SMUAC » YAC »EMUAC « STHAC , ELAC* THBC » WBC « SMUSC * YdC» EMUBC » 

4STH8C » EPBC * TEPBC ,£MBC » TANMA » TANM8 «EMC »PC ♦ CPC * EM « THCP» 

5C0S1 »TAN1 ,C0£2»TAN2«K0UNT,TAPE«ALIGHTfBEGS,SVC0N,SVC0Nl , 

6SVC0N2 , SVC0N3, SVC0N4 ,SVC0N5 < SVC0N6 » S VCON7 , SVC0N8 , SVCON9 1 SVCO 1 0 » 
7SVC01 1 , EEP*EEP1 . EEP2 t EEP3 % EEP4 t EEP5 » EEP6 t EEP7 , CNTPL t CNTPL 1 * 
8CNTRL2» CNTPL3, CNTPL4 , CNTRL5 , CNTRL6 * KNTRL7 * CNTRL8 . CNTRL9 , CNTR 1 0 » 
9KNTR1 1 ,CNTR12*CNTR13,CNTR14.CNTR1B,CNTR16.CNTR17.CNTR18,CNTR19. 
XCNTR54 , ERASE » ERAS 1 , ERAS2 . ERA S3 . ERAS4 . ERAS5 , ER AS6 » ERAS7» 
XERAS8,ERAS9,ERAS1 0,ERAS1 1 , ERASl 2 , ERAS 1 3 , ERAS 1 4 ,ER AS 1 5 , ER AS 1 6 
C0MM0NERAS17»ERAS18,ERAS19»ERAS20»ERAS21 . ERAS22 » ERAS23 » 

1 ERAS24 , ERAS25 ,ERAS26 i ERAS27 . ERAS28 , ERAS29 « ERAS30 » ERAS3 1 , 
2ERAS32«ERAS33 .ERAS34 . ERA S3B » DEC ♦ DEC 1 , DEC2 » DEC3 t DEC4 » 
3DEC5tDEC6.0EC7,DEC8*DEC9,DEC10»DECl 1 »KON,CONl tPRT. 

4PRT1 ,PRT2 tPRT3 »PRT4 , PRT5 ♦ PRT 6 » PRT7 . PRT8 . PRT9 * DEG . 

5 TEMP , TEMP 1 , TEMP2 , TEMP3 , TEMP4 » X I , Y I . TH I , 'ji/ 1 » SMU I » 

6CMUI ,SThI ,CTHI ,SI ,KER10»EX1 » EX2 ,C1 »C2»C3,C4»C5,C6« 
7C7.C8»C9,C1C»C1 1 »C12*C13»C14»C15*C16%C17.Cl8 , Cl 9, C20 » 

8C2 1 ,C22 , C23, C24 ,C25, C2 6, C2 7, C28, C29,C30 .SAVtSAVl . 
9SAV2»SAV3,SAV4,SAV5»SAV6» SA V7 * SA V8 , S AV9 . SCI « THCPR.ENEWST t 
XENEWCT,EKP,EKP1 , EKP2 » E KP3 , EKP4 ♦ EKP5 » ENEWTH , THC 1 ,ENEWM, 

XAA .KAC . NRET , SW2 . EN ( 9 ) , PRNQPT ♦ HOLD ♦ J , K . L . STOP ( 1000 ) 

THC=ATAN (CNTRl 6 ) 

ERASE =SMUB/CMUB 

ERASl =STHB/CTHB 

ERAS2=-ERAS1 *ERASE+DEC 1 

ERASE= (ERASl+ERASE )/ERAS2 

ERASl =ERASE-CNTR1 6 

XC= (ERASE*XB-YB+CNTR15)/ERAS1 

YC=XC*CNTR16+CNTR1 5 

RETURN 

END 

SIBFTC RDB2 REF 

SUBROUT I NERDB2 

COMMONXA, YA,THA,WA,SMUA,CMUA ♦ STH A , CT HA , S A , xB , YB ♦ THB « 

1 WB»SMUB»CMUB»STHB.CTHB,SB,EPB«TEPB»XC»YC *THC*WC. 

2SMUC 1 CMUC 1 ST HC » CTHC 1 SC » EPC « TEPC t EMUA . EMUB ♦ EMUC ♦ THAC * 

3WAC,SMUAC, YAC«EMUAC,STHAC,ELAC*THBC» WBC»SMU3C«YBC»EMUBC. 
4STHBC»EPBCiT£PBC«EMBC» TANMA , TANMB « EMC » PC » CPC ♦ EM , THCP* 

5C0S 1 t TANl , C0S2 » TAN2 1 KOUNT , TAPE » AL I GHT »BEGS , SVCON . SVCONl « 

6SVC0N2 » S VC0N3 . SVC0N4 » SVC0N5 , SVC0N6 « SVC0N7 . SVC0N8 » SVC0N9 . SVCO 1 0 « 
7SVC01 1 » EEP*EEP1 , EEP2 , EEP3 i EEP4 « EEP5 « EEP6 » EEP7 , CNTRL i CNTRL 1 « 
8CNTRL2 , CNTRL3 .CNTRL4 , CNTRLB « CNTRL6 « KNTRL7 « CNTRL8 , CNTRL9 , CNTR 1 0 » 
9:<NTR1 1 ,CNTR12»CNTR13,CNTR14.CNTR15*CNTR16«CNTR17.CNTR18»CNTR19« 
XCNTR54 , ERASE i ERASl , ERAS2 » ERA S3 « ERAS4 * ERAS5 , ERAS6 « ERAS7 ♦ 
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XERAS8,ERAS9,ERAS10*ERAS1 1 , ER AS 1 2 * ER AS 1 3 . ERAS 1 4 , ER A S 1 5 , ERAS 1 6 
COMMONERASl 7«ERAS1 8,ERAS1 9*ERAS20*ERAS21 , ERAS22 » ERAS23 . 

1 ERAS24 ,ERAS25«ERAS26 »ERAS27,ERAS28»ERAS29,ERAS30*ERAS31 , 
2ERAS32.ERAS33,ERAS34 , ERAS35 , DEC , DEC 1 » DEC2 , DEC3 . DEC4 « 
3DECS.DEC6.DEC7»DEC8*DEC9»DEC10*DEC1 1 .KON.CONl «PRT» 

4PRT1 ,PRT2,PRT3,PRT4 .PR T5*PRT6* PRT7.PRT8.PRT9, DEG. 

5TEMP, TEMPI ,TEMP2, TEMP3, TEMP4 .XI . YI ,THI .W1 . SMU 1 . 

6CMUI . STH 1 , CTH I ,SI .KERl 0 . EX 1 . EX2 .Cl , C2 , C3 , C4 . C5 , C6 , 
7C7.C8.C9.C10.C11,C12,C13.C14.C15.C16,C17.C18,C19,C20. 

8C21 ,C22.C23.C24.C25.C26.C27,C28,C29.C30.SAV.SAV1 . 

9SAV2 . SAV3. SAV4 . SAV5 . SAV6. SAV7.SAV8. SAV9. SCI . THCPR .ENEWST . 
XENEWCT .EKP.EKPl ,EKP2 . EKP3 , EKP4 . EKP5 . ENEWTH . THC 1 .ENEWM, 

XAA .KAC , NRET , SW2 .EN ( 9 > .PRNOPT .HOLD . J . K . L . STOR (1000) 

IF (STOR ( J+3 ) )2 . 1 . 2 
2 XB=STOR(J) 

YB=STOR ( J+1 ) 

THB=STOR ( J+2 ) 

WB=STOR (J+3 ) 

SMUB = STOR ( J + 4 ) 

CMUB = STOR ( J+5 > 

STHB=STOR ( J+6 ) 

CTHB = STOR ( J + 7 ) 

5B=ST0R ( J+8 ) 

J = J+9 

NRET=2 

RETURN 

1 KAC=KNTR1 1 +KON-KER1 0 
KERl 0=0 
NRET= 1 
RETURN 
END 

SIBFTC MOVCA REF 

SUBROUT I NEMOVCA 

COMMONXA , YA , THA , WA, SMUA , CMUA , STHA , CTHA , SA , X8 , YB . THB . 

1 WB . SMUB ,CMUB . STHB , CTH3 , SB , EPB . TEPB . XC . YC . THC . WC . 

2SMUC , CMUC .3THC . CTHC , SC . EPC . TEPC. EMUA , EMUd , EMUC , THAC . 

3WAC . SMU AC . YAC .EMU AC . STHAC , EL AC . T H8C . WBC . SMU3C . YbC . EMUBC . 

4STHBC , EPBC , TEPBC , EMBC. TANMA , TANMd .EMC .PC .CPC . EM , THCP. 

5C0S1 ,TAN1 .COS2.TAN2.KOUNT, tape. alight, begs, SVCON.SVCON l , 
6SVCON2. SVCON3 . SVCON4 .SVCON5 . SVCON6 . S VCON7 . SVC0N8 . SVC0N9 . S VCO 1 0 
7SVC01 1 . EEP.EEPl . EEP2 . EEP3 . EEP4 . EEP5 . EEP6 . EEP7 . CNTRL . CNTRL 1 . 
8CNTRL2 . CNTRL3. CNTRL4 .CNTRL 5. CNTRL6 . KNTRL7 . CNTRL8 . CNTRL9 , CNTR 1 0 
9KNTR1 1 . CNTR 1 2 . CNTR 1 3 , CNTR 1 4 . CNTR 1 5. CNTR 1 6 . CNTR 1 7 . CNTR 1 8 , CNTR 1 9 
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XCNTR54 t ERASE » ERAS 1 » ERASE « ERAS3 » ERAS4 » ERAS5 « ERAS6 » ERAS7 * 
XERAS8,ERAS9.£RAS1 OtERASl 1 « ERAS 1 2 , ERAS! 3 « ERAS 1 4 ,ERAS 1 5 . ERAS 1 6 
COMMONERAS 1 7 « ERAS 1 8 « ERAS 1 9 « ERAS20 » ERAS2 1 » ERAS22 « ER AS23 ♦ 

1 ERAS24 , ERAS25 « ERAS26 i ERAS27 , ERAS28 « ERAS29 , ERAS30 * ERAS31 » 
2ERAS32,ERAS33»ERAS34 ,ERAS35,DECf OECl ♦ DEC2 » OEC3 , DEC4 . 
3DEC5«DEC6»0EC7«DEC8»DEC9*DEC10»DEC1 1 «KON»CONl »PRT» 

4PRT1 ,PRT2»PRT3*PRT4,PRT5»PRT6«PRT7«PRT8«PRT9«DEG« 

5TEMP, TEMPI t TEMP2 , TEMPS , TEMP4 . X I » YI ,THI « W I » SMU I « 

6CMUI «STHI ,CTHI ,SI .KERlO»EXl » EX2 « C 1 , C2 « C3 « C4 . C5 , C6 » 
7C7.C8*C9»C10»C11.C12«C13«C14«C13»C1 6 . C 1 7 » C 1 8 , C 1 9 » C20 » 

3C2 1 . C22 . C23 « C24 * C25 « C26 , C2 7 , C23 « C29 . C30 « SAV * SA VI . 
9SAV2*SAV3«SAV4«SAV5»SAV6.SAV7.SAV8.SAV9.SC1 . THCPR « ENEWST « 
XENEWCT.EKP,EKP1 . EKP2 » EKP3 ♦ EKP4 » EKP5 « ENEWT H , T HC 1 *£NEWM, 

XAA ,KAC«NRET» SW2 .EN ( 9 ) * PRNOPT « HOLD , J , K , L » STOP (1000) 

XA = XC 

YA = YC 

THA=THC 

WA = WC 

SMUA=SMUC 

CMUA=CMUC 

STHA=STHC 

CTHA=CTHC 

SA = SC 

RETURN 

END 

SIBFTC PRINT REF 

SUBROUTINEPRINT 

COMMONXA, YA , THA ♦ WA « SMUA,CMUA * STHA ♦ CTHA « SA , XB » YB » THB » 

1 WB«SMUB»CMUB»STHB*CTHB, SB,EPB*TEPB«XC* YC.THCi WC« 

2SMUC *CMUC tSTHC » CTHC » SC , EPC * TEPC » EMUA «EMUB « EMUC » THAC « 

3WAC » SMUAC « YAC «EMUAC . ST HAC « EL AC ♦ THBC * WBC « SMuBC » YdC » EMUBC » 

4STHBC» EPBC«TEPBC ,EM3C« TANMA , TANMB , EMC « PC » CPC . EM , THCP « 

5C0S1 ,TAN1 ,C0S2»TAN2»K0UNT,TAPE<ALIGHT.BEGS»SVC0N.SVC0N1 « 
6SVC0N2»SVC0N3*SVC0N4 « S VC0N5 » SVC0N6 » S VC0N7 * S VC0N8 . S VC0N9 » S VCO 1 0 ♦ 
7SVC01 1 , EEP,EEP1 i EEP2 . EEP3 « EEP4 « EEP5 , EEP6 » EEP7 » CNTRL » CNTKL 1 » 
8CNTRL2 «CNTRl3»CNTRL4 , CNTRlS » CNTRL 6 i KNTRL7 » CNTRLB « CNTRL9 » CNTR 1 0 t 
9KNTR1 1 .CNTR12»CNTR13»CNTRl4»CNTR15»CNTR16«CNTR17»CNTR18,CNTR19i 
XCNTRB4 » ERASE 1 ERAS 1 , ERA S2 * ERA S3 » ERAS4 » ERAS5 » ERAS6 « ERAS7 « 
XERAS8»ERAS9*ERAS10iERASl 1 , ERASl 2 . ERAS 1 3 * ERAS 1 4 , ER AS 1 5 * ERAS 1 6 
COMMONERAS 1 7 »ERAS1 8 » ER AS 1 9 » ERAS20 « ERAS2 1 % ERAS22 , ERAS23 » 

1 ERAS24 ,ERAS25,ERAS26 «ERAS2 7,ERAS28,ERAS29,ERAS30«ERAS31 , 
2ERAS32,ERAS33«ERAS34 «ERAS35» DEC* OECl * DEC2 ♦ DECS , DEC4 « 
3DEC5,DEC6,DEC7*DEC8*DEC9*DEC10*DEC1 1 *<0N*C0N1 *PRT* 

4PRT1 ,PRT2 ,PRT3*PRT4 *PRT5*PRT6*PRT7« PRT8*PRT9*DEG« 

5TEMP. TEMPI . TEMP2 , TEMPS , TEMP4 * X I * YI *TH1 ,WI *SMUl , 
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6CMUI .STHI tCTHI ,SI »I<ER10»EX1 ♦ EX2 » C 1 « C2 » C3 » C4 » C5 , C6 « 
7C7»C8»C9«C10iCl 1 « C 1 2 « C 1 3 « C 1 4 » C 1 o » C 1 6 » C 1 7 » C 1 8 « C 1 9 , C20 ♦ 

8C21 «C22»C23»C24 i C25 « C26 * C27 , C28 » C29 » C30 . SAV » SA V 1 . 

9SAV2 , SAV3 . SAV4 fSAVSi SAV6 ^ SAV7<SAV8<SAV9<SC1 < THCPR < ENEWST < 
XENEWCT<E<P<EKP1 < EKP2 < E<P3 < EKP4 < EKP5 < ENEWTH < THC 1 <.ENEWM< 

XAA <KAC < NRET < SW2 < EN ( 9 ) < PRNOPT<HOLD<J< K <L < STOP (1000) 
ERASE=CNTRL3-DEC1 
ERAS2=ERASE/DEC2 
ERAS3=ERAS2*EEP5+DEC 1 
ERAS4=-CNTRL 3/ERASE 
ERAS5= (CNTRL3*EEP5)/DEC2 
EMC=DEC1/SMUC 

ERASE= ( EMC*EMC*ERAS2+DEC1 )/ERAS3 
ERASE = EXP ( ALOG ( ERASE )*ERAS4 ) 

ERASl =-SC/CNTRL5 
PC=EXP(ERAS1 )*ERASE 
CPC=(PC-DEC1 )/ERAS5 
THD=THC*DEG 

WRITE (6 < 100 ) XC<YC<THD<EMC<WC<SC<PC<CPC 
lOO FORMAT (8E16. 8 ) 

RETURN 

END 

SIBFTC MUSR REF 

SUBROUT I NEMUSR 

COMMONXA <YA<THA<WA<SMUA<CMUA<STHA<CTHA<SA<XB<YB<THB< 

1 WB < SMUB < CMUB <STHB<CTHB< S6< EPB< TEPB < XC < YC <THC < WC< 
2SMUC<CMUC<STHC<CTHC<SC<EPC<TEPC<EMUA<EMUB<EMUC<THAC< 

3WAC<SMUAC< YAC<EIMUAC<STHAC<ELAC<THBC< WBC<SMUBC <YBC<EMUBC< 

4STHBC <EPBC<TEPBC<EMBC<TANMA< TANMB<EMC<PC<CPC<EM < THCP< 

5C0S1 < TANl <C0S2<TAN2<K0UNT< TAPE<ALIGHT<BEGS<SVCON < SVCONl < 

6SVC0N2 < SVC0N3 < SVC0N4 < SVCON5 < SVCON6 < SVC0N7 < S VCON8 < SVCON9< SVCO 1 0< 
7SVC01 1 < EEP<EEP1 <EEP2<EEP3<EEP4<EEP5<EEP6<EEP7<CNTRL<CNTRL1 < 
8CNTRL2 < CNTRL3 <CNTRL4 <CNTRL5<CNTRL6<KNTRL7<CNTRL8<CNTRL9<CNTR10< 
9KNTR1 1 <CNTR12<CNTR13 <CNTR14<CNTR15<CNTR16<CNTR17<CNTR18<CNTR19< 
XCNTR54 < ERASE < ERASl <ERAS2<ERAS3<ERAS4<ERAS5<ERAS6<ERAS7< 
XERAS8<ERAS9<ERAS10<ERAS1 1 <ERAS12<ERAS13<ERAS14<ERAS15<ERAS16 
C0MMONERAS17<ERAS18<ERAS19<ERAS20<ERAS21 <ERAS22<ERAS23< 

1 ERAS24<ERAS25<ERAS26<ERAS27<ERAS28<ERAS29<ERAS30<ERAS31 < 

2ERAS32 < ERAS33 < ERAS34 <ERAS35 < DEC < DECl < DEC2 < DEC3 < DEC4 < 
3DEC5<DEC6<DEC7<DEC8<DEC9<DEC10<DEC1 1 <KON<CONl <PRT< 

4PRT1 <PRT2<PRT3<PRT4<PRT5<PRT6<PRT7<PRT8<PRT9<DEG< 

5TEMP< TEMPI <TEMP2<TEMP3<TEMP4 <X1 <YI <THI <Wl < SMU I < 

6CMU1 <STHI <CTHI ,SI <K£R1 0<EX1 <EX2<C1 <C2<C3<C4<C5<C6< 
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7C7 fCe* C9.C1 0 »C1 1 »C12»C13»C14 . Cl 5 » Cl 6 » C 1 7 . C 1 8 . C 1 9 « C20 * 

8C21 ♦C22»C23»C24»C25»C26*C27»C28»C29»C30»SAV«SAV1 « 

9SAV2 * SAV3 . SA V4 ♦ SAV5 » SA V6 « SAV7» SA V8 ♦ SAV9 « SC 1 » THCPR » ENEWST » 
XENEWCT«EKP*EKP1 , EKP2 » EKP3 ♦ EKP4 . EKP5 » ENEWTH , THC 1 .ENEWM, 

XAA » K AC t NRET , SW2 » EN ( 9 ) ♦ PRNOPT * HOLD , J . K » L ♦ STOR (1000) 

ERASE=WC*WC 

ERASE=DEC 1 /ERASE 

ERASl =CNTRL3-DEC1 

SMUC=( (ERASE-DECl )/DEC2 )*ERAS1 

IF(DEC1-SMUC)2, 1 . 1 

2 CMUC= (DECl-SMUC ) 

GO TO 3 

1 CMUC=SQRT (DECl -SMUC ) 

3 IF (SMUC )4 ,5* 5 

4 SMUC=SMUC 
GO TO 6 

5 SMUC=SQRT (SMUC ) 

6 ERASE=SMUC/CMUC 
EMUC=ATAN (ERASE ) 

RETURN 

END 

SIBFTC main REF 

C MA I NPROGRAM-SECONDSHOCKREG I ON 

COMMONXA,YA,THA»WA,SMUA,CMUA,STHA,CTHA,SA,xB, YB«THB, 

1 WB , SMUB » CMUB i STHB . CTH6 , SB , EPB « TEPB , XC » YC « THC « WC « 
2SMUC»CMUC.STHC,CTHC«SC»EPC*TEPCf£MUA,EMUB»EMUCfTHAC* 

3 WAC, SMU AC « YAC » EMUAC « STHAC » EL AC ♦ THBC<WBC<SMUSC<YOC<EMUBC< 
4STHBC<EPBC<TEPBC<EMBC<TANMA<TANMB<EMC<PC<CPC<EM<THCP< 

5C0S1 <TAN1 <C0S2<TAN2<K0UNT<TAPE<ALI GHT<BEGS<SVC0N<SVC0N1 < 

6SVC0N2 < SVC0N3 < SVC0N4 < S VC0N5 < SVC0N6 < S VC0N7 < SVC0N8 < SVC0N9 < SVCO 1 0 < 
7SVC01 1 < EEP<EEP1 < EEP2 < EEP3 < EEP4 < EEP5 < EEP6 < EEP7 < CNTRL < CNTRu 1 < 
8CNTRL2 < CNTRL3 < CNTRL4 <CNTRL5 < CNTRL6 < i<NTRL7 < CNTRLB < CNTRL9< CNTKl 0< 
9KNTR1 1 <CNTR12<CNTR13<CNTR14<CNTR15<CNTR16<CNTR17<CNTR18<CNTR19< 
XCNTR54 < ERASE < ERASl <ERAS2<ERAS3<ERAS4 <ERAS5<ERAS6<ERAS7< 
XERAS8<ERAS9<ERAS10<ERAS1 1 <ERAS12<ERAS13<ERAS14<ERAS15<ERAS16 
COMMONERAS 1 7 < ERAS 1 8 < ERAS 1 9 < ERAS20 < ERAS2 1 < ERA S22 < ERAS23 < 

1 ERAS24 < ERAS25 < ERAS26 < ERAS27 < ERAS28 < ERAS29 < ERAS30 < ERAS31 < 
2ERAS32<ERAS33<ERAS34<ERAS35<DEC<DEC1 <DEC2<DEC3<DEC4< 
3DEC5<DEC6<DEC7<DEC8<DEC9<DEC10<DEC1 1 <KON<CONl <PRT< 

4PRT1 <PRT2<PRT3<PRT4<PRT5<PRT6<PRT7<PRT8<PRT9<DEG< 

5TEMP< TEMPI <TEMP2<TEMP3<TEMP4<XI <YI<THI<WI <SMUl < 

6CMUI <STHI <CTHI <SI <KER10<EX1 < EX2 < C 1 <C2<C3<C4<C5<C6< 
7C7<C8<C9<C10<C1 1 <C12<C13<C14<C15<C16<C17<C18<C19<C20< 

8C21 <C22<C23<C24<C25<CE6<C27<C2B<C29<C30<SAV<SAV1 < 
9SAV2<SAV3<SAV4<SAV5<SAV6<SAV7<SAV8<SAV9<SC1 < THCPR < ENEwST < 
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XENEWCT «EKP ♦ EKP 1 ,EKP2 ^Ei<P3 , EKP4 ^ EKP 5 t ENE WTH , T HC 1 ,£NEWfvi , 

XAA ;kAC , NRET « SW2 i EN ( 9 ) ♦ PRNQPT i HOLD » J ^ K ♦ L « 5T0R (1000) 

TAPE=2# 

AL I GHT=0 
DEC = 0 
DECl = 1 • 

DEC2=2* 

DEC3= 1 • 570796327 

DEC4 =6. 2831 85308 

0EC5=0 

DEC6=0 

DEC7=7 

DEC8=0 

DEC9=9 

DECl 0=0 

DECl 1=11 

KOM= 1 

CONI =3, 

DEG=57. 2957795 
9 CALLSL I TE ( 0 ) 

WR I TE ( 6 ^ 98 ) 

9 8 FORMAT ( 1 HI , 1 HX 1 3X . 1 HY ♦ 1 3X ^ 9HTHETA DEG . 5X ♦ 1 HM , 1 3 X » 1 H W , 1 3 X ^ 1 H S » 1 3 X » 

1 6HPT/PT 1 1 X20H (PT/PT 1 ) M ( 1 - W ^*2 ) -«-* 3 1 X2 5H ( ( PT/PT 1 ) ) N ( 1 -W**2 ) *^<*3 ) 

DO 999 NN= 1^1000 
999 STOR(NN)=0 
HOLO= 1 • 

J= 1 
l<= 1 
L= 1 

CALLSL I TE ( 2 ) 

KERl 0=0 

CALL FINP(16^CNTRL^CNTRLlf CNTRL3 ^ CNTRL5 » CNTR 1 3 , CNTR 1 6 , EM 1 ^ Xb « YB , 
1XA»YA,EMB^THB^5LJ«EPB ,PRNOPT ) 

EPb=EPB/DEG 
THB=TH5/DEG 
CNTRl6=9999 • 

KNTRL7=3 
CNTRl8= 1 * 

CNTRL9= 1 « 

KNTR 11=0 
SW2 = C 

\a/B = SQRT (DECl / (DLCl+DEC2/{ ( CN TRL3-‘DEC 1 ) -5t-EMb-«-E Mb ) ) ) 

SMUB=DEC1 /EmB 
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SMUA=SMUB 

CMUB=SQRT ( 1 .-SMUB*SMU8 ) 
CMUA=CMUB 
STHB=SIN (THB ) 
CTHB=COS(THB) 

STHA=STHB 
CTHA=CTHB 
SA=SB 
WA = WB 
THA=THB 

TEPB=S IN (EPS I/COS (EPB ) 

BEGS=EPB 

CALLSH<PT 

IE (DEC )361 0 » 30 , 361 0 

30 CALL PRINT 
CALLWRTSK 
CALLMOVCB 
CALLBODY 

IF (DEC )361 0« 31 ,3610 

31 CALL Moves 
CALLMOVBA 
CALL I PT 

GOTO (22, 22, 22, 23) ,NRET 

22 CALLMOVIC 
CALLWRTGN 
CALLPRINT 

23 CALLMOVSC 
CALLWRTGN 
CALLPRINT 
CALLSPACE 
CALLWRTFL 

7 CALLRDB 

GOTO ( 10 » 1 1 ) ,NRET 

10 CALLEXIT 

1 1 CALLRDA 

GOTO (12,13), NRET 

12 CALLEXIT 

13 CALLSHKPT 

I F (DEC )361 0, 32 ,3610 

32 CALL PRINT 
CALLWRTSK 
CALLMOVCB 
CALLRDA 

GOTO( 14, 15) ,NRET 

14 CALLEXIT 
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15 CALLGEN 

IF (DEC >3610 *33 *3610 

33 GO TO 16 

17 CALLMOVBA 
CALLIPT 

GOTO (24 *24.24,20) *NRET 

24 CALLMOVIC 
CALLWRTGN 
CALLPR INT 

25 CALLMOVSC 
CALLWRTGN 
CALLPRINT 

4 CALLRDA 

GOTO( 18*19) *NRET 

18 GOTOl 

19 CALLMOVCB 
CALLGEN 

I F ( DEC ) 36 1 0 * 34 * 36 1 0 

34 CALL WRTGN 

IF (PRNOPT )3*2* 3 
3 CALLPRINT 
2 G0T04 
1 CALLMOVCB 
L = 1 

IF (EN(L+1 ) )20*21 *20 
21 XA=EN(L) 

YA=EN(L+1 ) 

THA=EN(L+2 ) 

WA=EN (L+3 ) 

SMUA=EN (L+4 ) 

CMUA=EN (L+5 ) 

STHA = EN (L + 6 ) 

CTHA=EN (L+7 ) 

SA = EN (L + 8 ) 

G0T020 

20 CALLBODY 

I F ( DEC ) 36 1 0 * 35 * 36 1 0 

35 CALL WRTGN 
CALLPRINT 

3610 DEC=Q 

CALL SPACE 
CALLWRTFL 

IF(XC-CNTRL1 )6*6*5 
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6 IF(SW2)8«7»8 
8 WRITE(6* 100 )TAPE, ALIGHT 
5 G0T09 
16 CALLMOVCS 
L=1 

EN (L ) =XA 
EN (L+1 ) =YA 
EN (L+2 ) =THA 
EN (L+3 ) =WA 
EN (L+4 ) =SMUA 
EN (L+5 ) =CMUA 
EN (L+6 ) =STHA 
EN (L+7 ) =CTHA 
EN (L+8 ) =SA 
GOTOl 7 

lOO F0RMAT(2E16.8) 

END 

SIBFTC NEWS REF 

SUBROUT I NENEWB 

COMMONXA.YA.THA*WA,SMUA,CMUA»STHA,CTHA,SA,XB»YB*THB« 

1 WB*SMUB»CMUB*STHB«CTHB,Sb,EP3«TEP6.XC*YC*THC*WC» 

2 SMUG . CMUC . STHC , CTHC . SC , EPC » TEPC » EMUA ♦ EMUri , EMUC , THAC * 

3WAC , SMUAC . YAC »EMUAC . STHAC , ELAC* THBC * WBC * SMUBC ♦ YdC * EMUBC » 
4STHBC*EPBC»TEPBC.EMBC«TANMA,TANMB,EMC*PC.CPC*EM»THCP. 

5C0S1 «TAN1 ,C0S2.TAN2«K0UNT*TAPE, alight » BEGS, SVCON.SVCON l » 

6SVC0N2 , SVCON3 , SVC0N4 , SVC0N5 , SVC0N6 , SVC0N7 , SVC0N8 , SVC0N9 « SVCO 1 0 < 
7SVC01 1 ♦ EEP,EEP1 ,'EEP2 , EEP3 , EEP4 , EEP5 , EEP6 , EEP7 ♦ CNTRL , CNTRL 1 , 
8CNTRL-2,CNTRL3,CNTRL4 , CNTRLb , CNTRL6 , i<NTRL7 , CNTRL8 , CNTRL9 ♦ CNTR 1 0 » 
9KNTR1 1 ,CNTR12,CNTR13,CNTR14,CNTR15,CNTR16,CNTR17,CNTR18,CNTR19, 
XCNTR54 , ERASE , ERAS 1 , ERA S2 , ERA S3 , ERAS4 , ERA S5 , ER AS6 , ERAS7 , 

XER AS8, ERAS9, ERAS 10, ERAS 1 1 , ERAS 1 2 , ERASl 3 , ERAS 1 4 , ER AS 1 5 , ERAS 1 6 
COMMONERAS17,ERAS18,ERAS19,ERAS20,ERAS21 , ERAS22 , ERAS23 , 

1ERAS24 ,ERAS25,ERAS26 ,ERAS27,ERAS28,ERAS29,ERAS30,ERAS31 , 
2£RAS32,ERAS33,ERAS34 , ER AS3b , DEC , DEC 1 , UEC2 , DEC3 , DEC4 , 
3DEC5,DEC6,DEC7,DEC8,DEC9,DEC10,DEC1 1 ,K0N,C0N1 ,PRT, 

4PRT1 ,PRT2,PRT3,PRT4,PRT5,PRT6«PRT7,PRT8,PRT9,DEG, 

5TEMP, TEMPI , TEMP2 , TEMP3 , TEMP4 , X I ,YI ,THI ,WI ,SMUl , 

6CMUI ,STHI ,CTH1 ,SI ,KER10,EX1 , EX2 , C 1 , C2 , C3 , C4 » C5 , C6 , 
7C7,C8,C9,C10iC1 1 , C 1 2 , C 1 3 , C 1 4 , C 1 5 , C 1 6 , C 1 7 , C 1 8 , C 1 9 , C20 , 

8C2 1 , C22 , C23 , C24 , C25 , C26 , C27 , C28 , C29 , C30 , SAV»SAVl » 
9SAV2,SAV3,SAV4,SAV5,SAV6,SAV7,SAV8,SAV9,SCl , THCPR , ENEWST , 
XENEWCT,EKP,EKP1 ,EKP2 , EKP3 , EKP4 , EKP5 , ENE WTH , THC 1 ,ENEWM, 
XAA,KAC,NRET,S';t/2,EN(9) , PRNOPT , HOLD , J ,<, L , STOP ( lOOO) 

SAV=CMUA/SMUA 
I F (STHA )5 ,4 , 5 
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4 SAV4=SMUA/CMUA 
G0T03 

5 SAVl =CTHA/STHA 
SAV3=DEC1 /SAVl 
IF (SAV )2* 1 

2 SAV2=DEC1 /SAV » 

SAV4=-SAV2*SAV3+nEC 1 
SAV4= ( SAV3 + SAV2 J/SAV4 
GOT03 

1 SAV4=-SAV1 

3 SAV5 = XB*TEPB-YB + YA. 

SAV5=-XA*SAV4+SAV5 

SAV6=TEPB-5AV4 

XC=SAV5/SAV6 

SAV5=XC-XB 

YC = SAV5^(-TEPB+Y3 
CALLSUBR 
ENEWTH=THC 1 +THA 
ENEWCT=C0S (ENEWTH ) 

ENEWST=S IN (ENEWTH ) 

EKP5 = ENEWM*e:nEWM 
EKPl =EI<P5 + EKP5 
EKP4=EKP1 +DEC1 
EI<P = F!<P5*FKP5 

EKP3=( (CNTRL3+DEC1 ) /DEC2 ) *El<P5 + DEC 1 

EKPl =-EKP3*EKP 1 

EKP2=- ( EKPB+DECa ) *EKPB 

RETURN 

END 

SIBFTC PRINT REF 

SUSROUT I NEPR I NT 

COMMONXA ,YA,THA»WA,SMUA,CMUAiSTHA»CTHA,SAiX3,YB,TMbi 
1 WB .SMUB , CMUB * STHB « CTHt3 , SB , EPB » TEPB , XC » YC » THC i WC i 
2SMUC t CMUC * STHC , CTHC , SC , EPC « TEPC » EMu A , E YUB » EMUC , T H AC i 
3WAC ,SMUAC * YAC ,EMUAC , STHAC , ELAC. THBC * WdC » SMuBC ^ Y-iC , EMUBC* 
4STHbC»EPBCiTEPBC,EMbC,TANMA,TANMB,FiV,c,PC«CPCiEM,ThCP* 

5C0S1 ,TAN1 »C0S2«TAN2iK0UNT,TAPE.AElGHT»bEGS, SVCON»SVCONl * 
6SVC0N2 * SVC0N3 , SVC0N4 , SVC0N5 , SVC0N6 «SVC0N7*SVC0N8, SVC0N9 i SVCO 1 0 
7SVC01 1 , EEP»EEP1 » EEP2 » EEP3 » EEP4 ♦ EEP5 » EEP6 . EEP7 , CNTRL » CNTRU 1 » 
8CNTRL2, CNTRL3 «CNTRL4 i CNTR lB , CNTRL 6 » KNTRL7 i CNTRL8 , CNTRL9 , CNTR 1 0 
9KNTR1 1 ,CNTR12.CNTR13.CNTRl4,CNTR15»CNTR16»CNTRl7,CNTRia.CNTRl9 
XCNTR54 » ERASE 1 ERASl * ERAS2 »ERAS3 » ERAS4 *£RAS5 ,ERAS6 * ERAS7 » 
XERAS8*ERAS9»ERAS10«ERAS1 1 »ERAS12.ERAS13.ERAS14,ERAS15iERAS16 
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C0MM0NERAS17*ERAS18«ERAS19*ERAS20tERAS21 ♦ ERAS22 » ERAS23 i 
lERAS24*ERAS25»ERAS26*ERAS27*ERAS28»ERAS29»ERAS30iERAS31 i 
2ERAS32«ERAS33»ERAS34 » ERAS35 i DEC * DEC 1 * DEC2 » DEC3 ♦ DEC4 i 
3DEC5 *DEC6.DEC7« DECS «DEC9,DEC10* DEC 1 1 ,K0N*C0N1 »PRT, 

4PRT1 iPRT2 «PRT3*PRT4 tPRT5»PRT6«PRT7* PRT8 ,PRT9 . DEGi 
5TEMP*TEMP1 . TEMP2 « TEMP3 » TEMP4 ♦ X I » YI iTHi ♦Wl »SMUI * 

6CMUI »STHI ♦ CTHI .SI .KERIO.EXI . EX2 . C 1 .C2.C3.C4.C5.C6. 

7C7.C8.C9. CIO. Cl 1.C12.C13.C14.C15.C16.C17.C18.C19.C20. 

8C21 .C22.C23.C24.C23.C26.C27. C28.C29.C30 .SAV.SA VI . 
9SAV2.SAV3.SAV4.SAV5.SAV6.SAV7.SAva.SAV9. SCI . THCPR . ENEWST . 
XENEWCT.EKP.EKPl .EKP2 . EKP3 . EKP4 . EKP5 . ENEWTH . THC 1 .ENEWM . 

XAA .KAC . NRET. SW2 .EN ( 9 ) . PRNOPT .HOLD.J.K.L.STORC 1000) 
EMC=DEC1/SMUC 
TEMP=-SC/CNTRL5 
TEMP=EXP(TEMP) 

TEMPI =-WC*WC+DECl 

TEMPI =TEMP1*TEMP1 *TEMP1*EMC* TEMP 

TEMP2=TEMP1 *TEMP 

THD=THC*DEG 

WRITE (6. 100 ) XC.YC.THD.EMC.WC.se. TEMP. TEMPI .TEMP2 
RETURN 

100 format (9E14. 6 ) 

END 

SIBFTC EPSR REF 

SUBROUT I NEEPSR 

COMMONXA. YA.THA.WA.SMUA.CMUA . STHA . CTHA . S A . XB . YB . THB . 

1 WB.SMUB.CMUB.STHB.CTHB.SB.EPB.TEPB.XC.YC.THC.WC. 
2SMUC.CMUC.STHC.CTHC.SC.EPC.TEPC.EMUA .EMUB.'EMUC.THAC. 

3WAC .SMUAC . Y AC .EMU AC . STHAC , EL AC . THBC . WBC . SMUBC . YdC. EMUBC. 
4STHBC.EPBC.TEPBC.EMBC.TANMA.TANMB.EMC.PC.CPC.EM.THCP. 

5C0S1 .TANl .C0S2.TAN2.K0UNT. TAPE. ALIGHT. BEGS. SVCON.SVCONl . 

6SVC0N2 . SVC0N3 . SVC0N4 . SVC0N5 . SVC0N6 . S VC0N7 . S VC0N8 . SVC0N9 . SVCOl 0. 
7SVC01 1 . EEP.EEPl . EEP2 . EEP3 . EEP4 . EEP5 . EEP6 . EEP7 . CNTRL . CNTRL 1 . 
8CNTRL2.CNTRL3.CNTRL4 .CNTRl5.CNTRL6.KNTRL7.CNTRL8.CNTRL9,CNTR10. 
9KNTR1 1 .CNTR12.CNTR13.CNTR14.CNTR15.CNTR16.CNTR17.CNTR18.CNTR19. 
XCNTR54 . ERASE. ERAS 1 . E RAS2 . ERA S3 . ERAS4 . ERAS5 . ERAS6 . ERAST . 
XERASS.ERAS9.ERAS10.ERAS1 1 . ER AS 1 2 . ERAS 1 3 . ERAS 1 4 . ER AS 1 5 . ERAS 1 6 
C0MM0NERAS17.ERAS18.ERAS19.ERAS20.ERAS21 . ERAS22 . ER AS23 . 

1 ERAS24.ERAS25.ERAS26 .ERAS27.ERAS28.ERAS29.ERAS30.ERAS31 . 

2ERAS32 . ERAS33 .ERAS34 .ERAS35 . DEC . DEC 1 . DEC2 . DECS . DEC4 . 
3DEC5.DEC6.DEC7.DEC8.DEC9.DEC10.DEC1 1 .KON.CONl .PRT. 

4PRT1 .PRT2.PRT3.PRT4.PRT5.PRT6.PRT7.PRT8.PRT9.DEG. 

5TEMP. TEMPI . TEMP2 . TEMP3 . TEMP4 . X I .YI .THI .Wl .SMUI . 

6CMUI .STHl .CTHI .SI .KERIO.EXI . EX2 . C 1 .C2.C3.C4.C5.C6. 

7C7.C8.C9. CIO. Cl 1.C12.C13.C14.C15.C16.C17.C18.C19.C20. 
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8C2 1 , C22 , C23 » C24 * C25 » C 26 . C27 , C 2 B , C29 » C30 * SAv » SA Vl » 
9SAV2»SAV3»SAV4.SAV5.SAV6»SAV7iSAV8*SAV9»SCl i THCPR , ENEWST i 
XENEWCT»EKP.EKP1 .EKP2 , E KP3 . EKP4 . EKP5 * ENEWT H , T HC 1 » ENEWM i 
XAA »KAC » NRET . SW2 »EN (9 ) , PRNOPT » HOLD , J » K « L » S TOR < 100 0 ) 
ERASE=ERASE*ERASE 
ERASl =ERASE*EKP+EKP 
ERAS2= ( ERASE*EKP1 +E<P2 >/CONl 
ERAS13= (EKP3*EKP3*ERASE+E<P4 )/CONl 
ERAS6=ERAS 1 3*ERAS2+E RA S 1 
ERAS9=ERAS2*ERAS2 
ERASl 5=ERAS1 *ERAS1 3-ERAS9 
ERASE=ERAS1 5+ERASl 5 
ERAS 1 0= ( ERAS 6/ERASE ) *ERAS1 +ERAS2 
ERAS6=ERAS6*ERAS6 
ERAS7=ERAS13*ERAS1 3+ERAS2 

AA= (SORT ( ABS ( (ERAS7 + ERAS7 ) *ERASE + ERAS6 ) > /ERASE )*ERAS1 
AA=ATAN2 (AA.ERASl 0 ) 

1 F ( AA )2 *2 , 1 
2 AA=AA+DEC4 
1 ERASE=AA/C0N1 
ERASE=COS (ERASE ) 

ERAS9 = SQRT ( -ERAS 1 5 ) 

ERASE = - ( (ERAS9 + ERAS9 )*ERASE + ERAS2 )/ERASl 
ERAS2 = SQRT ( DEC 1 -ERASE ) 

ERAS 1 =SQRT ( ERASE ) 

TEPC=ERAS1 /ERAS2 
EPC = ATAN (TEPC ) 

RETURN 

END 

*IBFTC BODY REF 

SUBROUT I NEBODY 

COMMONXA» YA,THA»WA,SMUA,CMUA * STHA , CTHA » S A , XB i YB » THB » 

1 WB *SMUB ,CMUBiSTHB ,CTHc , SB , EPB i T EPB * XC i YC . THC » WC » 

2SMUC * CMUC *STHC , CTHC » SC ,EPC » TEPC i EMUA » EMUB , EMUC , THAC * 

3WAC , SMUAC , YAC *EMUAC i STHAC , ELAC* THSC » WbC . SMUBC , YBC» EMUBC » 

4STHBC. EPBC* TEPBC , EMBC , TANMA , TANMB ,EMC tPC , CPC » EM , THCPi 
5C0S1 ,TAN1 *C0S2»TAN2»<0UNT,TAPE»ALIGHT»BEGS,SVC0N,SVC0N1 » 
6SVC0N2»SVC0N3«SVC0N4.SVC0N5»SVC0N6»SVC0N7,SVC0N8,SVC0N9,SVC01 0» 
7SVC01 1 * EEP»EEP1 » EEP2 » EEP3 « EEP4 t EEP5 i EEP6 « EEP7 , CNTRL , CNTRL 1 t 
8CNTRL2, CNTRL3 ♦CNTRL4 , CNTWL5 * CNTRL6 » KNTRL7 , CNTRL8 , CNTRL9 , CNTR 1 0 « 
9KNTR1 1 «CNTR12»CNTR13»CNTR14,CNTR15»CNTR16«CNTR17»CNTR18«CNTR19» 
XCNTR54, ERASE* ERASl * ERA S2 , ERA S3 « ER AS4 * ER A S5 * ER AS6 * ERAS 7 , 
XERAS8»ERAS9*ERAS10,ERAS1 1 , ER A S 1 2 , ERAS 1 3 * ERAS 1 4 , ER AS 1 5 , ER AS 1 6 
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C0MM0NERAS17»ERAS18 .ERAS19tERAS20»ERAS21 . ERAS22 « ERAS2 3 ♦ 

I ERAS24 , ERAS25 , ERAS26 » ERAS27 , ERAS28 t ERAS29 , ERAS30 « ERAS31 * 
2ERAS32* ERAS33 »ERAS34 *ERAS35 » DEC « DEC 1 »D£C2 » DEC3 ,DEC4 i 
3DEC5iDEC6.DEC7«DEC8 .DEC9.DEC10»DEC1 1 «i<ON»CONl »PRT* 

4PRT1 ,PRT2*PRT3 *PRT4 ,PRT5 *PRT6fPRT7 ,» PRT8 i PRT9 , DEG , 

5TEMP * TEMP 1 i TEMP2 , TEMP3 , TEMP4 i X I , Y I , TH I ♦ W I » SMU I i 
-6CMUI ,STHI,CTHI,SI .KER10«tXl » EX2 » C 1 , C2 , C3 , C4 i C5 , C6 * 

7C7 *C8 *C9*C10»Cll.C12,C13iC14»Clb»C16»C17,Cl8,C19»C20* 

8C21 » C22 * C2 3 , C24 . C25 » C26» C2 7, C2B « C29 . C30 «SAV«SAV1 f 

9SA V2 » SAV3 » SA V4» SAVE iSAV6»SAV7tSAV8.SAV9. SCI > THCPR » ENEWST » 

XENEWCT»EKP.EKP1 * EKP2 » E<P3 , EKP4 . EKP5 . ENEWTH , THC 1 .ENEWM, 

XAA «KAC » NRET , SW2 »EN ( 9 ) » PRNOPT *HOLD» J »K »L * STOR (1000) 

COMMON LTEST 

WBC=WB 

K0UNT=0 

ERAS2=SMUB/CMUB 
EMUB = ATAN (ERAS2 > 

SC = SA 

9 THCP=THC 
CALLBODYX 
IF (CNTRL8 ) 2» 1 »2 

2 ERASF=CTHB*CMUB 
ERAS 1 =STHB^<-SMUB 
E-RASE = ERAS1 +ERASE 
ERAS2=SMUB/CMUB 

ERASE= ( ERAS2*ERAS1 )/ERASE 
YBC= ( YB+YC )/DEC2 
ERASE= ( (XC-XB )/YBC )*ERASE 
G0T03 
1 ERASE=0 

3 CALLSL I TET ( 2 «LTEST ) 

GOTO (5»4 ) »LTEST 

5 CALLSLITE(2) 

G0T06 

4 AA=0 
G0T07 

6 I F (CNTRL9 ) 8, 7«8 

8 AA= ( { ( (SC-SB )/CNTRL3 )*SMUB )/CNTRL5 )*SMUB 

7 ERASE=AA-ERASE 

WC= ( (THB-THC)*ERAS2-ERAS£ )*WBC+WB 
CALLMUSR 
I F (DEC ) 21 » 20 i2 1 
20 EMUB= (EMUC + EMUB )/DEC2 
SMUB = S I N (EMUB ) 

CMUB = COS (EMUB ) 
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THBC= ( THB+THC ) /DEC2 
STHB=SIN(THBC ) 

CTHB = COS (THnC ) 

WBC= ( WB+WC )/DEC2 
KOUNT= ( KOUNT+KON ) 

1 F (K0UNT-KNTRL7 )9 « 1 0 « 1 0 

10 I F ( ABS < THC-THCP )-CNTRL6 ) I 1 « 1 1 

11 STHC=5IN(THC) 

CTHC=COS (THC ) 

21 RETURN 


END 

TIBFTC THASR REF 

SUBROUT INETHASR 

COMMONXA , YA . THA » WA , SMU A » CMUA 1 STHA , CTHA , SA , Xi3 » YB » THB * 

1 WB f SMUB 1 CMUB , STH6 » CTHB , SB , EPB » TEPB < XC » YC < THC « WC « 

2SMUC.CMUC,STHC,CTHC,SC,EPC.TEPC,£MUA,£MUB,EMUC,THAC, 

3WAC . SMUAC » YAC * EMU AC , STHAC , EL AC » THBC » WUC » SMUBC ’ ’ EMUBC « 

4STHBC « EPBC , TtPBC , £MdC » TANMA « TANMB , £MC » PC « CPC ♦ EM ♦ THCP » 

5COS1 ,TAN1 ,C0S2,TAN2.K0UNT, TAPE, ALIGHT, BEOS, BVC0N,SVC0N1 , 

6SVCON2,SVCON3,5VCON4,SVCON5,SVCON6,SVCON7,SVCON6,SVCON9,SVC010, 

7SVC01 1 , EEP,EEP1 , EEP2 , EEP3 , EEP4 , EEP5 , LEP6 , ElP7 , CNTRl , CNT RL 1 , 
8CNTRL2,CNTRL3, CNTRL4 ,CNTRLB,CNTRL6,KNTRL7, CNTRLt) ,CNTRL9,CN R 0, 

9KNTR1 1 ,CNTR12,CNTR13,CNTR14,CNTR15,CNTR16,CNTR17,CNTR18,CNTR19, 
XCNTR54, ERASE, ERAS 1 ,ERAS2, ERA S3, ERAS4, ERA Sb,ERAS6,ERAS7, 

XERASB ,ERAS9 ,ERAS1 0 ,ERAS1 1 , tR A S 1 2 , ERA S 1 3 , ERAS 1 4 , ER AS 1 U , n.R A S 1 6 
COMMONERAS17,ERAS18,ERAS19,ERAS20,ERAS21 , ERA S22 , ERASp , 

1 ERAS24 , ERAS25 ,ERAS26 ,ERAS27 , ERAS28 , eRAS 29 , ERAS30 , ERAS31 , 
2ERAS32,ERAS33,ERAS34 , ERAS35, DEC, UE Cl ,UEC2,0EC3,0EC4, 
3D£C5,DEC6,DEC7 ,uEC 8,DEC9,DEC10,DEC1 1 ,KON,CONl , PRT , 

4PRT1 ,PRT2,PRT3 ,PRT4 ,PRT5,PRT6,PRT7,PRT8, PRT 9 , DEG , 
g-pEMP, TEMPI ,T£MP2, TEMPS, T£ MPA, XI ,YI ,TH1 ,'iA/I ,SMUl , 

6CMU I , STH I , CTH 1 , SI ,KER1 0 , EX 1 , EX2 ,CI ,C2,C3,C4,C5,C6, 
7C7,C8,C9,C10,C1 1 ,C12,C13,C14,C15,Cl&,C17,Cia,C19,C20, 

8C2 1 ,C22,C23,C24,C25,C26,C27,C28,C29,C30,SAV,SAV1 , 

9SAV2 , SAV3 , SAVA , SA VB , SA V6 , SAV7 , SA V8 , SA V9 , SC 1 , TriCPR , lNEWST , 
XENEWCT , EKP , EKP 1 , EKP2 , EKP3 , EKP4 , £KP5 , ENEWTH , THC 1 , ENEWM , 

XAA,KAC,NRET,SW2,EN(9 ) ,PRNOPT,HOLU, J,K,L,STOR C 1000 ) 

g-RAS6 = EKP5*ERAS£ 

PrAS 7=CNTRL3-DFC 1 
frRASa = ERAS7+0£C2 
ErAS9=S IN (ENEWTH ) 

ERAS10=C0S (ENEWTH ) 
eras 1 1 =£RAS9/ERAS1 0 


43 



APPENDIX 


ERAS 1 2 = ERA36-DEC 1 
ERAS12=EKP5/ERAS1 2 
ERASI 3=DEC2*CNTRL3 
ERAS13=ERAS7/ERAS1 3 
ERAS14=ERAS2*ERAS2 
ERAS 1 5 = ERAS1 /ERAS2 
ERAS16=ERAS1*ERAS10 
ERAS16=-ERAS2*ERAS9+ERAS16 
ERASI 7 = ERAS16*ERAS1 6 
ERAS 1 8 = ERAS2*ERAS1 0 
ERASI 8 = ERAS1 *ERAS9 + ERAS 1 8 
ERAS19 = ERAS18*ERAS1 8 
ERAS20 = ERAS 16/ERAS 18 
ERAS21 =ERAS6-ERAS1 3 
ERAS21 =DEC1/ERAS21 
ERAS22=D£C 1/CNTRL3 

ERAS22= (ERAS6*ERAS1 3+ERAS22 ) *ERAS6 
FRAS23=-DEC1 /ERAS22 

ERAS21 = ( (ERAS23 + ERAS21 ) /ER A S 7 ) *DEC2 •*S'CNTRL5*Ei<P5*ERA S 1 *ER AS2 

ERAS23= ( (ERAS8/DEC2 ) *ERASl 2-DEC 1 ) /ERAS 1 4 

FRAS23=-FRAS 1 2*ERAS 1 2*ER AS8*ERASE+ERAS2 3 

CALLSSR 

SC=SC+SC1 

ERAS22=ERAS1 1*ERAS1 1 
ERAS24=-DEC1 /ERAS22 
ERAS24= (ERAS24-UEC1 ) /ERAS23 
ERAS5=ERAS24*ERAS21 

ERAS?6= ( (ERAS8/ERAS7)*ERAS15 1/ERAS20 
ERAS26= (ERAS26-UEC1 ) *ERAS1 7+DECl 
ERAS28=ERAS1 6*ERAS 1 8 
FRAS27= ( ERAS 28/ERAS 1 4 ) •»ERAS24 
ERAS29=ERAS24-DEC 1 

ERAS27= ( ( (ERASI 9-ER A S 1 7 ) *FRA S 1 5*ERAS29 + ERAS2 7 ) /DEC2 )*ERAS8 

ERAS27=-ERAS29*ERAS28*ERA57+ERAS27 

ERAS35=DEC1 /ERAS26 

ERAS3 = SQRT ( ERAS35 ) 

ERAS4= ( (-ERAS35*ERAS3)/ERAS7 )*ERAS27 

RETURN 

END 

SIBFTC MUSR REF 

SUBROUT I NEMUSR 

COMMON/ A ,YA,THA*WA,SMUA,CMUA,STHA«CTHA»SA,XB«YB»THB» 

1 WB»SMUB,CMUB*STHaiCTHB,SB,EPB»TEPB*XC»YC,THC»WC« 

2SMUC ,CMUC ,STHC , CTHC * SC ,EPC » TEPC « EMUA , EMUB , EMUC , THAC » 

3WAC , SMUAC , YAC »EMUAC » S TH AC , EL AC » T HBC , WbC i SMUBC » YBCi EMUBC , 
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4STHBC»EPBC,TEPBC»EME3C*TANMA.TANMBiEMCiPCiCPC»EM»THCP, 

5C0S1 ,TAN1 ,C0S2»TAN2»<0UNT»TAPE»ALIGHT«BEGS«SVC0NiSVC 0N1 « 
6SVCON2«SVCON3*SVCON4»SVC0N5iSVCON6»SVCON7iSVCON8,SVCON9,SVCO10 
7SVC01 1 . EEP,EEP1 , EEP2 » EEP3 * EEP4 » ££P5 » EEP6 « EEP7 i CNTRL « CNTRL I i 
8CN'TRL2»CNTRL3.CNTRL4 * CNTRL5 » CNTRL6 » KNTRL7 . CNTRL8 , CNTRL9 , CN T'R 1 0 
9KNTR1 1 iCNTR12.CNTR13*CNTR14,CNTR15.CNTR16.CNTR17,CNTR18,CNTRl9 
XCNTR54» ERASE* ERAS 1 , ERA S2 * ERA S3 * ERAS4 * BRASS * ERAS6 * ERAS7 » 
XERAS8*ERAS9.ERAS10*ERAS1 1 . ER AS 1 2 . ERAS 1 3 * ERAS 1 4 , ER AS 1 b * ERAS 1 6 
COMMONERASI 7*ERAS18*ERAS1 9*ERAS20*ERAS21 . ERAS22 . ERAS23 . 

1ERAS24 ,ERAS25*ERAS26 *ERAS27, ERAS28*ERAS29*ERAS30 *ERAS31 * 
2ERAS32 .ERAS33.ERAS34 , ERAS35 * DEC , DEC 1 * DEC2 . DEC3 * DEC4 , 
3DEC5*DEC6*DEC7,DEC8*DEC9.D£C10*DEC1 1 *KON*CONl *PRT* 

4PRT1 ,PRT2.PRT3*PRT4*PRT5*PRT6«PRT7*PRT8.PRT9*DEG* 

5TEMP, TEMPI , TEMP2 * TEMP3 * TEMP4 * X I * YI *THI *WI , SMU I * 

6CMUI *STHI ,CTHI .SI .KERIO.EXI * £X2 , C 1 * C2 . C3 * C4 * C5 . C6 * 

7C7.C8.C9 .CIO. Cl 1 .C12.C13.C14.Clb.C16.C17.C18.C19.C20. 

6C2 1 .C22 . C23 . C24 .C2S . C2 6 . C27 . C28 . C29 .C30 . SAV . SAVl . 

9SAV2 . SAV3. SAV4 . SAVE . SAV6 . SAV7 . SAV8. SAV9 . SCI . THCPR . ENEWST . 
XENEWCT.EKP.EKPl .EKP2 . EKP3 . EKP4 . £i<P5 . ENE WT H . THC 1 .ENtWM . 

XAA .KAC . NR£T . SW2 .EN ( 9 ) .F’RNOPT .HOLD. J . I< . L . STOR ( 1 000 ) 

ERASE=WC*WC 
ERASE=DEC1 /ERASE 
ERAS 1 =CNTRL3-DEC 1 

SMUC=( (ERASE-DECl )/DEC2)*ERASl 
IF(DEC1-SMUC)3610. 1 . 1 
1 CMUC = SQRT (DECl-SMUC ) 

SMUC = SQRT (SMUC ) 

ERASE=SMUC/CMUC 
EMUC=ATAN (ERASE ) 

RETURN 
3610 DEC=1 . 

RETURN 

END 

SI5FTC IPT REF 

subrout 1 NEI PT 

COMMONXA . YA . THA . WA .SMUA .CMUA . S T HA . C THA . S A . X3 , YB . THB . 

1 WB .SMUB . CMUB .STHB.CTHb.SB.EPb.TEPB.XC.YC.THC.WC. 

2SMUC . CMUC . STHC . CTHC . SC .EPC . TEPC . EMU A . EMUB . EMUC . THAC . 

3WAC . SMU AC . YAC . EMUAC . STHAC . ELAC. THBC . WuC . SMUBC . YoC . EMUBC . 
4STHBC.EPBC.TEPBC , EMBC . TANM A . TANMB . EMC . PC . CPC . EM . THCP . 

5C0S1 .TANl . C0S2.TAN2.K0UNT. TAPE. ALIGHT. BEGS. 3 VCON.SVCONl . 
6SVC0N2 . SVC0N3 . SVC0N4 . SVC0N5 . SVC0N6 . SVC0N7 . SVC0N8 . SVC0N9 . SVCO 1 0 
7SVC01 1 . EEP.EEPl . EEP2 . EEP3 . EEP4 . EEPb . EEP6 . EEP7 . CNTRL . CNTRL 1 . 
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8CNTRL2 » CNTRL3 » CNTRL4 tCNTRLS . CNTRL6 « KNTRL7 . CNTRL6 » CNTRL9 , CNTR 1 
9KNTR1 1 «CNTR12»CNTR13»CNTR14,CNTR15»CNTR16»CNTR17»CNTR18»CNTR1 
XCNTR54 ♦ ERASED ERAS 1 » ERAS2 * ERA S3 » ERAS4 i ERAS3 , ERAS6 » ERAS7 . 
XERAS8, ERAS9, ERAS 10. ERAS 1 1 . ER AS 1 2 . ERASl 3 . ERAS 1 4 , ER AS 1 5 , ERAS 1 6 
COMMONERAS 1 7 . ERASl 8 . ERAS 1 9 . ERAS20 . ERAS2 1 . ERAS22 . ERAS2 3 . 

1 ERAS24 .ERAS25.ERAS26.ERAS27.ERAS28.ERAS29.ERAS30.ERAS31 . 
2ERAS32.ERAS33.ERAS34 . ERAS35 . DEC . l/EC 1 . DEC2 . DEC3 , DEC4 . 
3DEC5.DEC6.DEC7.DEC8.DEC9.DEC10.DEC1 1 . KON , CON f . PRT . 

4PRT1 .PRT2.PRT3.PRT4 .PR To. PRT6.PRT7.PRT8.PRT9, DEG. 

5TEMP , TEMPI . TEMP2 , TEMP3 . TEMP4 . X 1 . Y I , TH I , W I . SMU I . 

6CMUI .STHI .CTHI ,SI .KERIO.EXI .EX2.C1 .C2.C3.C4.C5.C6. 

7C7.C8.C9. CIO. Cl 1.C12.C13.C14.C15.C16.C17.C18.C19.C20. 

8C21 .C22.C23.C24.C25.C26.C27.C28.C29.C30.SAV.SAV1 . 
9SAV2.SAV3.SAV4 . S A V5 . SA V6 . S A V7 , S A V8 . S AV9 . SC 1 . THCPR . ENEWST . 
XENEWCT.EKP.EKPl .EKP2 . EKP3 . EKP4 . EKP5 . ENEWTH , THC 1 .ENEWM. 

XAA .KAC .NRET .SW2 .EN ( 9 ) . PRNOPT . HOLD . J . I< . L . STOR (1000) 

IF (XC-XA-,01 )2 , 1 . 1 
2 NRET=4 
RETURN 
1 XI=XA+.01 

TEMP=XC-XA 
TEMP=. 01/TEMP 
Y I = ( YC-YA ) *TEMP+YA 
THI = (THC-THA )*TcMP+THA 
W I = ( WC-WA ) *TEMP+WA 
SI = (SC-SA )*T£MP+SA 
CTHI=COS(THI ) 

STH I =S I N ( TH I ) 

TEMP1=W1*WI 
TEMPI =DEC1 /TEMPI 
TEMP2=CNTRL3-DEC1 

SMUI = ( ( TEMPI -DEC 1 )/DEC2 )*TEMP2 
CMUI=SQRT (DECl-SMUI ) 

SMUI =SQRT (SMUI ) 

NRET=1 

RETURN 

END 

SIRFTC SHKPT REF 

SUBROUT I NESHKPT 

COMMONXA.YA.THA.WA.SMUA.CMUA.STHA.CTHA.SA.XB. YB.THB. 

1 WB.SMU6.CMUB.STHB.CTHir..Sb.EPB.TEPB,XC.YC.THC.WC. 
2SMUC.CMUC.STHC.CTHC.SC.EPC.TEPC.EMUA.EMUB.EMUC.THAC. 

3WAC . SMU AC . YAC . EMUAC . STHAC . EL AC. THBC . WBC . SMUBC . YoC . EMUBC . 
4STHBC.EPBC.TEPBC.EMBC.TANMA.TANMB.EMC.PC.CPC.EM.THCP. 

5C0S1 .TANl . C0S2. T AN2 . KOUNT, TAPE. alight. BEGS. 3 VCON.SVCONl . 
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65VC0N2 « SVC0N3 . SVC0N4 » SVC0N5 » SVCON6 « S VC0N7 , SVC0N8 , SVC0N9 , 8 VCO 1 0 
7SVC01 1 , EEPtEEPl , EEP2 « EEP3 « EEP4 » EEP5 ♦ EEP6 « EEP7 , CNTPL . CNTRL 1 » 
8CNTRL2» CNTRL3«CNTRL4 iCNTRLb « CNTRL6 » KNTRL7 » CNTRL8 » CNTRL9 , CNTR 1 0 
9KNTR1 1 «CNTR12*CNTR13«CNTR14.CNTR15.CNTR16»CNTR17,CNTR18,CNTR19 
XCNTR54 < ERASE t ERA SI « ERA S2 « ERA S3 » ERAS4 « ERA S5 , ER AS6 » ERAS7 » 
XERAS8,ERAS9,ERAS10,ERAS1 1 ♦ ER AS I 2 , ERAS 1 3 . ERAS 1 4 , ER AS 1 b , ER AS 1 6 
COMMONERAS17,ERA.<;i8,ERAS19,ERAS20»ERAS2 1 * ERAS22 t ERAS23 * 

1 ERAS24 «ERAS25 ,ERAS26 »ERAS27«£RAS28,ERAS29.cRAS30 «ERAS31 « 
2ERAS32»ERAS33»ERAS34 . ERAS35 * DEC » OEC 1 « 0EC2 . DEC3 , DEC4 , 
3DEC5»DEC6»DEC7«DEC8.DEC9.DEC10»D£C1 1 .KON,CONl «PRT» 

4PRT1 »PRT2.PRT3.PRT4»PRT5«PRT6.PRT7.PRT8,PRT9,DEG« 

STEMPt TEMPI « TEMP2, TEMPS, TEMP4 , X I , Y I , TH I , W I , SMU I , 

6CMUI ,STHI ,CTHI ,S1 ,KER10,£X1 , EX2 , C 1 , C2 , C3 , C4 , C5 , C6 , 
7C7,C8,C9,C10,C1 1,C12,C13,C14,C15,C16,C17,C18,C19,C20, 

8C2 1 ,C22,C23,C24 , C25 , C26 , C27 , C28 , C29 , C30 , SAV , SA VI , 
9SAV2,SAV3,SAV4,SAV5,SAV6,SAV7,SAV8,SAV9,SC1 , THCPR , ENE WST , 
XENEWCT,EKP,EKP1 , EKP2 , EKP3 , EKP4 , EKP5 , ENEWTH , THC 1 ,ENEWM, 

XAA ,KAC , NRET , SW2 ,FN ( 9 ) , PRNOPT ,HOLO , J , K ,L , STOR ( 1 000 ) 

CALLNEWB 
HOLO=- 1 . 

THAC=THA 

WAC=WA 

ERASE=SMUA/CMUA 
EMUA=ATAN (ERASE ) 

EMUAC=EMUA 

EP6C=EPB 

TEPBC=TEPB 

K0UNT=0 

ERASE=ENEWST/ENtWCT 
CALLEPSR 
CALLTHASR 
6 ERAS6=SIN(THAC) 

ERAS7=C0S ( THAC ) 

ERAS8=S IN (EMUAC ) 

ERAS9=C0S (EMUAC ) 

ERAS10 = COS(EPBC > 

ERAS 10= (SIN(EPBC) ) /ERAS 10 
ERASl 1 =ERAS8/ERAS9 
ERASl 3=ERAS6*ERAS8 
ERASl 3=ERAS7*ERAS9-ERAS1 3 
ERAS12=ERAS6*ERAS9 
ERASl 2=ERAS7*ERAS8+ERAS1 2 
ERAS 1 4 = ERAS1 2/ERAS 1 3 



APPENDIX 


XC = ERAS1 4-ERAS 1 0 
YC=-XB*ERAS1 0+YB 
XC= (XAi«-ERASl 4-YA + YC ) /XC 
YC=XC*ERAS1 0+YC 
YAC= ( YC+YA )/DEC2 
IE(CNTRL8) 1 ,2* 1 

2 ELAC=0 
G0T031 

1 ELAC=( ( (ERAS8*ERAS8 ) /ERAS9 )*ERAS6 5/ERASl 3 
ELAC=( ( XC-XA )/YAC )*ELAC 
31 ERAS15=ERAS4/WAC 

ERAS 1 5 = ERAS 1 5-ERASl 1 
ERAS 16= (ERAS3-WA )/WAC 
ERASl 6=ELAC-ERAS16 
IF(CNTRL9)4«3»4 

4 ERAS 1 8 = CNTRL3*CNTRL=5 
ERAS18= (ERAS8*ERAS8 ) /ERASl 8 
ERASl 5=ERAS1 8*ERAS5+ERAS1 B 
ERAS16= (SA-SC)*cRAS18+ERAS16 

3 ERAS15=ERAS16/ERAS15 
THCP=THC 
THC=THA+ERAS 1 5 
THCPR=THC+THC1 
WC=ERAS4*ERAS1 5+ERAS3 
CALLMUSR 

IF(DEC )21 «20»21 
20 EiV|UAC= ( EMUC + EMUA ) /DEC2 
WAC= ( WA+WC ) /DEC2 
THAC= (THA+THC )/UEC2 
ERAS7=C0S ( THCPR ) 

ERAS6=S IN ( THCPR ) 

ERASE =ERAS6/ERAS7 
CALLEPSR 

EPBC= (EPC+EPB-THCl 1/DEC2 
K0UNT = K0UNT + l<0N 
IF (K0UNT-KNTRL7 )6 »5 » 5 

5 IF(A6S(THC-THCP)-CNTRL6)7,7,6 

7 IF (CNTRL9)9,a,9 
9 CALLSSR 

GOTOlO 

8 SC = SA 
GOTOl 1 

10 SC=SC+SC1 

11 STHC=SIN(THC) 

CTHC = COS ( THC ) 
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EPC=EPC-THC1 
SAV = COS (EPC ) 

TEPC = S I N (EPC )/SAv 
HOLD= 1 . 

21 RETURN 
END 

SIBFTC SODYX REF 

SUBROUT I NEBODYX 

COMMONXA » YA « THA < WA « SMUA » CMUA * STHA * CTHA » SA ♦ XB » YB » THB « 
lWB*SMU3»CMUB«STHB*CTHd*S3* EPB »TEP3»XC»YC*THC«WC» 
2 SMUC«CMUC»STHC.CTHCiSC»EPC*TEPC.EMUA»EMU3»EMUC.THAC* 

3WAC , SMUAC « YAC <EMUAC » STHAC , ELAC< TH6C « WBC » SMUBC * YoC « EMUBC » 

4STHBC * EP3C ♦ TEPBC «EM3C * TANMA » TANMB » EMC ♦ PC * CPC « EM » THCP « 

5C0S1 ,TAN1 , C0S2 ♦ T AN2»'kOUNT,T A PE « Alight* BEGS «SVCON,SVCONl « 

6SVC0N2 » SVC0N3, SVC0N4 ♦ S VC0N5 * SVC0N6 * S VC0N7 ♦ S VCONB , S VC0N9 , SVCO 1 0 * 
7SVC01 1 * EEP*EEP1 . EEP2 * EEP3 * EEP4 * EEP5 ♦ EEP6 * EEP7 * CNTHL * CNTRL 1 * 
eCNTRL2.CNTRL3«CNTRL4 , CNTHL5 * CNTRL6 * i<NTRL7 . CNTHL8 * CNTRL9 * CNTR 1 0 ♦ 
9KNTR1 1 «CNTR12*CNTR13 , CNTR 1 4 * CNTR 1 5 * CNTH 1 6 * CNTH 1 7 * CNTR 1 8 * CNTR 1 9* 
XCNTR54* ERASE* ERAS 1 * ERAS2 ♦ ERA S3 » ERAS4 * BRASS * ERAS6 * ERAS7 * 

XERAS8 *ERAS9*ERAS10*ERAS11«ERAS12*ERAS13*ERA514*ERAS15*ERAS16 
COMMONERAS 1 7*ERAS18*ERAS19«ERAS20*ERAS21*ERAS22*£RAS23* 

1 ERAS24 * ERAS2 5 *ERAS26 *ERAS27 * ERAS28 * ERAS29 *ERAS30*ERAS31 * 

2ERAS32 * ERAS33 *ERAS34 *ERAS35 * DEC* DECl *DEC2 * DEC3*DEC4 * 
3DEC5*DEC6*DEC7*DEC8*DEC9*DEC10«DEC1 1 *KON*CONl *PRT. 

4PRT1 *PRT2*PRT3*PRT4*PRT5*PRT6*PRT7*PRT8.PRT9*DEG* 

5 TEMP * TEMPI * TEMP2 * TEMPS * TEMP4 *XI*YI*THI*WI*SMUl* 

6CMU1 *STHI *CTHI *SI *KER10*EX1 * EX2 * C 1 * C2 * C3 * C4 * C5 * C6 * 
7C7*C8*C9*C10*C1 1*C12*C13*C14*C15*C16*C17*C18*C19*C20» 

8C21 *C22*C23*C24*C25*C26*C27*C2d*C29*C30»SAV«SAVl * 
9SAV2*SAV3*SAV4*SAV5*SAV6*SAV7*SAV8*SAV9*SC1*THCPR* ENEwST * 
XENEWCT * EKP * EKP 1 * EKP2 * EKP3 * EKP4 * EKP5 * ENEWT H * T HC 1 *ENEWM* 

XAA *KAC *NRET * SW2 *EN ( 9 ) * PRNOPT *HOLD * J * K * L * STOR ( 1 00 0 ) 

YC = 0 
THC=0 

ERASE=SMUB/CMUB 
ERAS 1 =STHB/CTHB 
ERAS2=ERAS 1 *ERASE+DEC 1 
ERASE= ( ERAS 1 -ERASE ) /ERAS2 
XC= (ERASE*XB-YB )/ERASE 
RETURN 

C CENTERL I NEBODY 

END 

SIBFTC SUBR REF 
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SUBROUT I NESUBR 

COMMONXA. YA.THA»WA,SMUA»CMUA » STHA , CTHA , SA , XB ♦ YB « THB , 

1 WB t SMUB « CMUB f STHB « CTHti * SB , EPB « TEPB < XC » YC » THC i WC « 

2SMUC « CMUC » ST HC » C THC * SC , EPC • TEPC f EMUA « EMUB « EMUC » THAC » 

3 WAC « SMUAC . YAC « EMUAC « STHAC , EL AC « TH8C * WbC . SMUBC » YBC » EMUBC « 

4STHBC » EPBC « TEPBC , EMBC « TANMA » TANMB « EMC « PC ♦ CPC * EM * THCP » 

5C0S1 ,TAN1 ,C0S2,TAN2»K0UNT,TAPE,ALIGHT<B£GS«SVC0N,SVC0N1 « 
6SVC0N2«SVC0N3«SVC0N4»SVC0N5f SVCON6«SVCON7,SVCON8,SVCON9»SVCOlO» 
7SVC01 1 ♦ EEP«EEP1 » EEP2 * EEP3 « EEP4 * EEP5 » EEP6 « EEP7 » CNTRL » CNTRL 1 . 
8CNTRL2» CNTRL3»CNTRL4 . CNTRL5 * CNTRL6 * KNTRL7 ♦ CNTHL8 « CNTRL9 * CNTR 1 0 ♦ 
9KNTR1 1 «CNTR12«CNTR13»CNTR14,CNTR15»CNTR16.CNTR1 7 » CNTR 1 8 i CNTR 1 9 » 
XCNTR54 » ERASE « ERAS 1 » ERAS2 * ERAS3 *ERAS4 » ERAS5 » ERAS6 . ERAS7 « 
XERAS8»ERAS9*ERAS10,ERAS1 1 ♦ ER AS,1 2 » ER AS 1 3 » ERAS 1 4 , ER AS 1 b , ERAS 1 6 
C0MM0NERAS17»ERAS18.ERAS19.ERAS20.ERAS21 ♦ ERA S22 ♦ ERAS23 « 

1 ERAS24 *ERAS25 ,ERAS26.ERAS2 7,ERAS28fERAS29,ERAS30«ERAS3l « 

2ERAS32 *ERAS33.ERAS34 ,ERAS35» OEC« OECl « DEC2 » DEC3 , 0EC4 , 
3DEC5.DEC6»DEC7fDEC8*DEC9»DEC10«DECl 1 tl<ON,CONl «PRT» 

4PRT1 ,PRT2 «PRT3 »PRT4 ♦PRT5«PRT6«PRT7,PRT8»PRT9»0EG ♦ 

5TEMP, TEMPI , TEMP2 , TEMPS « TEMP4 , X I « YI ,THI .WI ♦SMUI « 

6CMUI ,STHI ,CTHI ,SI .KERlOtEXl » EX2 « C 1 , C2 » C3 , C4 « C5 « C6 « 
7C7.C8«C9,C10«C1 1 « C 1 2 . C 1 3 » C 1 4 , C 1 5 . C 1 6 . C 1 7 , C 1 8 , Cl 9 . C20 » 

8C2 1 . C22 ♦ C23 . C24 » C25 . C26 , C2 7 , C2B » C29 , C30 ♦ S AV « SA VI ♦ 
9SAV2«SAV3«SAV4,SAV5*SAV6»SAV7»SAV8,SAV9*SC1 , THCPR ♦ ENEWST » 
XENEWCT.EKP«EKP1 « EKP2 . E<P3 , EKP4 » E<P5 ♦ ENE WTH , THC 1 .ENEWM, 

XAA »KAC * NRET » SW2 *EN ( 9 ) «PRNOPT .HOLD ♦ J , K »L . STOR (1000) 

ENEWM=2 .685947+. 725327*YC+SQRT (-1 7. 506468* YC*YC 
1 +7. 870908* YC-. 1 52477 )-( 1 .6~2.8o7*YC ) * ( XC- ( YC+ . 74 79 ) / . 277 ) 

THCl =- .0325/ ( YC**.647b ) -1 . 65*YC**2 . 795+ . 0424* ( XC- 
1 (YC+.7479)/.277)/YC**l .068 
THCl =-THCl 
SCI =153. /YC**. 675 


RETURN 

END 

SISFTC MOVIC REF 

SUBROUTINEMOVI C 
C 

THESE SUBROUTINES USED I IV BOTH FIRST AND SECOND SHOCK REGIONS 


COMMONXA , YA , THA . WA , SMUA. CMUA , STHA , CTHA , SA . XB « VB . THB ♦ 

1 WB. SMUB. CMUB. STHB. CTHB, SB. EPB .TEPB. XC.YC. THC. WC. 

2SMUC . CMUC . STHC . CTHC . SC . EPC . TEPC . EMUA . EMUB . EMUC . THAC . 

3 WAC . SMUAC . YAC .EMUAC . STHAC .EL AC . THBC . WBC . SMuBC . YBC .EMUBC . 

4 STH3C. EPBC. TEPBC .EMBC. TANMA . TANMB . EMC . PC . CPC . EM . THCP . 

5C0S1 .TANl .C0S2.TAN2.K0UNT. TAPE. ALIGHT. BEGS. SVCON.SVCONl . 

6SVC0N2 . SVC0N3 . SVC0N4 . SVC0N5 . SVC0N6 . S VC0N7 . SVC0N8 . S VC0N9 . S VCO 1 0 . 
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7SVC01 1 » EEPfEEPl » EEP2 » EEP3 * EEP4 . EEP5 » EEP6 » EEP7 » CNTPL i CNTRL 1 » 
8CNTRL2»CNTRL3»CNTRL4 . CNTRL b . CNTRL6 i KNTRL7 » CNTRL8 , C.NTRL9 » CNTR 1 0 
9KNTR1 1 «CNTR12»CNTR13*CNTR14,CNTR15»CNTR16.CNTR17,CNTR18*CNTR19 
XCNTR54, ERASE* ERAS 1 . ERA S2 * ERA S3 * ERAS4 * ER A S5 , ER AS6 » ERAS7 * 
XERAS8«ERAS9«£RAS1 0 ♦ ERA S 1 1 * ER AS 1 2 * ERAS 1 3 * ERAS 1 4 * ER AS 1 5 , ER AS 1 6 
COMMONERASI 7 *ERAS1 8 *ERAS1 9 .ERAS20 *ERAS2 1 *ERAS22 i ERA^23, 

1 ERAS24 *ERAS25 .ERAS26 . £R A S27 . ERAS28 * ERAS29 * ERAS30 i ERAS31 * 
2ERAS32 , ERAS33 .ERAS34 *ERAS35 , DEC ♦ DECl * DEC2 , DEC3 ,DEC4 * 
3DECS,D£C6f DEC7*DEC8*DEC9.DEC10*DEC1 1 *KON*CONl *PRT* 

4PRT1 *PRT2*PRT3*PRT4,PRT5*PRT6*PRT7*PRT8.PRT9*DEG^ 

5TEMP . TEMPI * TEMP2 , TEMP3, TEMP4 . X I * Y I , TH I , W I , SMU I * 

6CMUI ,STHI ,CTHI ,SI ,i<ER10«EXl iEX2*Cl i C2 * C3 , C4 i C5 , C6 * 

7C7 *C8 *C9,C10*Cll.C12*C13*C14.C15*C16iC17*Cl8 . Cl 9, C20 * 

8C2 1 . C22 . C23 , C24 * C25 * C26 * C2 7 . C2B * C29 * C30 * S AV * SAV 1 * 
9SAV2.SAV3*SAV4.SAV5*SAV6*SAV7iSAV8,SAV9,SCl * THCPR i ENEWST * 
XENEWCT*EKP*EKP1 *EKP2 . EKP3 , E.<P4 i EKP5 ♦ ENEWTH , THC 1 lENEWM, 

XAA ,KAC .NRET * SW2 *EN ( 9 ) * PRNOPT *HOLD , J,l< *L * STOR (1000) 

XC = X I 
YC = Y I 
THC = TH I 
WC = W I 
SMUC=SMU I 
CMUC=CMU I 
STHC=STH I 
CTHC=CTH I 
SC = S I 
RETURN 
END 

SIBFTC Moves REF 

SUBROUT I NEMO VC S 

COMMONXA , YA,THA*WA,SMUA, CMUA *STHA,CTHA*SA,Xd,YB*THB* 

1 WB *SMUB .CMUB *STHB , CTHB ♦ SB , EPB , TEPB * XC * YC * THC * WC * 

2SMUC .CMUC *STHC *CTHC * SC *EPC i TEPC * EMUA .EMUB , EMUC . THAC * 

3WAC , SMUAC ,YAC,EMUAC,STHAC,ELAC*THBC«WBC*SMu 6C*YBC*EMUBC* 
4STHBC*EPBC«TEPDC ,EMBC* TANMA , TANMB,EMCiPC .CPC .EM.THCP, 

5C0S1 .TANl ,C0S2.TAN2,K0UNT, TAPE. ALIGHT. BE GS.SVCON.SVCONl . 
6SVCON2 . SVC0N3 . SVCON4 .SVC0N5 . SVC0N6 . S VC0N7 . SVCONB . S VC0N9 . SVCOl 0 
7SVC01 1 . EEP.EEPl .EEP2.EEP3.EEP4.EEP5.EEP6.EEP7.CNTRL.CNTRL1 . 
8CNTRL2.CNTRL3,CNTRL4.CNTRLb.CNTRL6.KNTRL7.CNTRL8.CNTRL9.CNTK10 
9KNTR 1 1 . CNTRl 2 . CNTR 1 3 . CNTR 1 4 . CNTR 1 5. CNTR 1 6 .CNTRl 7 . CNTRl 8 . CNTRl 9 
XCNTR54 .ERASE . ERASl , ERAS2 . ERA S3 . ERAS4 .ERAS5 .ERAS6 . ERAS7. 

XERAS8 . ERAS9 . ERAS 1 0 , ERAS 1 1 . ER AS 1 2 . ER A S 1 3 . E RAS 1 4 . ER A S 1 5 . ER AS 1 6 
C0MM0NERAS17 . ERA S 1 8 . ER AS 1 9 . E RAS20 . ER A.S2 1 .ERAS22 .ERAS23. 
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1ERAS24 ,ERAS25,ERAS26»ERAS27,ERAS28f ERAS29»ERAS30 *ERAS31 , 

2ERAS32.ERAS33»ERAS34 f ERAS35 « DEC » DEC 1 t DEC2 « DEC3 * DEC4 » 

3DEC5.DEC6, DECT, DECS. DEC9,DEC10»DEC1 1 .KON.CONl .PRT. 

4PST1 ,PRT2,PRT3,PRT4 .PRT5.PRT6.PRT7.PRT8.PRT9 .DEG. 

5TEMP . TEMP 1 . TEMP2 . TEMP3 . TEMP4 . X I . Y I . T.H I . W I . SMU I . 

6CMU1 .STHI .CTHl .SI .<ER10.EX1 . EX2 . C 1 .C2.C3.C4.C5.C6. 

7C7.C8.C9. CIO. Cl 1 . C 1 2 . C 1 3 . C 1 4 . C 1 5 . C 1 6 . C 1 7 . C 1 8 . C 1 9 . C20 . 

8C21 .C22.C23.C24 . C25 . C26 . C27 . C28 . C29 . C30 . SAV . SAV 1 . 

9SAV2.SAV3.SAV4.SAV5.SAV6.SAV7.SAV8.SAV9.se! . THCPR . ENEWST . 

XENEWCT.EKP.EKPl . EKP2 . £<P3 . EKP4 . EKP5 . ENE WTH . THC 1 .ENEWM. 

XAA .KAC . NRET .SW2 «EN ( 9 ) . PRNOPT . HOLD . J . K . L . STOR (1000) 

SAV=XC 

SAVl =YC 

SAV2=THC 

SAV3=WC 

SAV4=SMUC 

SAV5=CMUC 

SAV6=STHC 

SAV7=CTHC 

SAV8=SC 

RETURN 

END 

S.IBFTC MOVSC REF 

SUBROUT I NEMOVSC 

COMMONXA . YA . THA . WA . SMUA , CMUA . STHA . CTHA . SA . XB . Y3 . THB . 

1 WB .SMUB .CMU8 . STHB . CTHB . SB . EPB . TEPB . XC . YC . THC . WC . 

2SMUC.CMUC.STHC.CTHC.SC.EPC.TEPC.EMUA.EMUB.EMUC.THAC. 

3WAC . SMU AC . YAC .EMUAC . STHAC . EL AC . THBC . WbC . SMUBC . YBC . EMUBC . 

4STHBC.EPBC.TEPBC.EMBC.TANM A. TANMB. EMC. PC. CPC. EM. THCP. 

5C0S1 .TANl . C0S2. TAN2.K0 UNT. TAPE. Alight. BEGS. SVCON.SVCON l . 

6SVC0N2. SVC0N3. SVC0N4 . S VC0N5 . SVC0N6 . SVC0N7 . S VC0N8 . SVC0N9 . SVCO 1 0 . 

7SVC01 1 . EEP.EEPl . EEP2 . EEP3 . EEP4 . EEP5 . EEP6 . EEP7 . CNTRL . CNTRL 1 . 

8CNTRL2.CNTRL3.CNTRL4 .CNTRlB . CNTRL6 . KNTRL7 . CNTRL8 . CNTRL9 . CNTR 1 0 . 

9KNTR1 1 .CNTR12.CNTR13.CNTk’14.CNTR15.CNTR16.CNTR17.CNTR18.CNTK19. 

XCNTR54 .ERASE. ERAS 1 . ERAS2 . ERA S3 . ERAS4 . BRASS . ERAS6 . ERAST . 

XERAS8.ERAS9.ERAS1 0. ERAS I 1 . ERAS 1 2 . ERAS 1 3 . ERAS I 4 . ERAS I 5 . ERAS 1 6 
C0MM0NERAS17.ERAS18.ERAS19.ERAS20.ERAS2 1 . ERAS22 . ERAS23 . 

1 ERAS24 . ERAS25 . ERAS26 . ERAS2 7 . ERAS28 . ERAS29 . FRAS30 . ERAS31 . 

2ERAS32.ERAS33.ERAS34 . ERAS35 . DEC . DEC 1 . DEC2 . DEC3 . DEC4 . 

3DEC5.DFC6.0EC7.DEC8.DEC9.DEC10.DEC1 1 .KON.CONl .PRT. 

4PRT1 .PRT2.PRT3.PRT4.PRT5.PRT6.PRT7.PRT8.PRT9.DEG. 

5TEMP. TEMPI . TEMP2 . TEMP3 . TEMP4 . X I .YI .THI .Wl . SMU I . 

6CMU I . STH I . CTH 1 .SI .KERl 0. EX 1 . EX2 .Cl .C2.C3.C4.C5.C6. 

7C7.C8.C9.C1 0. Cl 1 . C 1 2 . C 1 3 . C 1 4 . C 1 5 . C 1 6 . C 1 7 . C 1 8 . C 1 9 . C20 . 

8C21 . C22 .C23. C24 .C25 . C26. C27. C28 . C29 .C30 . SAV . SAVl . 
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9SAV2*SAV3»SAV4.SAV5»SAV6.SAV7i3AV8»SAV9iSC1 « THCPR * ENE WST » 
XENEWCT«E<P»EKP1 . EKP2 » EKP3 , EKP4 . E<P5 » ENEWTH , THC 1 lENEWM, 

XAA ,KAC » NRET » SW2 » EN ( 9 ) » PRNOPT «HOLD ,JiK»Li5TOR(1000) 

XC=SAV 

YC=SAV1 

THC=SAV2 

WC=SAV3 

SMUC=SAV4 

CMUC=SAV5 

STHC=SAV6 

CTHC=SAV7 

SC=SAVB 

RETURN 

END 

SIdFTC MOVBA ref 

SUBROUT I NEMOVBA 

COMMONXA »YA»THAiWA.SMUA,CMUA»STHA,CTl-iA*SA*x3»YB»THB» 
1WB.SMUB«CMUB.STH8»CTHB,SB,EP6 iTEPB,XC«YC«THC.WC» 

2SMUC 1 CMUC 1 STHC ,CTHC » SC ,EPC » TEPC i EiYu A i EMUB , EMUC , THAC » 

3WAC , SMUAC »YAC,EMUAC«STHAC,ELACiTHBCi WBC » SMUBC * YBC » EMUBC < 

4STHBC ,EPBC»TEPBC,EMBC»TANMA,TANM3,EMC»PC» CPC . EM » THCP » 

5C0S1 ,TAN1 , C0S2 , TAN2 ,KOUNT , TAPE » AL I GHT » BEGS » SVCON , SVCONl , 
6SVC0N2« SVC0N3. SVCON4 » SVCONB » SVC0N6 * SVC0N7 » SVCONS . SVC0N9 » SVCOl 0 
7SVC01 1 . EEP.EEPl i EEP2 » EEP3 » EEP4 » EEP5 , EEP6 i EEP7 . CNTRL i CNTRL 1 » 
8CNTRL2 » CNTRL3 » CNTRL4 tCNTRLb » CNTRL 6 1 KNTRL7 * CNTRL6 , CNTRL9 » CNTRl 0 
9KNTR 1 1 ,CNTR12»CNTR13,CNTRl4.CNTRl5*CNTR16»CNTR17,CNTRia.CNTRl9 
XCNTR54 , ERASE < ERAS 1 » ERAS2 .ERAS3 » ERAS4 »ERAS5 * ERAS6 1 ERAS7 1 
XERAS8 . ERAS9 , ERAS 1 0 , ERA SI 1 , ERAS 1 2 , ERAS 1 3 , ERAS 1 4 ,ER ASl 5 , ERAS 1 6 
COMMONER AS 17«ERAS18»ERASl'9»ERAS2 0iERAS21»ERAS22«ERAS2 3i 
1 ERAS24 ,ERAS25 ,ERAS26 .ERAS27 , ER AS28 , ERAS29 , ER A S30 i ER ASS 1 i 
2ERAS32 ,ERAS33 ,ERAS34 »ERAS35 * DEC » UECl .DEC2 » DECS » DEC4 » 
3DEC5,DEC6»DEC7,DEC8,DEC9.DEC10*DEC1 1 ,i<ON,CONl »PRT* 

4PRT1 ,PRT2,PRT3,PRT4,PRT5.PRT6.PRT7,PRT8,PRt9,DEG« 

5TEMP, TEMPI , TEMP2 , TEMPS . TEMP4 « X Y I , T H I » W I i SMU I * 

6CMUI .STHI ,CTHI ,S1 ,KER1 OiEXl » EX2 » C 1 , C2 . C3 i C4 » C5 , C6 » 

7C7 »Ca»C9»C10»Cl 1 ,C12,C13*C14»ClbiC16»C17,Cl8,C19,C20» 

8C2 1 . C22 » C23 » C24 » C25 . C26 . C2 7 i C28 i C29 , C30 i SAV « SA VI » 
9SAV2«SAV3»SAV4»SAV=>»SAV6.SAV7,SAV8»SAV9.SC1 « THCPR * ENEWST » 
XENEWCTiEKP.EKPl ,EKP2 , EKP3 , EKP4 » EKP5 i ENEWTH , THC 1 «ENEWM, 

XAA iKAC» NRET, £‘a/ 2 »EN (9 ) . PRNOPT i HOLD i J i K . L » STOR ( 1000) 

XA = XB 
YA = YB 
THA=THB 
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WA = WB 

SMUA=SMUB 

CMUA=CMUB 

stha=sthb 

CTHA=CTHB 
SA = SB 
RETURN 
END 

SIBFTC MOVCB REF 

SUBROUT I NEMOVCB 

COMMONXA * YA , THA * WA , SMUA, CMUA ♦ STHA , CTHA « 3 A « XB , YB » THB ♦ 

1 WBiSMUBiCMUB.STHB»CTHB,SB,EPS.TEPB«XC.YC«THC»WC» 
2SMUCiCMUC,STHC,CTHC,SC ♦ EPC » TEPC ♦ EMUA , EMUB « EMUC « THAC . 

3WAC , SMUAC , YAC , EMUAC 1 STHAC » EL AC» THBC » WBC » SMUBC i YBC « EMUBC » 

4STHBC* EPBC 1 TEPBC « EMBC» TANMA « TANMB «EMC «PC « CPC « EM « THCP* 

BCOSl 1 TANl , C0S2 » TAN2 «KOUNT ♦ TAPE. AL I GHT . BEGS . S VCON . SVCONl . 

6SVC0N2 . SVC0N3. SVC0N4 .SVC0N5. SVCON6. SVC0N7.SVC0N8 .SVC0N9.SVC01 0. 
7SVC01 1 . EEP.EEPl . EEP2 . £EP3 » EEP4 . EEP5 . EEP6 . FEP7 . CNTRU » CNTRL 1 » 
8CNTRL2 . CNTRL3 .CNTRL4 , CNTRL5 . CNTRU6 . KNTRL7 . CNTRL6 . CNTRL9 ♦ CNTR 1 0 . 
9KNTR1 1 .CNTR12.CNTR13.CNTR14.CNTR15.CNTR16.CNTR17.CNTR18.CNTR19. 
XCNTR54 . ERASE . ERAS 1 . ERAS2 » ERA S3 « ERAS4 .ERAS5 » ERAS6 « ERAS7 « 
XERAS8.ERAS9.ERAS1 O.ERASl 1 , ERAS 1 2 . ERAS 1 3 . ERAS 1 4 , ER AS 1 5 . ERAS 1 6 
C0MM0MERAS17.ERAS18.ERAS19.ERAS20.ERAS21 ♦ ERAS22 . ERAS23 . 

1 ERAS24 .ERAS25 .ERAS26 iERAS27.£RAS23. ERAS29.ERAS30 .ERAS31 . 
2ERAS32.ERAS33.ERAS34 . ERAS35 . DEC . DEC 1 . DEC2 . DECS . DEC4 . 

3DEC5 , DEC6 . DEC7 . decs . DEC 9 . DEC 1 0 . DECl 1 .KON . CONI .PRT . 

4PRJ1 ,PRT2.PRT3.PRT4 .PRT5.PRT6.PRT7.PRT8.PRT9.DEG. 

5TEMP. TEMPI , TEMP2 . TEMPS . TEMP4 . X I .YI ,THI .Wl .SMUI . 

6CMUI .STHI ,CTHI ,SI .KERIO.EXI . EX2 .Cl . C2 . C3 . C4 . Co , C6 . 

7C7.ce.C9. CIO. Cl 1 .Cl2,C13.C14.C15,C16.C17.Clb,Cl9.C20. 

8C2 1 . C2? . C23 . C24 . C25. C26. C27, C28 . C29 . C30 . SAV . SAVl . 
9SAV2.SAV3.SAV4.SAV5.SAV&.SAV7.SAV8.SAV9.3C1. THCPR .ENEWST » 
XENEWCT . EKP . EKP 1 . EKP2 . EKP3 . EKP4 . EKP5 . ENEWTH , THC 1 . ENEWM . 

XAA .KAC . NRET . SW2 .EN ( 9 ) . PRNOPT . HOLD , J . K . L . STOR (1000) 

XB = XC 

YB = YC 

TH6=THC 

WB = WC 

SMUB=SMUC 

CMUB=CMUC 

STHB=STHC 

CTH3=CTHC 

SB = SC 

RETURN 

END 
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SIBFTC WRTSK REF 

SUBROUT I NEWRTSK 

COMMONXA » YAiTHA»WA,SMUA*CMUA.STHA»CTHA»SA,Xd»YB«THB» 
1 WB.SMUB.CMUB»STHB*CTHB,SB,EPB»TEPBiXC*YCiTHCiWC» 

2SMUC .CMUC . STHC . CTHC * SC «EPC t TEPC » EMU A lEMUB .EMUC » THAC » 
3WAC,SMUAC, YAC«EMUAC,STHAC,ELAC»THBC»WBCi SMUBC t YdC » EMUBC , 
4STHBC.EPBC.TEPBC.EMBC»TANMA,TANMB,EMC.PCiCPC»EM,THCP« 

5C0S1 ,TAN1 , C0S2. T AN2 1 KOUNT, TAPE. ALIGHT* BEGS iSVCON,SVCONl . 
6SVCON2.SVCON3.SVCON4.SVCON5.SVCON6.SVCON7.SVCON8,SVCON9,SVC010 
7SVC01 1 . EBP.EEPl . EEP2 . EEP3 . EEP4 . EEP5 . EEP6 . EEP7 , CNTRL . CNTRL 1 . 
8CNTRL2 . CNTRL3 .CNTRL4 .CNTRL5. CNTRL6 .KNTRL7 . CNTRL8 , CNTRL9 , CNTR 1 0 
9KNTR1 1 .CNTR12.CNTR13.CNTRl4,CNTR15.CNTR16.CNTR17.CNTRia,CNTR19 
XCNTR54. ERASE .ERAS 1 . ERA S2. ERA S3 . £R AS4 . ER AS5 . ER AS6 . ER AS7 . 
XERAS8.ERAS9.ERAS10.FRAsri . ER AS 1 2 . ER AS 1 3 . ERA S 1 4 . ER AS 1 5 . ER AS 1 6 
COMMONERAS 1 7.ERAS1S.ERAS19.ERAS20.ERAS21.ERAS22.ERAS23. 

1 ERAS24.ERAS2B.ERAS26.ERAS2 7.ERAS28.ERAS29.ERAS30 .ERAS31 . 
2ERAS32 . ERAS33 .ERAS34 .ERAS35 . DEC . OECl .0EC2 . DEC3 .0EC4 . 
3DEC5.DEC6.OEC7.DECa.DEC9.0EC10.DECl 1 .KON.CONl .PRT. 

4PRT1 .PRT2 .PRT3 .PRT4 . PRT 5 . PR T6 . PR T7 . PRT8 . PRT9 . DEG , 

5 TEMP . TEMPI . TEMP2 . TEMP3 . TEMP4 . X I . Y I . TH I . W I . SMU I . 

6CMU I .S TH I . CTH I .SI .KERl 0 . EX 1 . EX2 .Cl .C2.C3.C4.C5.C6. 
7C7.ca.C9.C10.Cl 1 .C12.C13.C14.C15.C16.C17. CIS. C19.C20. 

8C2 1 . C22 .C23 . C24 . C25 . C26. C27 . C28 . C29 . C30 . SAV . SAVl . 
9SAV2.SAV3.SAV4.SAV5.SAV6.SAV7.SAva.SAV9.SCl . THCPR . ENEWST . 
XENEWCT.EKP.EKPl .EKP2 . E KP3 . EKP4 . EKP5 . ENE WT H . T HC 1 .ENEWM. 

XAA .KAC . NRET . SW2 .EN ( 9 ) . PRNOPT .HOLD.J.K.L.STORI 1000) 

EXl =FPC 

EX2=TEPC 

CALLWRTGN 

RETURN 

END 

SIBFTC WRTGN REF 

SUBROUT I NEWRTGN 

COMMONXA . yA.THA.WA.SMUA.CMUA.STHA.CTHA.SA.XB.YB.THB. 

1 WB . SMUB . CMUB . STHB . CTHB . SB , EPS . TEPB . XC . YC . THC . WC . 

2SMUC . CMUC . STHC . CTHC . SC . EPC . TEPC .EMUA.EMUB.EMUC.THAC. 

3WAC .SMUAC . YAC .EMUAC . ST H AC. EL AC. THBC . WBC . SMUBC . YBC . EMUBC . 
4STHBC.EPBC. TEPBC .EMBC. TANMA . TANM6 .EMC .PC .CPC . EM.THCP. 

5C0S1 .TANl . CO S2. T AN2, KOUNT. TAPE. ALIGHT. BEGS. SVCON.SVCONl . 
6SVCON2 . SVCON3 . SVC0N4 . SVC0N5 . SVC0N6 . SVC0N7 . SVC0N8 . S VC0N9 . SVCOl 0 
7SVC01 1 . EEP.EEPl . EEP2 . EEP3 , EEP4 . EEP5 . EEP6 . EEP7 . CNTRL . CNTRL 1 . 
8CNTRL2 . CNTRL 3 . CNTRL4 . CNTRL5 . CNTRL6 . KNTRL7 . CNTRL8 . CNTRL9 . CNTR 1 0 
9KNTR1 1 .CNTR 12. CNTR 13. CNTR 14. CNTR 15. CNTR 16. CNTR 17. CNTR 18. CNTR 19 
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XCNTR54f ERASE* ERAS 1 ♦ ERASE * ERA S3* ERAS4 ♦ ERASE ,ERAS6 * ERAS? * 
XERAS8*ERAS9*ERAS1 0*ERAS1 1 * ERAS 12 * ERAS 1 3 * ERAS 1 4 * ERAS 1 5 * ERAS 1 6 
COMMONERAS17*ERAS18*ERAS19*ERAS20*ERAS21 *ERAS22*ERAS23* 
1ERAS24*ERAS25*ERAS26*ERAS27*ERAS28*ERAS29*ERAS30*ERAS31 ♦ 
2ERAS32*ERAS33*ERAS34*ERAS35*DEC*DEC1 * DEC2 * DECS * DEC4 * 
3DEC5*DEC6*DEC7*DEC8*DEC9*DEC10*DEC1 1 *K0N*C0N1 *PRT* 

4PRT1 *PRT2*PRT3*PRT4*PRT5*PRT6*PRT7*PRT8*PRT9*DEG* 

5TEMP* TEMPI * TEMP2 * TEMPS * TEMP4 * X I *YI *THI *WI * SMU I * 

6CMUI *STH I * CTH I ,S l.*KERl 0 * EX 1 * EX^ * Cl * C2 * C3 * C4 * C5 * C6 * 
7C7*'C8*C9*C10*C1 1 * C 1 2 * C 1 3 * C 1 4 * C I 5 * C 1 6 * C 1 7 * C 1 8 * C 1 9 * C20 * 

8C21 *C22*C23*C24*C25*C26*C27*C23*C29*C30*SAV*SAV1 * 

9SAV2 * SAV3*SAV4* SAVE *SAV6*SAV7*SAV8*SAV9* SCI * THCPR * ENEWST * 
XENEWCT*EKP*EKP1 *EKP2 * EKP3 , EKP4 * EKPS * ENEWTH * THC 1 *£NEWM* 

XAA *KAC*NRET * SW2 *EN (9 ) *PRNOPT *HOLD* J *K *L * STOR < 1000 ) 

STOR ( K ) =XC 
STOR (K+l )=YC 
STOR (K+2 ) =THC 
STOR (K+3 ) =WC 
ST0R(K+4)=SMUC 
STOR ( K+E ) =CMUC 
STOR (K+6 ) =STHC 
ST0R(K+7)=CTHC 
STOR (K + 8 ) =SC 
K = K + 9 
RETURN 
END 

SIBFTC SSR REF 

subrout I NESSR 

COMMONXA*YA * TH A * W A * SMUA * CMUA * STHA * CTHA * S A * XB * YB * THB * 

1 WB*SMUB*CMUB*STHB*CTHB* SB*EPB*TEPB*XC*YC*THC*WC* 

2SMUC * CMUC *STHC * CTHC * SC * EPC * TEPC * EMUA * EMUB , EMUC * THAC * 

3WAC*SMUAC* YAC*EMUAC*STHAC*ELAC*THBC*WBC*SMUBC* YBC*EMUBC* 

4STHBC * EPBC * TEPBC , EMBC* TANMA * TANMB * EMC * PC * CPC * EM * THCP* 

ECOSl *TA'N1 *C0S2*TAN2*K0UNT*TAPE*AL1GHT*BEGS*SVC0N*SVC0N1 * 
6SVCON2*SVCON3*SVCON4*SVCON3*SVCON6*SVCON7*SVCON8*SVCON9*SVCO10* 
7SVC01 1 * EEP*EEP1 * EEP2 * EEP3 * EEP4 * EEPE * EEP6 * EEP7 * CNTRL * CNTRL 1 * 
8CNTRL2 *CNTRl 3*CNTRL4 *CNTRL5*CNTRL6*KNTRL7*CNTRL8*CNTRL9*CNTR10* 
9KNTR1 1 *CNTR12*CNTR13*CNTR14*CNTR1S*CNTR16*CNTR17*CNTR18*CNTR19* 
XCNTR54* ERASE* ERAS I * ERA S2 * ERA S3 * ERAS4 * ERASE * ERAS6 * ERAS? * 
XERAS8*ERAS9*ERAS10*ERAS1 1 *ERAS12*ERAS13*ERAS14*ERAS1S*ERAS16 
COMMONERAS17*ERAS18*ERAS19*ERAS20*ERAS21 * ERAS22 * ER AS23 * 
1ERAS24*ERAS2S*ERAS26*ERAS27,ERAS28*ERAS29*ERAS30*ERAS31 * 
2ERAS32*ERAS33*ERAS34 *ERAS3E*DEC*DEC1 * DEC2 * DECS * DEC4 * 
3DEC5*DEC6*DEC7*DEC8*DEC9*DEC10*DEC1 1 *<ON*CONl *PRT* 

4PRT1 *PRT2*PRT3*PRT4 * PR TE*PRT6* PR T7*PRT8*PRT9 * DEG* 
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5 TEMP, TEMPI , TEMP2 ♦ TEMP3 , TEMP4 . X I , YI ,THI ,WI , SMU I , 

6CMUI , STH I , CTH I ,S I ,K£R10,EX1 ,EX2,C1 , C2 , C3 , C4 , C5 , C6 , 
7C7,C8,C9,C10,C1 1 , C 1 2 , C 1 3 , C 1 4 , C 1 5 , C 1 6 , C 1 7 , C 1 B , C 1 9 , C20 , 

8C21 ,C22,C23,C24 ,C25, C26, C27, C28,C29,C30,SAV,SAV1 , 
9SAV2,SAV3,SAV4 . S A V5 , SA V6 , S A V7 , S A V8 , S A V9 , SC 1 , THCPR , ENEWST , 
XENEWCT,EKP,EKP1 ,EKP2 , EKP3 , EKP4 , EKP5 , ENEWTH , T HC 1 lENEWM, 

XAA ,KAC,NRET ,SW2 ,EN ( 9 ) , PRNOPT ,HOLD , J ,K ,L , STOR ( 100 0 > 

ERASl =CNTRL3+CNTRL3 
ERAS2=CNTRL3+DEC 1 
ERAS3=CNTRL3-DEC 1 
ERAS4=ERAS3/ERAS1 
1F(H0LD)3,2,2 

3 ERAS5=EKP5*ERASE 
GOTOl 

2 ERAS5=EEP5*ERASE 

1 ERAS4= ( (ERAS5-ERAS4 ) /ERAS2 ) *ERAS I 
ERAS 1 =ERAS3/DEC2 
ERAS5=DEC1 /ERAS5 

ERASl = ( (ERAS5+ERAS1 1/ERAS2 )*DEC2 

ERASl =ALOG (ERASl )*CNTRL3 

SC= ( ( ALOG (ERAS4 1+ERASl ) /ERAS3 ) *CNTRL5 

RETURN 

END 

SIBFTC SPACE REF 

SUBROUTINE SPACE 

COMMONXA , YA , THA , WA , SMUA, CMUA , STHA , CTHA , SA , Xd , YB , THB , 

1 WB,SMUB,CMUB,STHB,CTHB,Sb,EPB,TEPB,XC, YC,THC,WC, 
2SMUC,CMUC,STHC,CTHC,SC , EPC , TEPC , EMUA , EMUB , EMUC , THAC , 

3WAC ,SMUAC , YAC ,EMUAC , S TH AC , EL AC , T HBC , WBC , SMUBC , YdC , EMUBC , 

4STHBC , EPBC , TEPBC ,EM8C,TANMA, TANMB ,EMC , PC , CPC , EM , THCP , 

5C0S1 ,TAN1 ,COS2,TAN2,KOUNT, TAPE, ALIGHT, BEGS, SVCON,SVCONl , 
6SVC0N2, SVC0N3, SVC0N4 , SVC0N5 , SVC0N6 , S VC0N7 , S VC0N8 , S VC0N9 , SVCO 1 0 , 
7SVC01 1 , EEP,EEP1 , EEP2 , EEP3 , EEP4 , EEP5 , EEP6 , EEP7 , CNTRL , CNTRL 1 , 
8CNTRL2,CNTRL3,CNTRL4,CNTRL5,CNTRL6,I<NTRL7,CNTRL8,CNTRL9,CNTR10, 
9KNTR1 1 ,CNTR12,CNTR13,CNTR14,CNTR15,CNTR16,CNTR17,CNTR18,CNTR19, 
XCNTR54, ERASE, ERASl , ERAS2 , ER A S3 , ERAS4 , ER A S5 , ER AS6 , ERAS7 , 

XERASa ,ERAS9, ERASl 0 , ERA S 1 1 , ER AS 1 2 , ERAS 1 3 , ERAS 1 4 , ER AS 1 b , ER AS 1 6 
C0MM0NERAS17,ERAS18 , ER AS 1 9 , ERAS20 , ERAS2 1 , ER A S22 , ER AS23 , 

1 ERAS24 , ERASES , ERAS26 ,ERAS27 , ERAS 28, ERAS29 , ERAS30 , ERAS31 , 

2ERAS32 ,£RAS33,ERAS34 , ERA S35 , DEC , DEC 1 ,DEC2 , DEC3 ,DEC4, 
3DEC5,DFC6,DEC7,DEC8,DEC9,DEC10,DEC1 1 ,KON,CONl ,PRT, 

4PRT1 ,PRT2,PRT3 ,PRT4 ,PRT5,PRT6,PRT7,PRTb,PRT9 ,DEG, 

5TEMP , TENJPl , TEMP2 , TEMPS , TEMP4 ,X I , Y I , THI , W I , SMU I , 
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6CMUI .STHI ,CTHI ,S1 .KER10«EX1 .EX2»C1 , C2 » C3 » C4 « C5 » C6 » 
7C7»Ca»C'5*C10.Cll*C12.C13.C14»C15«C16«C17*Cia»C19*C20i 
8C2 1 , C22 1 C23 , C24 tC25» C26 . C27 » C28 i C29.C30 » SAV» SAVl » 
9SAV2»SAV3*SAV4iSAV5*SAV6«SAV7»SAV8»SAV9iSCl » JHCPR i ENEWST » 
XENEWCT,EKP*EKP1 , EKP2 . EKP3 , EKP4 . EKP5 » ENEWTH , THC 1 lENEWM, 

XAA ,KAC » NRET , SW2 »EN { 9 ) » PRNOPT » HOLD » J « l< iL » STOP (1000) 

WR I TE ( 6 . 1 00 ) 
lOO FORMAT (IH // ) 

RETURN 

END 

ilEFTC GEN REF 

SUBROUT I NEGEN 

COMMONXA, YA,THA*WA,SMUA,CMUA . STHA , CTHA , SA , XB , YB » THB » 

1 WB .SMUa , CMUB » STHB i CTHB , So , EPB * TEPB » XC « YC » THC » WC » 
2SMUC*CMUC.STHC.CTHC»SC',FPc.TEPC*EMUA»ElviUB.EMUCiTHAC. 

3 WAG * SMUAC » YAC »EMUAC * STH A C , EL AC » T H8C » WBC * SMUBC < YBC »EMUBC i 
4STH6C»EPBC»TEPEC,EMi3C<TANMA,TANMB,EMC»PC*CPCiEM»THCP» 

BCOSl *TAN1 ,C0S2.TAN2»K0 UNT,TAPE»ALIGHTiBEGS , SVCON » SVCONl i 
6SVC0N2 * SVC0N3 » SVC0N4 , S VCONB » SVC0N6 » SVC0N7 » SVC0N8 , SVC0N9 . SVCO 1 0 ♦ 
7SVC0 1 1 » EE P 1 EE P 1 » EEP2 * EEP3 » EEP4 » EEP5 * EEP6 , EEP7 » CNTRL » CNTRL 1 « 
8CNTRL2*CNTRL3»CNTRL4 * CNTRL B . CNTRL6 » KNTRL7 , CNTRLB , CNTRL9 , CNTR 1 0 ♦ 
9KNTR1 1 ,CNTR12 iCNTR13.CNTR14iCNTR15.CNTR16.CNTR17iCNTR 18,CNTR19» 
XCNTR54, ERASE* ERA SI , ERA S2 * ERA S3 * ERA S4 * ERA S5 , ERAS6 « ERAST . 
XERAS8.ERAS9.ERAS1 OtERASl 1 , ERAS 12 * ER AS 1 3 * ERAS 1 4 * ER AS 1 5 . ERAS 1 6 
COMMONER AS 1 7 * ERAS 1 8 , ER AS 1 9 i ERAS2 0 * ERAS2 1 * ERAS22 * ERAS23 » 

1 ERAS24 ,ERAS2B*ERAS26 * ER A S2 7 * ER AS28 , ERAS29 * ER AS30 * ERAS3 1 * 

2ERAS32 , ERAS33 *ERAS34 .ERAS 35 * DEC . DEC 1 * DEC2 * DEC3, DEC4 * 
3nEC5,CFC6,DFC7,DFC8,DEC9.DEC10*DFCl 1 *K0N*C0N1 *PRT* 

4PRT1 *PRT2 *PRT3 ,PRT4 * PRT5 * PRT6 * PRT7 * PRT8 * PRT9 , DEG * 

5 TEMP .TEMPI .TEMP2.TEMP3.TEMP4.XI .YI .THI .WI .SMUI . 

6CMU1 .STHI ,CTHI .SI .KERIO.EXI .EX2.C1 .C2.C3.C4.C5.C6. 
7C7.C8.C9.C10.C11.C12.C13.C14.C15.C16.C17.C18.C19.C20. 

8C2 1 . C22 . C2 3 . C2 4 . C2 5 . C26. C2 7. C28 . C29 . C30 . SAV . SAVl . 
9SAV2.SAV3.SAV4.SAV5.SAV6.SAV7.SAV8.SAV9.SC1 . THCPR . ENEWST . 
XENEWCT .EKP.EKPl . EKP2 . EKP3 . EKP4 . E KPS . ENEWTH . TH C 1 .ENEWM. 

XAA .KAC . NRET . SW2 .EN ( 9 ) . PRNOPT .HOLD . J .:< .L . STOR (1000) 

COMMON LTEST 

THAC=THA 

THCP=0 

THBC=THB 

KOUNT=0 

WBC = WB 

WAC = WA 

TANMA=SMUA/CMUA 
EMUA=ATAN ( TANMA ) 
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EMUAC=EMUA 
TANMB=SMUB/CMUB 
EMUB = ATAN ( TANMB ) 

em'ubc = emub 

19 ERASE=THAC+EMUAC 
COSl =COS (ERASE ) 

TANl =S I N (ERASE ) /COS 1 
erase =THBC-EMUBC 
C0S2 = C0S (ERASE ) 

TAN2=S IN (ERASE )/COS2 
ERASE=TAN1-TAN2 
ERASl =TAN1*XA 

XC= ( -TAN2*XB+YB+ERAS 1-YA ) /ERASE 

ERAS4=XC-XA 

YC=TAN1*ERAS4+YA 

YAC= ( YC+YA 1/DEC2 

YBC= ( YB+YC ) /DEC2 

ERASE=COS (EMUAC ) 

SMUAC = S IN (EMUAC ) 

TANMA=SMUAC/ERASE 
ERASE = COS (EMUBC ) 

SMUBC = S IN (EMUBC ) 

TANMB=SMUBC/ERASE 
CALLSL I TET (2 »LTEST ) 

GOTO ( 1 » 3 ) ,LTEST 
1 CALLSLITE(2) 

I F (CNTRL9 )4 ,3,4 

4 ERASF=ERAS4*SMUAC 
ERASF=COS 1 /ERASE 
ERAS 1 =XC-XB 
FRASl =SMUBC*ERAS1 
ERASl =COS2/ERASl 
ERAS3=ERASE+ERAS1 
ERASE=ERASE*SA 

SC=(ERAS1*SB+ERASE)/ERAS3 
G0T05 
3 SC=SB 

5 ERASE=WAC*TANMA 
ERASl =WBC*TANMB 
ERAS5=CNTRL3*CNTRL5 
ERAS3=SC-SA 

ERAS2=( (ERAS3*SMUAC)/ERAS5 )*SMUAC 
ERAS3=SC-SS 
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ERAS3= ( (ERAS3*SMUBC ) /ERASE ) *SMUBC 

ERAS6=TANMA*THA 

ERAS7=TANMB*THB 

IF (CNTRL8 )7«6*7 

7 STHAC=S IN ( THAC ) 

STHBC=SIN(THBC ) 

ELAC=( (SMUAC*TANMA )/C0S1 )*STHAC 
EMBC= ( (SMUBC*TANiMB )/C0S2)*STHBC 
IF (YAC + YBC )9*8.9 
9 I F ( YAC ) 1 1 , 1 0 * 1 1 
1 1 I F ( YBC ) 1 3 » 1 2 » 1 3 

12 ERAS8=FRAS 1 +ERAS1 +ERASE 

ERAS9= (ERAS7+ERAS7-ERAS3 )*WBC+WB 

ERAS4= ( ( XC-XA )/YAC ) *ELAC 

ERASl 0=-ERAS2-ERAS6+ERAS4 

THC= (- (ERASl 0*WAC + WA I+ERAS9 ) /ERAS8 

WC= (THC*TANMA + ERAS1 0 )*WAC+WA 

G0T014 

10 ERAS8.= ERASE + ERASE + ERAS 1 
ERAS5= ( (XC-XB )/YBC )*EMBC 
ERAS9= ( ERAS7+ERAS5-ERAS3 )*WBC+W8 
ERASl 0=-ERAS6-ERAS6-ERAS2 
THC= (-(ERAS10*WAC+WA ) +ER AS9 ) /ERAS8 
ERAS8=THC*TANMA 

WC= (FRAS8+ERASS+ERAS10 )*WAC+WA 
GOTOl 4 

8 ERAS8=ERASE+ERAS1+ERASE+ERAS1 
ERAS9= { ERAS7+ERAS7-ERAS3 )*WBC+WB 
ERASl 0=-ERAS6-ERAS6-ERAS2 

THC= ( - ( ERAS 1 0*WAC + WA )+ERAS9 ) /ERAS8 
ERAS8=THC*TANMA 

WC= (ERAS8+ERAS8+ERAS10 )*WAC+WA 
G0T014 
6 ERAS4=0 
ERAS5=0 
G0T015 

13 ERAS4= ( (XC-XA )/YAC )*ELAC 
ERAS5= ( (XC-XB )/YBC )*EMBC 

15 ERAS8=ERASE+ERAS1 

ERAS9= (ERAS7+ERAS5-ERAS3 )*WBC+WB 
ERAS 1 0 = ERAS4-ERAS6-ERAS2 
THC= (- ( ERASl 0*WAC+WA )+ERAS9 ) /ERAS8 
WC= (THC*TANMA+ERAS1 0 )*WAC+WA 

14 CALLMUSR 

IF (DEC )21 «20»21 
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20 K0UNT=K0UNT+K0N 

IF (K0UNT-KNTRL7 > 1 7, 1 6i 16 

16 IF ( ABS (THC-THCP )-CNTRL6 )18» 1 7* 17 

17 THBC= ('THB + THC 1/DEC2 
WBC= ( WB+WC 1/DEC2 
EMUBC= (EMUB+EMUC ) /DECS 
THAC= ( THA+THC 1/DEC2 
WAC= ( WA+WC 1/DEC2 
EMUAC= ( EMUA+EMUC ) /DECS 
THCP=THC 

G0T019 

18 STHC=S I N (THC ) 

CTHC=C0S (THC ) 

21 RETURN 
END 

SIBFTC WRTFL REF 

SUBROUT I NEWRTFL 

COMMONXA » YA » THA » WA , SMUA, CMUA *STHA »CTHA« SA, XB « YB « THB » 

1 WB » SMUB 1 CMUB , STHB » CTHB , SB , EPB» TEPB . XC « YC » THC « WC * 
2SMUC«CMUC»STHC*CTHC«SC,EPC»TEPC*EMUA.EMUBiEMUC«THAC« 
3WAC,SMUAC»YAC,CMUAC»STHAC» ELAC *THBC * WBC » SMUBC « YoC »EMUBC « 

4STHBC * EPBC« TEPBC , EMBC« TANMA • TANMB » EMC ♦ PC ♦ CPC » EM » THCP* 

5C0S1 *TAN1 »C0S2.TAN2»K0UNT,TAPE»ALIGHT.BEGS»SVC0N«SVC0N1 t 
6SVC0N2 ♦ SVC0N3 , SVC0N4 « SVC0N5 , SVC0N6 « S VCON7 « SVCONS « SVC0N9, SVCO 1 0 
7SVC01 1 « EEP.EEPl , EEP2 • EEP3 * EEP4 t EEP5 . EEP6 « EEP7 . CNTRL « CNTRL 1 ♦ 
8CNTRL2,CNTRL3»CNTRL4 , CNTRLB , CNTRL6 » KNTRt_7 , CNTRL8 « CNTRL9 ♦ CNTR 1 0 
9KNTR1 1 ,CNTR12»CNTR13,CNTR14,CNTR16,CNTR16*CNTR17«CNTR18,CNTR19 
XCNTR54 .ERASE , ERA SI . ERA S2 . ERA S3 . ERAS4 . ER ASS . ERAS6 . ERAS7 . 

XERAS8 , ERAS9 , ERAS I 0 . ERAS 1 I . ERASl 2 .ERASl 3 .ERAS 1 4 . ER AS 1 5 . ERAS 1 6 
COMMONERASl 7 . ERASl 8 . ERAS 1 9 . ERAS20 .ERAS2 1 . ERAS22 . ERAS23 . 

1 ERAS24 .ERAS25.ERAS26 . ER A S2 7 . ERAS28 . ERAS29 . ERAS30 . ERAS3 1 . 
2ERAS32.ERAS33.ERAS34 .ERAS3S . DEC . DECl .DEC2 . DEC3 .DEC4 . 
3DEC5.DEC6.DEC7.DEC8.DEC9.DEC10.DEC1 1 .KON.CONl .PRT. 

4PRT1 .PRT2 .PRT3 .PRT4 . PR T5 . PRT6 . PRT7 . PRT8 . PRT9 , DEG. 

5TEMP .TEMPI . TEMP2 . TEMPS . TEMP4 .XI . Y 1 . THI . W I . SMU I « 

6CMUI .STHI .CTHI .SI .KERIO.EXI . EX2 . C 1 .C2.C3.C4.C5.C6. 

7C7 .C8 .C9.C10.C1 1 .C12.C13.C14. Cl 5 . C 1 6 . C 1 7 . C 1 8 . C 1 9 . C20 . 

8C2 1 . C22 . C23 . C24 . C25 . C26. C27 . C28. C29 . C30 . SAV.SAVl . 

9SAV2 . SAV3 . SAV4 . SAV5 . SAV6.5AV7.SAV8.SAV9 . SCI . THCPR.ENEWST . 
XENEWCT.EKP.EKPl .EKP2 . EKP3 . EKP4 . EKP5 . ENEWTH . THC 1 .ENEWM. . 

XAA .KAC . NRET . SW2 . EN (9 ) . PRNOPT . HOLD . J . K . L . STOR ( 1000 ) 

COMMON LTEST 
K= 1 
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J= 1 

TAPE=-TAPE 

CALLSL 1 TET ( 2 * LTEST ) 

GOTO ( 1 » 3 ) »LTEST 

3 RETURN 

1 CALLSLITE(2)- 
Al IGHT=DEC1 
RETURN 

END 

SIBFTC RDA ref 

SUbROUT I NERDA 

COMMONXA » YA » THA * WA « SMUA. CMUA 1 STHA * CTHA < SA , X3 ♦ YB t THB » 

1 WB<SMUB«CMUB*STHB«CTHB,s'B,EPB.TEPB,XC»YC, THCiWC* 
2SMUC.CMUC»STHC*CTHC»SC»EPC*TEPC»EMUA»£MUB,EMUC«THAC» 

3WAC * SiMUAC » YAC »EMUAC 4 ST'hAC * ELAC* THBC» WBC » SMUBC ♦ YbC » EMUBC » 

4STHBC » EPBC. TEPBC ,EMBC« TANMA » TANMB , EMC * PC * CPC * EM » THCP» 

5C0S1 .TANl ♦C0S2,TAN2»K0UNT*TAPE«ALIGHT.BEGS*SVC0N»SVC0N1 * 

6SVC0N2 . SVC0N3 » SVC0N4 « SVC0N5 * SVC0N6 i SVC0N7 » SVC0N8 . S VC0N9 . SVCOl 0 ♦ 
7SVC01 1 , EEP,EEP1 . EEP2 « EEP3 * EEP4 » EEP5 * EEP6 . EEP7 ♦ CNTRL » CNTRL 1 * 
8CNTRL2 . CNTRL3«CNTRL4 .CNTRL5 . CNTRL6 . KNTRL7 , CNTRL8 , CNTRL9 » CNTR 1 0 « 
9KNTR1 1 ,CNTR12»CNTR13.CNTRl4,CNTR15»CNTR16»CNTRl'7«CNTRia«CNTRl9* 
XCNTR54 , ERASE ♦ ERAS 1*ERAS2»ERAS3<ERAS4»ERAS5,ERAS6« ERAS7 » 
XERAS8,ERAS9,ERAS10,ERAS1 1 » ERAS 12, ERAS 13. ERAS 14, ERAS 15, ERAS 16 
C0MM0NERAS17,ERAS18,ERAS19,ERAS20,ERAS21 , ERAS22 , ERAS23 , 

1 ERAS24 .ERAS25 .ERAS26 ,ERAS27, ERAS28, ERAS29 ,ERAS30 , ERAS31 , 

2ERAS32 , ERAS33 ,ERAS34 ,ERAS35 , DEC, DEC 1 ,DEC2 , DECS ,DEC4 , 
3DEC5,DEC6,DEC7,DEC8,D£C9,DE-C10,DEC1 1 ,KON,CONl ,PRT, 

4PRT1 ,PRT2,PRT3,PRT4,PRT5,PRT6,PRT7,PRT8,PRT9,DEG, 

5TEMP, TEMPI , TEMP2 , TEMPS , TEMP4 , X I , YI ,THI ,WI , SMU I , 

6CMUI ,STHI ,CTHI ,SI .KERIO.EXI , EX2 ♦ C 1 , C2 , C3 , C4 , C5 , C6 , 

7C7,C8,C9,C1 0,C1 1 , C 1 2 , C 1 3 , C 1 4 , C 1 5 , C 1 6 , C 1 7 , C 1 8 , C 1 9 , C20 , 

8C21 ,C22,C23,C24,C25,C26,C27,C28,C29,C30,SAV,SAV1 , 

9SAV2.SAV3, SAV4 , S A V5 , SA V6 , S AV7 ♦ SAV8 ♦ S A V9 , SC I , THCPR , ENEWST ♦ 
XENEWCT.EKP.EKPl , EKP2 , E<P3 , EKP4 , EKP5 , ENEWTH , THC 1 .ENEWM, 

XAA ,KAC ,NRET , SW2 , EN ( 9 ) .PRNOPT , HOLD , U , K , L , S TOR (1000) 

IF (STOR ( J + 3 ) )2 , 1 , 2 

2 XA=STOR(U) 

YA = STOR (J+1 ) 

THA=STOR (J+2) 

WA=STOR (U+3) 

SMUA=STOR (U+4) 

CMUA=STOR (U+5) 

STHA=STOR (U+6) 

CTHA=STOR (U+7) 

SA=STOR (U+8) 
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J = J+9 

NRET=2 

RETURN 

1 KAC=KNTR1 1+KON-KERl 0 
KER10=0 

NRET= 1 
RETURN 
END 

SIBFTC RDB REF 

SUBROUT I NERDB 

COMMONXA, YA* THA » WA, SMUA, CMUA . STHA , CTHA « SA , X3 . YB« THB » 

1 WB «SMUB ,CMUB «STH0 ,CTHB, SB,EPB«TEPB, XC, YC, THC , WC. 

2SMUC « CMUC » STHC , CTHC » SC , EPC « TEPC . EMUA » EMUB , EMUC . THAC . 

3WAC , SMUAC ,YAC, EMUAC , STHAC , EL AC . THBC f WBC . SMUBC < YdC , EMUBC < 
4STHBC« EPBC * TEPBC t EMBC » T ANM A « TANMB <EMC » PC « CPC . EM « THCP* 

5COS1 «TAN1 , COS2, T AN2. KOUNT, TAPE « ALIGHT. BEGS. SVCON.SVCONl . 
6SVC0N2 . SVCON3. SVCON4 . SVCON5 . SVC0N6 . SVCON7 . SVCON8 . S VCON9 . SVCO 1 0 
7SVC01 1 . EEP.EEPl . EEP2 . EEP3 . EEP4 . EEP5 . EEP6 . EEP7 . CNTRL . CNTRL 1 . 
8CNTRL2.CNTRL3.CNTRL4 . CNTRL5 . CNTRL6 . KNTRL7 . CNTRL8 . CNTRL9 . CNTR 1 0 
9KNTR1 1 .CNTR 12. CNTR 13. CNTR 14. CNTR 15. CNTR 16. CNTR 17. CNTR 18. CNTR 19 
XCNTR54. ERASE. ERAS 1 . ERAS2 . ERA S3 . ERAS4 . ER AS5 . ERAS6 . ERAS7 . 
XERAS8.ERAS9.ERAS10.ERAS1 1 . ER AS 1 2 . ER AS I 3 . ERAS 1 4 . ER AS 1 5 . ERAS 1 6 
COMMONERAS17.ERAS18.FRAS19.ERAS20.ERAS21 . ERAS22 . ERAS23 . 

1 ERAS24 .ERAS25 .ERAS26 .ERAS27. ERAS28.ERAS29.ERAS30 .ERAS31 . 
2ERAS32.ERAS33.ERAS34 . ER AS38 . DEC . DEC 1 . DEC2 . DEC3 . DEC4 . 
3DEC5.DEC6.DEC7.DEC8.DEC9.DEC10.0EC1 1 .KON.CONl .PRT. 

4PRT1 .PRT2 .PRT3 .PRT4 .PRT5.PRT6.PRT7. PRT8.PRT9.DEG. 

5TEMP. TEMPI . TEMP2 . TEMP3 . TEMP4 . X I . YI .THI .WI . SMU I . 

6CMU1 .STHI .CTHI .SI .KERIO.EXI .EX2.C1 .C2.C3.C4.C5.C6. 

7C7.C8.C9. CIO. Cl 1 . C 1 2 . C 1 3 . C 1 4 . C 1 5 . C 1 6 . C 1 7 . C 1 8 . C 1 9 . C20 . 

8C2 1 .C22.C23.C24. C25 . C26 . C27 . C28 . C29 . C30 . S AV . SA V 1 . 

9SAV2 . SAV3 . SAV4 . SAV5 . SAV6 . S AV7 . SAV8 . SAV9 . SC 1 . THCPR . ENEWST . 
XENEWCT.EKP.EKPl . EKP2 . EKP3 . EKP4 . EKP5 . ENE WTH , T HC 1 .ENEWM. 

XAA .KAC .NRET . SW2 .EN (9 ) . PRNOPT . HOLD . J .K .L . STOP ( 1000) 

EPB=EX1 

TEPB=EX2 

KERl 0=KER1 O+KON 

IF(ST0R(J+3) )2. 1 .2 

2 XB=STOR(J) 

YB=STOR( J+1 ) 

THS=STOR ( J+2 ) 

WB=STOR ( J+3 ) 

SMUB = STOR( J+4 ) 



APPENDIX 


CMUB=STOR( J+5 ) 
5THB=ST0R ( J+6 ) 
CTHB=STOR{ J+7 ) 
SB=STOR { J+8 ) 

J = J + 9 
G0T04 

1 KAC=KNTR1 l+KON-KERl 0 
KER10=0 
G0T03 
4 NRET=2 
RETURN 
3 NRET=1 
RETURN 
END 
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Figure 2 - Inlet model with starting flap open. 
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Figure 7.- Variation of total- pressure ratio with mass-flow ratio. 
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Figure 8.- Variation of kinetic-energy efficiency with mass-flow ratio. 
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^'The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere aiid space. The Administration 
shall provide for the widest practicable ajid appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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